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Preface 
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Cheryl  Moreno  and  Shirley  Utz,  for  their  patience  and  under- 
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Abstract 

A military  contingency  planning  model  is  designed, 
tested,  and  evaluated  using  goal  programming  analysis. 

This  prototype  is  designed  to  aid  the  achievement  of  Con- 
tinental United  States  (CONUS)  and  overseas  theater  require- 
ments for  personnel  of  different  categories  over  an  arbi- 
trary length  of  time  subject  to  resource  and  time  con- 
straints. The  two  alternate  approaches  of  the  research 
in  goal  programming  are  the  linear  goal  programming  simplex 
method  and  the  pattern  search  method.  Tests  are  run  on 
somewhat  small  examples — a three- job  skill,  six-time  period 
case  and  six  priority  levels  of  projected  goals  and  a three- 
job  skill,  twenty-time  period  contingency  of  seven  priority 
levels  of  CONUS  and  overseas  theater  goals.  The  pattern 
search  method  (tailored  to  fit  to  the  particular  policy 
goal  structure)  appears  to  provide  the  analyst  with  reason- 
able results  and  computer  resource  conservation,  whereas 
the  simplex  method  provides  for  a mathematical  global 
optimum  solution  at  increased  expense  in  computer  resources. 


Vll 


I.  Introduction 


Background 

Peacetime  planning  for  wartime  requirements  has  been  a 
part  of  the  military  establishment  throughout  history, 
and  this  thesis  focuses  on  that  planning  for  both  peace- 
time and  combat  contingencies.  The  considerations  which 
are  involved  in  this  planning  are  indeed  numerous  and 
tend  to  complicate  the  process.  However,  the  development 
of  mathematical  programming  from  its  beginnings  in  the 
Second  World  War  has  been  effective  in  streamlining  these 
planning  procedures;  moreover,  the  recent  growth  in  the 
computational  power  of  digital  computers  has  made  feasible 
the  increased  application  of  mathematical  programming 
techniques  to  improve  the  planning  process. 

Linear  programming  is  probably  the  most  widely  applied 
and  carefully  studied  technique  of  mathematical  program- 
ming. Its  practical  nature  stems  directly  from  G.  B. 
Dantzig's  development  of  the  simplex  algorithm  in  1947 
which  solves  the  general  class  of  linear  programming  prob- 
lems. Fundamentally,  linear  programming  requires  the 
optimization  of  a single  linear  objective  function  subject 
to  a set  of  linear  resource,  process,  and  policy  constraints. 
The  idea  behind  this  formulation  is  the  determination  of 
how  to  allocate  scarce  resources  among  competing  activities. 


Unfortunately,  the  problem  for  the  contingency  planner 
may  not  be  realistically  reduced  to  such  a formulation.  A 
contingency  planner  often  must  make  a tradeoff  among  several 
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conflicting  measures  of  effectiveness,  and,  as  a result, 
he  normally  is  reluctant  to  determine  a single  criterion 
to  use  as  an  objective  function  for  the  linear  programming 
format.  Instead,  he  determines  several  goals  to  be  achieved 
and  tries  to  satisfy  them  in  the  face  of  both  goal  conflicts 
and  limited  resources.  Ideally,  resolution  comes  when 
the  decision  maker  develops  a priority  scheme  or  ranking 
of  these  goals  based  upon  his  particular  value  system. 

Once  he  is  satisfied  with  this  ordering  of  objectives, 
the  planner  tries  to  fulfill  goals  of  higher  priority  before 
lower-priority  goals;  as  a result,  he  finds  that  some  of 
the  lower-priority  goals  cannot  be  fully  achieved,  and 
tradeoffs  must  be  introduced. 

In  real  world  planning  situations,  however,  this  trade- 
off process  may  not  be  fully  developed  and  explicitly  fol- 
lowed. The  Air  Force,  for  example,  has  developed  its  own 
professional  planning  staffs  to  assure  that  viable  solu- 
tions are  available  to  answer  operational  questions.  How- 
ever, 

because  of  the  limited  functional  and  organizational 
horizons  of  each  staff,  potential  tradeoffs  cannot 
be  recognized — and  the  information  required  to  make 
tradeoffs  is  often  filtered  out  before  it  is  passed 
up  (the  hierarchy  to  coordinating  and  consolidating 
staff  agencies)  [Ref  18:27]. 

A course  of  action  that  a contingency  planner  may 
consider  in  order  to  generate  solutions  to  his  problem  is 
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the  use  of  goal  programme,  g.  In  1961  A.  Charnes  and  W.  W. 
Cooper  formally  introduced  goal  programming  as  an  extension 
to  linear  programming  (Ref  19:16).  Prior  to  this  develop- 
ment, more  than  one  measure  of  effectiveness  involved  in 
the  objective  forced  the  decision  maker  to  reformulate 
his  entire  problem.  All  measures  of  effectiveness  were 
usually  combined  into  a single  measure — utility;  with  this 
, single  measure  and  a familiarity  with  linear  programming, 

f 

the  decision  maker  solved  the  problem  with  linear  program- 
ming methods.  Unfortunately,  transforming  the  measures 
into  utility  is  difficult  and  clumsy;  reality  offers  no 
such  single  objective  for  every  occasion. 

To  accommodate  this  reality,  the  development  of  goal 
programming  stemmed  from  the  foundations  of  linear  program- 
ming. Goal  programming  tries  to  consider  all  objective 
measures  that  a planner  thinks  important  and  allows  a simul- 
taneous solution  to  a system  of  conflicting  goals.  Without 
a prior  reformulation  of  a single  measure,  this  solution 

* procedure  treats  the  tradeoff  process  and  requires  the 
establishment  of  an  ordinal  hierarchy  of  importance  among 
goals.  (The  hierarchy  stipulates  that  lower-order  goals 

* 

* are  considered  only  after  higher-order  goals  are  fulfilled 

to  the  maximum  extent  possible.) 

Purpose  and  Organization 

The  purpose  of  this  thesis  is  to  develop  and  test 
a proposed  model  for  military  contingency  planning.  The 

3 


plan  of  attack  is  to  use  the  techniques  of  goal  program- 
ming in  order  to  develop  and  evaluate  various  planning 
policies.  This  model  will  aid  in  the  meeting  of  Continental 
United  States  (CONUS)  and  combat  theater  requirements  for 
personnel  for  a specific  time  duration,  subject  to  various 
resource  and  time  constraints. 

A general  representation  of  the  model  for  personnel 
planning  is  shown  in  Figure  1-1.  The  model  is  designed 
to  chart  the  flows  of  personnel  resources  within  the  mili- 
tary to  support  contingency  requirements  for  an  overseas 
theater  engagement.  Supplies  of  human  capital  come  from 
the  national  labor  market,  the  military  reserves,  and  the 
quantity  of  personnel  already  in  the  system  (stationed 
in  the  CONUS,  in  the  overseas  theater,  in  basic  training, 
and  in  technical  training  as  students  or  instructors) . 
Planning  for  an  engagement  involves  consideration  of  the 
demand  for  particular  job  skills  needed  for  operations, 
the  duration  of  the  contingency,  and  several  other  factors 
under  and  not  under  the  planner's  control.  An  example 
of  a factor  under  the  planner's  control  is  the  desired 
ratio  of  students  to  instructors,  while  an  example  of  a 
factor  not  normally  under  the  planner's  control  is  the 
size  of  the  national  labor  market.  Human  resources  flow 
into  the  CONUS  force  and  then  to  the  theater  or  to  the 
instructor  pool  as  conditions  warrant.  Return  of  person- 
nel from  the  theater  or  instructor  pools  is  also  possible. 
Attritions  from  training  and  the  theater  are  provided  to 
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Figure  1-1.  Personnel  Contingency  Planning  Model 
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reflect  the  losses  that  would  occur  during  an  engagement. 

The  model  may  consider  several  job  categories  simultane- 
ously, as  well  as  consider  requirements  for  several  time 
periods . 

By  incorporating  goal  programming  procedures  into 
the  model,  the  planner  may  identify  and  distinguish  between 
goals  of  varying  priorities.  As  the  model  is  tested  for 
more  and  more  time  intervals  and  job  categories,  the  number 
of  equations  and  decision  variables  is  expected  to  grow, 
necessitating  the  use  of  an  efficient  computerized  algorithm. 

The  thesis  is  divided  into  six  chapters.  After  the 
introductory  chapter.  Chapter  Two  deals  with  the  fundamental 
elements  involved  in  the  goal  programming  method  and  dis- 
cusses previous  applications.  Chapter  Three  describes 
the  Contingency  Planning  Model  of  this  research  in  detail 
with  emphasis  on  goal  programming  application.  Chapter 
Four  presents  the  goal  programming  simplex  computer  algo- 
rithm, whereas  Chapter  Five  presents  the  pattern  search 
computer  algorithm;  both  chapters  include  model  testing 
procedures.  Chapter  Six  reviews  the  results  of  the  tests 
and  establishes  conclusions  about  the  usefulness  of  the 


model  and  the  method. 


II.  The  Goal  Programming  Method 


Basic  Elements  of  a Goal 
Programming  Model 

The  basic  elements  of  a goal  programming  model  are  a 
set  of  goal  equations  and  an  objective  function  (or  achieve- 
ment function)  which  evaluates  how  well  the  overall  goal 
structure  is  met.  As  the  structure  of  the  goal  equations 
and  the  objective  function  is  similar  to  the  format  of 
their  counterparts  in  linear  programming,  we  discuss  these 
basic  elements  in  the  vocabulary  of  linear  programming 
and  assume  that  the  reader  has  a fundamental  knowledge 
of  linear  programming  (Ref  13:15-97). 

The  goal  equations  are  similar  to  the  constraint  equa- 
tions of  a linear  programming  model  in  that  the  equations 
are  written  as  linear  combinations  of  the  various  decision 
variables.  Unlike  constraint  equations,  the  righthand- 
side  values  of  the  goal  equations  represent  various  goal 
levels  which  the  decision  maker  wants  to  meet.  The  deci- 
sion maker  can  determine  the  extent  to  which  a set  of  deci- 
sion variable  values  satisfies  a particular  goal  level 
by  evaluating  the  goal  equation  with  these  values  and  com- 
paring this  linear  combination  with  the  rig  thand-side 
value . 

Every  goal  equation  also  includes  both  a positive 
and  a negative  deviation  variable  which  act  similarly  to 


slack  and  surplus  variables  in  linear  programming.  Each 
represents  a measure  of  the  discrepancy  between  the  goal 
level  and  the  actual  level  of  goal  achievement  by  these 
particular  decision  variables.  If  the  goal  is  underachieved, 
the  negative  deviation  variable  of  that  corresponding  goal 
equation  will  be  positive  (and  the  corresponding  positive 
deviation  variable  will  be  zero) . If  the  goal  is  over- 
achieved, the  positive  deviation  variable  will  be  positive 
(and  the  corresponding  negative  deviation  variable  will 
be  zero) . If  the  goal  is  achieved  exactly  so  that  the 
linear  combination  of  the  decision  variables  for  the  equa- 
tion equals  the  corresponding  righthand-side , both  devia- 
tion variables  of  the  equation  are  zero. 

The  second  basic  element  of  a goal  programming  model 
is  the  objective  function.  To  adapt  a real  world  situation 
to  the  class  of  goal  programming  problems,  we  must  create 
a hierarchy  of  importance  among  goals.  The  hierarchy  is 
necessary  in  order  that  the  algorithm  consider  lower-order 
goals  once  the  higher-order  goals  have  been  satisfied  to 
the  maximum  possible  extent.  The  objective  function  is 
not  constructed  in  terms  of  decision  variables  as  it  is 
in  linear  programming;  instead,  the  objective  function 
is  described  in  terms  of  deviations  (the  deviation  vari- 
ables mentioned  earlier)  between  goals  and  the  levels  of 
achievement  within  a given  set  of  goal  equations  (Ref  19: 

22) . To  deviate  more  and  more  from  a particular  goal  implies 
the  tendency  not  to  achieve  that  goal  fully;  hence,  the 


decision  maker's  purpose  is  to  minimize  these  deviations 
from  the  specific  goals  of  the  problem. 

The  model's  objective  function  is  written  in  terms 
of  deviation  variables,  and  those  deviations  from  the  high- 
est priority  goals  are  minimized  first.  If  two  or  more 
goals  have  identical  priority  rankings  (the  decision  maker 
finds  these  goals  to  bear  equal  importance  to  the  organiza- 
tion) , the  decision  maker  imposes  a modified  ordering  scheme. 
Each  goal  at  this  priority  level  within  the  objective  func- 
tion is  weighted  by  having  a numerical  quantity  multiplied 
by  each  deviation  variable.  Weighting  these  deviation 
variables  does  not  create  differing  priorities  for  one 
goal  over  another  at  this  level;  the  algorithm  is  con- 
structed to  handle  the  deviations  from  goals  within  a prior- 
ity class  as  commensurable  but  with  different  tradeoff 
weights.  This  technique  of  weighting  along  with  preemptive 
priorities  allows  for  very  fine-tuned  analysis  of  goal 
achievement  (Ref  19:26-30). 

Fundamentals  of  Goal  Programming 

The  concepts  and  solution  techniques  for  goal  program- 
ming are  best  explained  through  some  examples.  The  follow- 
ing examples  illustrate  the  setup  conditions  and  the  graphi- 
cal solution  approaches  as  well  as  the  simplex  approach. 

Our  first  example  involves  the  training  and  alloca- 
tion of  two  skill  categories  of  enlisted  men  by  a technical 
training  unit.  The  decision  maker  is  required  to  establish 


goals  concerning  training  loads  and  the  utilization  of  train 
ing  resources.  The  following  information  and  data  are 
available  for  his  consideration: 

1.  Personnel  requirements  forecast  for  the  CONUS 
force  and  for  the  time  period  of  interest  are  six  thousand 
persons  and  eight  thousand  persons  with  skills  1 and  2, 
respectively . 

2.  Ten  thousand  of  the  individuals  (either  all  with 
skill  1,  all  with  skill  2 , or  a combination  of  both)  is 
the  total  training  capability  available  to  support  the 
CONUS  personnel  quotas.  This  capability  refers  only  to 
regular,  not  overtime,  working  hours  (necessary  day  shift 
conditions) . 

3.  Requirements  for  individuals  with  skill  type  1 
are  deemed  twice  as  urgent  as  those  demands  for  skill 
type  2.  For  the  alignment  of  his  objectives,  the  decision 
maker  believes  that  he  must  give  primary  consideration  to 
maintaining  full  capacity  utilization  during  regular  work- 
ing hours.  Secondly,  he  desires  to  meet  all  training 
requirements.  Finally,  he  wants  to  minimize  any  overtime 
operations  (use  of  crews  on  shifts  other  than  daytime 
shifts) . 

The  above  considerations  are  incorporated  into  the 
goal  programming  problem  which  follows: 


Minimize 

Z = P d“  + P2(2d2 

+ dj 

+ P3dt 

Subject  to 

X1  + x2  + d" 
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= 10000 

X1  + d2 

-d2 

= 6000 

x2  + d3 

"d3 

= 8000 

x . > 0 , dt  > 0 , d . 
i - 3 ~ 1 

> 0 

i = 1 / 2 ; j = 1/2/3 

where 

x^  is  the  number  of  individuals  of  skill  type  i 
is  the  priority  level  associated  with  a set  of 
deviation  variables  in  the  objective  function 

d^  is  the  underachievement  (negative  deviation)  of 
the  total  capability  goal 

d£  and  d^  are  the  underachievement  of  the  predicted 
personnel  requirements  goals  for  skill  types  1 
and  2,  respectively 

d^  is  the  overachievement  (positive  deviation)  of  the 
total  capability  goal 

The  remaining  deviation  variables  d*  and  d^  are  the 
levels  of  overachievement  of  the  predicted  person- 
nel requirements  goals. 

The  choice  variables  are  and  the  training  levels 

of  the  two  skills.  The  first  goal  equation  states  that 
there  is  a regular-working-hours  utilization  capability 
of  ten  thousand  individuals;  the  second  and  third  equations 
express  the  personnel  forecasts  for  both  skills;  each  goal 
equation  includes  variables  which  measure  the  amount  of 


positive  or  negative  deviation  from  the  values  on  the  right- 
hand-side  of  the  equation.  (Positive  deviation  implies 
goal  overachievement;  negative  deviation  implies  goal  under- 
achievement . ) 

The  objective  function  is  constructed  from  the  devia- 
tion variables  relating  to  goals  in  three  priority  levels-- 
Pl'  P2'  anc^  P3*  T*le  hi9hest  priority  goal  is  to  minimize 
d^;  this  is  equivalent  to  stating  that  utilization  capabil- 
ity of  at  least  full  regular  working  hours  is  desired. 

The  second  priority  goal  is  to  minimize  the  quantity 
(26.^  + d^)  and,  thus,  to  effect  maximum  levels  of  personnel 
in  training;  moreover,  the  training  for  skill  type  1 is 
considered  twice  as  important  as  that  for  skill  type  2. 

The  third  priority  goal  is  to  minimize  d+  and,  hence,  keep 
the  amount  of  overtime  operations  minimal. 

The  priorities  of  the  problem  are  defined  ordinally, 
and  the  priority  factors  P^'s  are  not  to  be  interpreted 
as  multiplicative  factors.  The  goal  programming  objective 
function  is  to  be  minimized,  and  these  priority  factors 
which  are  associated  with  the  deviation  variables  at  dif- 
ferent priority  levels  have  the  property  that  the  devia- 
tion variables  with  P^  are  minimized  to  the  best  possible 
extent  prior  to  the  variables  with  Pj+^.  In  this  example 
problem  with  P^,  P and  P^,  d^  is  minimized  to  the  best 
possible  extent  prior  to  the  sum  2d^  + d^;  the  sum  is  mini- 
mized to  the  best  possible  extent  (without  violating  the 
previous  minimization)  prior  to  d^.  This  definition 
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assures  that  each  succeeding  goal  is  optimized  to  the  extent 
that  preceding  goals  already  attained  are  not  violated. 

This  constitutes  a "preemptive  priority  structure." 

The  second  priority  term  of  the  objective  function 
includes  a relative  weight  of  2 assigned  to  d^  and  a rela- 
tive weight  of  1 assigned  to  d^.  This  weight  of  2 enables 
the  algorithm  to  minimize  d ^ more  quickly  than  to  minimize 
d^,  reflecting  the  decision  maker's  ranking  of  importance 
of  skill  types  within  the  second  priority. 

This  example  may  be  solved  by  either  a graphical  method 
or  a simplex  method.  Both  techniques  are  similar  to  their 
linear  programming  counterparts.  The  graphical  method 
(Ref  19:64-92)  will  be  presented  since  it  demonstrates 
how  deviations  and  the  priority  structure  of  goal  program- 
ming functions  are  handled.  The  simplex  method  will  be 
discussed  in  a subsequent  example. 

Figure  2-1  shows  the  goal  equations  with  possible 
deviations;  Figure  2-2  shows  the  feasible  region  for  the 
first  priority  goal.  d^  is  minimized  in  the  shaded  area, 
as  its  minimization  is  of  highest  importance  to  the  deci- 
sion maker.  The  new  feasible  region,  with  the  addition 
of  the  second  priority  goals,  is  shown  in  Figure  2-3.  This 
new  region  does  not  violate  the  attainment  of  the  first 
goal.  In  the  shaded  region  d^,  d^,  and  d^  are  zero.  The 
feasible  region  for  only  the  third  priority  goal  is  the 
triangular  area  enclosed  by  the  coordinate  axes  and  the 
line  + x2  = 10000 » however,  a solution  in  this  region 


would  violate  previously  attained  goals.  Therefore,  d* 
is  minimized  at  the  closest  possible  point  to  the  line 
+ x2  = 10000,  which  does  not  violate  the  goals  already 
met.  The  point  is  (x^Xj)  = (6000,8000)  (see  Figure  2-4). 
The  optimal  solution  is  (x1#x2)  = (6000,8000)  with  d^  = 
d2  = d^  = d2  = d^  = 0 and  d*  = 4000.  Only  the  third  goal 
is  not  completely  attained.  Modification  of  the  priority 
structure  or  the  introduction  of  new  goals  could  result 
in  different  solutions. 

It  is  significant  to  recognize  that,  in  no  case,  were 
both  the  positive  and  negative  deviation  variables  for 
any  constraint  equation  simultaneously  nonzero.  This  situa- 
tion is  not  possible  and  will  never  happen  in  a properly 
formulated  goal  programming  problem  (Ref  19:24). 

Goal  programming  simplex  is  discussed  in  the  second 
example.  This  method  differs  from  its  linear  programming 
counterpart  in  that  the  objective  function  is  always  mini- 
mized, the  deviation  variables  take  the  place  of  slack, 
surplus,  and  artificial  variables,  and  the  preemptive 
priority  structure  is  ensured.  The  sequence  of  required 
actions  in  this  simplex  method  is  explained  as  the  example 
problem  is  solved.  (References  15  and  19  provide  detailed 
discussions  of  goal  programming  simplex.) 

The  second  example  comes  from  Ref  19:97-104. 
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Minimize  Z = P-^d"  + P.,d4  + P3(5d2  + 3d~)  + P4d* 


Subject  to  + x2  + d^ 


= 80 


= 70 


= 45 


X1  + X2 


= 90 


x^  , x2  i d^ , d2  # d^ , d4  / d^  , d4  0 


The  second  and  third  goal  equations  do  not  contain  positive 
deviation  variables,  for  while  underachievement  of  the 
righthand-side  values  is  possible,  overachievement  is  not 
in  this  case.  If  the  positive  deviation  variables  were 
included  in  the  problem  formulation,  they  would  have  to 
be  incorporated  into  the  objective  function  as  priority 
one  goals.  Consequently,  the  original  terms  in  the  objec- 
tive function  would  be  assigned  to  the  next  lower  priority 
levels.  Then  the  first  priority  would  be  to  drive  these 
new  positive  deviations  to  zero  before  concentrating  on 
the  lower-order  goals. 

Before  the  first  simplex  tableau  is  presented,  a few 
relevant  points  are  considered.  First,  the  purpose  of 
the  objective  function  is  to  minimize  the  total  unattained 
goals.  This  minimizing  of  the  associated  deviation  vari- 
ables occurs  by  using  preemptive  priorities  and  weights, 
which  take  the  place  of  the  c^'s  of  linear  programming 
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(the  linear  constants  used  as  coefficients  of  decision 
variables  in  the  linear  objective  function) . 

Second,  since  preemptive  priority  factors  are  ordinal 

rather  than  cardinal  values,  they  are  not  commensurable 

at  different  levels.  As  a result,  the  "z.-c."  criterion 

1 1 

used  in  linear  programming  cannot  be  expressed  by  a single 
row  in  the  goal  programming  simplex  tableau.  Instead, 
the  simplex  criterion  becomes  a matrix  of  size  m by  n with 
m as  the  number  of  preemptive  priority  levels  and  n as 
the  total  number  of  both  choice  and  deviation  variables. 

Third,  since  the  criterion  for  determining  the  enter- 
ing basic  variable  is  now  a matrix  and  not  a row,  a new 
procedure  is  required  to  identify  the  entering  variable. 

As  Pj>>>Pj  + ^,  Pj  always  takes  precedence  over  Pj+^.  T*le 
selection  procedure  for  the  optimum  entering  variable  column, 
therefore,  must  consider  the  level  of  priorities  (Ref  19: 
98-99)  . 

Table  II-l  presents  the  initial  tableau  of  the  example. 
As  in  linear  programming,  the  basic  assumption  in  establish- 
ing this  tableau  is  that  the  initial  solution  occurs  at 
the  origin.  At  this  position,  values  of  all  choice  vari- 
ables are  zero.  This  condition  characterizes  the  maximum 
underachievement  of  the  righthand-side  values,  and,  there- 
fore, only  the  negative  deviation  variables  are  the  initial 
basic  variables.  Since  all  the  negative  deviation  vari- 
ables are  nonzero,  the  corresponding  positive  deviation 
variables  must  be  zero.  The  basic  variables  are  listed  in 
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the  first  four  rows  under  column  V in  the  table  with  their 
corresponding  values  under  column  C.  It  is  always  the 
rule  that,  in  the  initial  tableau,  the  negative  deviation 


variables  (d^)  will  appear  in  the  solution  base  (Ref  19:99) 


Zeroes  are  suppressed  to  make  the  table  more  readable. 


The  simplex  criterion  (z^-Cj)  is  a 4 x 8 matrix  because 


of  the  four  priority  levels  and  eight  variables — two  choice 
plus  six  deviation — in  the  model.  The  solution  procedure 
begins  to  achieve  the  most  important  goal  to  the  fullest 
possible  extent  and  then  considers  the  next  highest  order 
goal.  The  unit  contribution  rate  of  each  nonbasic  variable 
in  achieving  the  most  important  goal  is  the  basis  for 
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selecting  the  optimum  entering  variable  column.  Once  the 
first  goal  is  attained  to  the  maximum  extent  then  the  opti- 
mum column  selection  criterion  is  based  on  the  achievement 
rate  for  the  second  goal.  Hence,  when  the  order  goal 
is  attained,  this  criterion  is  based  on  the  unit  contribu- 
tion  rate  for  the  (j+l)'"  goal.  Preemptive  priority  fac- 
tors are  listed  from  lowest  to  highest  in  the  table  so 
that  the  optimum  column  can  be  identified  easily  at  the 
bottom  of  the  tableau. 

In  linear  programming  minimization,  the  righthand- 
side  value  of  the  (z^-c^)  row  represents  the  total  cost 
of  the  solution.  In  goal  programming,  the  values  of 
in  the  constant  column  represent  the  unattained  portion 
of  each  goal.  For  example,  x^  = x2  = 0 in  the  initial 
solution,  and  therefore,  d2  = 70  and  d^  = 45.  The 
unattained  portion  of  the  priority  three  goal  is  5d2  + 3d^  = 
5(70)  + 3(45)  = 485.  Similarly,  if  one  unit  of  x^  is  added 
to  the  solution,  d2  decreases  by  one  unit  and  5d2  decreases 
by  five  units.  Such  an  action  would  result  in  moving  the 
solution  five  units  closer  towards  the  satisfaction  of 
the  third  priority  goal.  By  adding  one  unit  more  of  x2 , 
the  problem  shows  a decrease  in  d^  by  one  unit  and  effects 
a shift  of  three  units  towards  meeting  the  third  goal. 

These  entries  of  5 and  3 in  the  respective  x^  columns  of 

the  P,  row  of  the  (z.-c.)  matrix  are  calculated  directly 

J J D 

from  rows  2 and  3 of  the  basic  variable  section  of  uhe 

tableau.  (From  row  2,  5P^  under  c.  is  multiplied  by  1 
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under  x1  with  the  result  positioned  under  Xj^  in  the  P3  row; 
likewise,  from  row  3,  3P3  under  c..  is  multiplied  by  1 under 
x2,  and  the  result  is  put  under  x2  in  the  P3  row).  Nothing 
is  subtracted  from  the  c..  row  at  the  top  of  the  tableau,  for 
entries  of  zero  exist  above  x1  and  x2-  The  entries  in  the 
P^  row  of  the  (z^-Cj)  matrix  under  x.^  and  x2  are  found 
similarly. 

Another  point  to  consider  is  the  set  of  column  entries 
for  d^  in  the  (z^-c^)  matrix.  With  the  initial  solution  at 
the  origin,  a unit  increase  in  d^  yields  an  increase  in  d+, 
moving  the  solution  one  unit  away  from  both  the  first 
priority  goal  and  the  fourth  priority  goal.  Therefore,  the 
entries  in  that  column  for  d*  are  reasonable.  Again,  this 
event  will  never  happen  because  d+  and  d~  cannot  be  nonzero 
simultaneously.  The  P^  row  entry  in  the  (z^-Cj)  matrix  is 
found  by  multiplying  -1  under  d+  in  the  first  row  of  the 
basis  matrix  by  the  P^  value  under  c^.  (The  P^  row  entry  is 
found  similarly.)  As  is  the  case  with  initial  basic  vari- 
ables of  linear  programming  simplex,  all  coefficients  of 

the  basic  variables  in  the  z.-c.  matrix  are  zero. 

1 1 

Once  the  initial  tableau  is  constructed,  the  entering 
and  leaving  basic  variables  are  determined.  The  column 
with  the  largest  positive  value  at  the  P^  level  in  (z^-Cj) 
is  selected  as  the  optimum  column.  In  Table  II-l  there 
are  two  identical  entries  in  the  x^  and  x2  columns  at  the 
P^  level.  In  order  to  break  this  tie,  the  next  order 
priority  level  is  checked.  Since  a greater  value  in  the 
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x^  column  exists  at  level  P^  than  that  value  corresponding 
to  X£ » the  x^  column  is  chosen.  As  with  linear  programming 
simplex,  the  programmer  determines  the  departing  row  to 
be  the  one  with  the  minimum  value  when  the  righthand-side 
constants  are  divided  by  the  coefficients  of  the  optimum 
column.  The  coefficient  of  1 is  circled  in  Table  II-l  to 
indicate  that  it  lies  at  the  intersection  of  the  optimum 
column  and  the  departing  row.  The  entrance  of  x^  into 
the  basis  causes  a reduction  in  d^  and  d2  and  an  improve- 
ment in  goal  levels  one  and  three  simultaneously. 

With  the  use  of  the  linear  programming  simplex  algo- 
rithm, the  first  tableau  is  revised  to  yield  the  second 
tableau  in  Table  II-2.  After  the  elements  in  the  basis 
matrix  are  calculated,  the  elements  in  the  (z^-Cj)  matrix 
are  determined  as  previously  explained.  From  Table  II-2 
we  see  that  x^  = 70  and  x2  = 0.  The  unattained  portion 
of  the  first  goal  has  been  reduced  to  10  units,  and 
unattained  portion  of  the  third  goal  to  135  units. 

Since  a positive  entry  remains  in  the  row  of 

(Zj-Cj),  it  is  possible  for  more  improvement  in  goal  1 

to  occur.  x2  becomes  the  entering  variable  and  d^  the 

leaving  variable.  Table  II-3  presents  the  third  stage 

solution.  x^  = 70,  x2  = 10,  and  the  first,  second,  and 

fourth  priority  goals  are  completely  achieved;  moreover, 

since  each  entry  in  the  P.  and  P_  rows  of  (z.-c.)  is  non- 

z 3 D 

positive,  no  further  improvements  in  either  of  these  two 
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goals  is  possible.  Since  there  is  a positive  entry  in 
the  row,  achievement  of  the  third  goal  may  be  improved. 

d*  becomes  the  entering  variable,  and  d^  becomes  the 
leaving  variable.  Our  present  objective  is  to  improve 
the  level  of  satisfying  this  third  goal,  but  accomplishing 
this  objective  must  not  degrade  the  existing  achievements 
of  higher-order  goals.  Since  no  negative  entries  lie  in 
the  and  P2  rows  of  the  (z^-c^)  matrix  under  d+,  the 
entry  of  d*  into  the  basis  will  not  pose  any  degradation. 

The  fourth-stage  tableau  appears  in  Table  II-4.  By 
examining  the  (z^-c^)  matrix  we  see  that  the  third  priority 
goal  achievement  can  be  improved  but  only  at  the  expense 
of  worsening  the  achievement  level  of  the  second  priority 
goal.  Similarly,  the  fourth  priority  goal  achievement 
can  be  improved  but  only  by  effecting  a degradation  in 
the  third  priority  goal.  Thus,  no  further  achievement 
is  possible;  the  optimal  solution  is  x^  = 70,  X2  = 20, 
d^  = 25,  d+  = 10,  and  other  goal  deviations  are  zero.  The 
first  two  goals  are  completely  attained,  while  underachieve- 
ment of  the  third  and  fourth  goals  by  25  units  and  10  units, 
respectively,  are  the  closest  possible  attainments  we  can 
expect. 

As  with  linear  programming,  procedures  have  been 
developed  to  handle  complications  in  goal  programming, 
such  as  negative  righthand-side  values,  ties  in  entering 
or  departing  variables  and  alternate  optimal  solutions. 
References  19  and  15  detail  such  procedures  as  well  as 
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Fourth  Tableau 


post-optimality  analysis  and  duality.  Ignizio  (15)  dis- 
cusses the  topics  of  integer  goal  programming  and  non- 
linear goal  programming  but  neither  is  incorporated  into 
the  model  development  in  this  thesis. 

Applications 

Goal  programming  has  been  applied  to  several  decision 
areas  including  problems  of  resource  allocation,  sched- 
uling, and  policy  analysis  (Ref  15:2).  While  the  number 
of  published  applications  is  small,  the  recent  evolution 
of  computer  codes  for  solving  practical  problems  has  created 
a rapid  growth  in  applications. 
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One  of  the  earliest  applications  of  goal  programming 
was  formulated  by  A.  Charnes,  W.  W.  Cooper,  and  R.  J.  Niehaus 
in  1967  for  civilian  manpower  planning  in  the  U.S.  Navy 
(Ref  4) . Subsequently,  several  modifications  were  made 
to  the  original  model,  resulting  in  the  Recruiting  Require- 
ments Model  in  1971,  which  Niehaus  depicted  as  the  "work- 
horse of  the  operational  modeling  system  [Ref  21]."  Vari- 
ations of  this  model  are  in  use  today  for  manpower  planning 
by  the  Navy.  A discussion  of  the  Recruiting  Requirements 
Model  is  provided  in  Appendix  A. 


III.  Contingency  Planning  Model 


Introduction 

This  chapter  describes  in  detail  the  Contingency  Plan- 
ning Model  of  this  research  with  emphasis  on  use  of  the 
goal  programming  solution  method.  The  model  is  developed 
under  a scenario  which  requires  a buildup  and  sustainment 
of  forces  in  an  overseas  theater  to  meet  combat  requirements 
over  an  extended  length  of  time.  Sources  of  personnel 
include  the  civilian  labor  pool,  the  reserves,  the  initial 
training  pipeline,  the  CONUS  force,  and  personnel  initially 
in  the  overseas  theater.  Within  the  model,  provisions 
are  made  for  basic  training,  technical  training,  instructor 
requirements  for  technical  training,  maintenance  of  a 
required  CONUS  force  level,  and  buildup  and  sustainment 
of  a desired  overseas  theater  force.  Consideration  is 
also  given  to  initiation  of  the  military  draft  system, 
the  limitations  on  the  availability  of  civilian  labor  and 
the  reserves,  and  various  attritions  in  the  system.  Several 
job  categories  are  considered  simultaneously  over  several 
time  periods. 

Model  Development 

The  development  of  the  model  begins  with  a considera- 
tion of  its  various  components  describing  levels,  flows, 
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and  restrictions.  The  discussion  which  follows  is  keyed 
to  the  model  network  diagram  in  Figure  3-1. 

The  symbols  used  in  the  network  diagram  and  in  the 
model  equation  formulations  which  follow  are  described 
below. 

Subscripts:  i - Skill  type  indicator 

t - Time  period  indicator.  For  stock 
type  variables,  t designates  a spe- 
cific point  in  time.  For  flow  type 
variables,  t indicates  the  time  inter- 
val (t-l,t) . t will  assume  only 
integer  values. 

- Total  number  of  skill  levels 

- Total  number  of  time  periods  (or  planning  horizon) 

- Number  of  people  in  the  overseas  theater 

- Required  number  of  people  in  the  overseas  theater 

- Flow  rate  of  people  from  the  CONUS  to  the  overseas 
theater 

- Flow  rate  of  people  from  the  overseas  theater  to 
CONUS 

- Combat  attrition  rate  (fraction  of  on-hand  per- 
sonnel lost  in  one  time  period) 

- Number  of  people  in  the  CONUS 

- Required  number  of  people  in  the  CONUS 

- Flow  rate  of  people  from  the  CONUS  to  the  instruc- 
tor pool 


n 

h 

T 

T 

x 

c 

a 

K 

K 

V 


28 


- Ira. > 


wtmmm 


u - Flow  rate  of  people  from  the  instructor  pool  to 
the  CONUS 

y - Flow  rate  of  people  from  the  reserves  to  CONUS 

z - Flow  rate  of  people  from  technical  training  to 

the  CONUS 

R - Number  of  people  in  the  reserves 
R - Desired  residual  reserve  force  size 
Y - Activation  limits  on  reserve  forces 
A - Length  of  technical  training 
X - Attrition  rate  in  technical  training 
w - Flow  rate  of  people  into  technical  training 

b - Flow  rate  of  people  out  of  basic  training 

e - Flow  rate  of  people  into  basic  training 

3 - Survival  rate  in  basic  training 
t - Length  of  basic  training 
E - Desired  trainee  per  instructor  ratio 

L - Number  of  people  in  recruit  supply  pool 

I - Number  of  people  in  the  instructor  pool 

p - Maximum  fraction  of  instructors  in  the  instructor 

pool  permitted  to  return  to  the  CONUS  force  during 
the  next  time  period 

A basic  assumption  is  that  all  transit  times  are  less  than 
one  time  interval. 

Goal  Equation  Formulation 

The  system  of  goal  equations  for  the  model  is  now 
presented.  Each  equation  set  will  be  fully  developed  and 
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discussed  prior  to  a presentation  of  the  next  equation  set. 
Finally,  the  entire  system  of  goal  equations  will  be  sum- 
marized . 


'■ 

i 


People  Available  in  the  Overseas  Theater  We  want 

the  people  available  in  the  overseas  theater,  T.  , to 

1 , t 

equal  the  requirements  specified  by  T.  .in  each  skill 

1 , t 

type  i and  time  period  t.  This  is  formulated 


T . . + d.  . 
i,t  1 , i , t 


d 


+ 

l,i,t 


T 


i/t 


(3.1) 


However,  the  number  of  people  currently  in  the  theater 
force  equals  the  number  available  at  the  previous  time  plus 
inflows  minus  outflows. 


a.  , is  the  combat  attrition  rate  in  the  time  interval 
i/t 

(t-l,t).  Thus,  we  have 


(1  - a . . ) T . , . + x . . - c . . + d.  . . - d.  . . = T . . 

i,t  i,t-l  i,t  i,t  l,i, t l,i, t i,t 


Vi;  t = 2,3,... ,h 


(3.3) 


An  initial  value  for  T.  n is  required  for  t = 1,  and  the 

1 , u 

resulting  goal  equation  is  expressed  as 


x.  i - c . . + d.  . . - d.  . . = T. 
i,l  i/1  1/i/l  1/i/l  i/t 


- Vo'1'  ai,l» 


Vi 


(3.4) 


In  equation  (3.2)  the  number  of  people  currently  over- 
seas is  expressed  in  terms  of  the  number  of  people  at  the 
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last  point  in  time  plus  subsequent  inflows  minus  outflows. 
Likewise,  the  number  of  people  currently  overseas  can  be 
expressed  in  terms  of  the  number  of  people  initially  over- 
seas plus  all  subsequent  inflows  minus  outflows.  This  is 

t 

demonstrated  below  for  skill  type  i at  t = 3.  Using  equa- 
tion (3.2)  we  find 
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In  general,  we  can  say  for  time  t ^ 4 

T.  . = T.  „ n1:  (1  - a . .) 
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Therefore,  the  first  set  of  goal  equations  is  represented  by 


T + d - d = T. 

i,t  + al , i , t al,i,t  i,t 


(3.7) 


where  T.  is  expressed  in  equations  (3.5)  and  (3.6)  as 

1 # U 

appropriate.  The  desired  overseas  goal  T.  and  the  initial 

overseas  strength  T.  n are  assumed  to  be  known.  If  values 

1 f u 

for  the  combat  attrition  rate  a.  . are  also  assumed,  then 

l , c 

equation  (3.7)  represents  a linear  set  of  goal  equations. 

People  Available  in  the  CONUS  Base ■ We  desire  also 

that  the  number  of  people  in  the  CONUS  base  K.  satisfy 

a requirement  K.  that  will  fulfill  the  CONUS  mission. 

1 / t 


K.  . + d_  . . - d0  . . = K.  . 
l,t  2,x, t 2,i,t  i,t 


(3.8) 


Again,  the  number  of  people  on  station  at  time  t is  the 
number  on  station  at  the  last  point  in  time  plus  subsequent 
net  flows. 


K.  . = K.  . . + z.  . + c.  . + u,  . + y.  . 

i,t  i,t-l  i,t  i,t  i,t  Jx,t 


' Xi,t  ' Vi,t 


(3.9) 


Thus , we  have 


. + c.  . + u.  . + y.  . - x.  . - v.  . 

i , 1 l r 1 i f 1 i r 1 i f 1 i r 1 


+ d2,i,l  " d2 , i , 1 
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1,1  1/0 


Vi;  t = 1 


(3.10) 


and 


Ki,t-1  + zi,t 

+ ci,t  + 

u . . + y . 
i,t  *i,t 

- x . 

l 

Vi,t 

+ d2,i,t 

- dt  • 

2 t i 1 1 

= K. 

l 

, t 

(3.11) 


Vi;  t > 2 


As  in  the  first  goal  equation  set,  the  number  of 
people  currently  in  the  CONUS  can  be  expressed  in  terms  of 
the  number  of  people  initially  in  the  CONUS  plus  all  subse- 
quent net  flows. 


K.  . = K.  + I (z.  . + c.  . + u.  . 

i,t  1,0  o=l'  i,o  i,0  1,0 


+y.  .-x.  .-v.  .) 

*1,0  i/D  i / 0 


(3.12) 


At  the  beginning  of  the  contingency  we  assume  that  all 
normal  personnel  separations  are  terminated  for  the  dura- 
tion of  hostilities.  Therefore,  no  losses  due  to  attrition 
in  the  CONUS  are  assumed. 

The  flow  rate  of  people  from  technical  training  z.  . 

1 f H 

is  determined  by  the  flow  rate  of  people  into  technical 
training  and  the  technical  training  attrition  rate. 


‘i,t 


= (1  - 


Xi)wi,t-A 


(3.13) 


is  the  attrition  rate  from  technical  training  and  A^  is 
the  length  of  the  technical  training  program.  (A^  is  an 
integer  number  of  time  periods.)  Substituting  the  expres- 
sion for  z.  above  into  equation  (3.12),  we  find 

ift 
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+ E.  ,(c.  . + u.  . + y.  . - x.  . - v.  .) 
3=1  1,3  1,3  1,3  1,3  1,3' 


(3.14) 


Therefore,  the  second  set  of  goal  equations  is  represented 


K.  . + d_  . . - d_  . = K . , 

i,t  2,i,t  2,i,t  i,t 


(3.15) 


where  K.  .is  expressed  in  equation  (3.14)  above.  The 

1 9 t 

desired  CONUS  goal  K.  and  the  initial  CONUS  strength 

1 f U 

K.  n are  assumed  to  be  known.  Additionally,  the  flow  rates 
of  people  into  technical  training  prior  to  the  beginning  of 


hostilities  (w  , w , 

l i 


, w.  n)  and  the  technical 
l , u 


training  attrition  rate  X.  are  also  assumed  to  be  known. 

Movement  from  the  Reserve  Forces.  We  want  to  activate 
the  reserves  as  needed  and  assume  that  no  new  people  are 
placed  in  the  reserves  for  the  duration  of  the  contingency. 
Thus,  there  is  an  upper  limit  on  the  reserve  total  drawdown 
The  third  set  of  goal  equations  is  represented  by 


E^  . y . . + d_  . - dt  i = R.  n - R. 
t=l  'i,t  3,i  3,1  i,0  l 


(3.16) 


where  R.  _ is  the  initial  size  of  the  reserve  and  R.  is 

i,0  i 

the  desired  residual  reserve  force  size.  We  assume  that 


both  R.  n and  R.  are  known. 

We  will  also  assume  that  there  is  a limitation  on  the 
number  of  reserves  which  can  be  activated  in  any  given  time 
period.  The  fourth  set  of  goal  equations  is  represented  by 


(3.17) 


yi,t  + d4,i,t 


d 


+ 

4,i,t 


Vi;  t 


where  Y^^  is  the  activation  limits  on  the  reserves  in  skill 
level  i during  the  time  period  (t-l,t) . We  assume  that  Y 

i,t 

is  known. 

Inputs  into  Technical  Training.  The  flow  rate  of 
people  into  technical  training  in  the  various  skill  levels 
is  limited  by  the  flow  rate  of  people  out  of  basic  training. 


w 


i/t 


t-1 


(3.18) 


(One  time  period  delay  is  provided  for  administrative  pur- 
poses and  leave  between  basic  training  and  technical  train- 
ing.) But 


Dt  " pet-T  (3.19) 

where  B is  the  basic  training  survival  rate  and  e is  the 

t 

flow  rate  of  people  into  basic  training,  x is  the  length 
of  basic  training  (an  integral  number  of  time  periods) . 

The  fifth  set  of  goal  equations,  therefore,  is  represented 
by 


Ei=l  wi,t  " 6et-x-l  + d5,t  ~ d5,t  = 0 Vt  <3-20> 

We  assume  that  8 and  the  flow  rates  of  people  into  basic 
training  prior  to  the  beginning  of  hostilities  (e  ,e, 

-T  1 — T 

. . . ,eg)  are  known. 
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Technical  Training  Instructors . We  would  like  to  main- 


tain a fairly  fixed  ratio  between  the  number  of  trainees 
and  the  number  of  instructors.  This  is  represented  by 


w 


i/t 


E . I . . . + d,  . , 

1 i,t-l  6,i,t 


d 


+ 

6,i,t 


0 


(3.21) 


where  is  the  desired  trainee  per  instructor  ratio  for 
skill  type  i. 

The  number  of  instructors  at  time  t equals  the  number 
at  t-1  plus  net  flows. 


i »t 
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i,t 


(3.22) 


Thus,  we  reformulate  equation  (3.21)  in  the  following  man- 
ner: 


w.  . - E.I.  . „ — E . (v . . . - u.  . .) 

i,t  l i,t-2  l i,t-l  i,t-l 


+ d,  . , - d. 

6 , l , t 6 , l , t 


= 0 (3.23) 


Our  sixth  set  of  goal  equations  for  the  first  time  period 
is  represented  by 


w.  . - E.I.  . - E.  (v.  n - u.  _) 

1,1  i 1,-1  l 1,0  1/0 


+ d,  . , - d,  . , 

6,1,1  6,1,1 


= 0 


(3.24) 


We  assume  that  E.,  I.  .,  v.  n,  and  u.  are  all  known. 

l i,-l  i/O'  i,0 

We  can  also  say  that  the  number  of  instructors  at 
time  t equals  the  number  initially  plus  net  flows. 


Based  on  equation  (3.25)  we  can  reformulate  equation  (3.21) 
for  t > 1 in  the  following  manner: 
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i/t 


u . . ) + d, 

l , j 6 , x , t 
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6,i,t 


E.I.  _ 

l 1,0 


(3.26) 


Equation  (3.26)  represents  our  sixth  set  of  goal  equations 

for  t > 1.  We  assume  that  I.  n is  known. 

l , u 

Flow  Into  Basic  Military  Training.  The  flow  rate  of 
people  into  basic  training  is  limited  by  the  supply  of 
recruits  available  at  time  t.  The  seventh  set  of  goal 
equations  is  represented  by 

e.  + dl  . - dl  . = L Vt  (3.27) 

L / / L / | L L 

We  assume  that  the  supply  of  recruits  is  known  and  is 
related  to  the  timing  on  the  activation  of  the  draft  system. 

Flow  of  Instructors  to  CONUS . In  order  to  provide 
stability  to  the  instructor  force,  the  flow  rate  of  instruc- 
tors from  the  instructor  pool  to  the  CONUS  should  not  exceed 
some  predesignated  percentage  (p)  of  instructors  in  the  pool 
at  any  time  t.  This  policy  may  be  expressed  by 
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i , t 


pIi , t-] 
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8 , i , t 
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+ 

u,  i , t 
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(3.28) 


The  eighth  set  of  goal  equations  for  t=l  is  represented  by 


(3.29) 


ui, 1 + d8,i,l  d8,i,l  pIi,0 


By  using  equation  (3.25)  we  can  derive  an  expression 

for  I . ...  in  terms  of  I.  and  net  instructor  flows. 
i,t-l  i,0 


Vt-l  - "i.O  + «vi,J  - “i,j» 


(3.30) 


Combining  equation  (3.30)  with  equation  (3.28),  we  find  the 
eighth  set  of  goal  equations  for  t > 1. 


ui,t  - ezf-i  (vi,j  ‘ ui,j>  + 

- d8,i,t  “ (>1i, o ,3-31» 


Summary  of  Goal  Equations 

The  eight  sets  of  goal  equations  are  summarized  below. 
Terms  containing  only  constants  (involving  initial  condi- 
tions and  assumed  values  for  certain  parameters)  have  been 
combined  with  the  righthand-side  goal  values  of  the  equa- 
tions. This  action  is  necessary  in  order  to  make  the  goal 
equations  conformable  to  a simplex  solution  procedure. 

First  Equation  Set — People  Available  in  the  Overseas 
Theater. 


x.  ,-c.  ,+d.  , - d.  . , 

i , 1 i , 1 1 , i , 1 1 , i , 1 
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dl,i,3  Tl,3  " Ti , 0 (1  ~ ai,l)  (1  " ai,2] 
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Vi;  t = 3 (3.34) 
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Vi;  t > 4 (3.35) 

Second  Equation  Set — People  Available  in  the  CONUS 


Base. 
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(3.37) 


+ d"  . - d,  . . = K.  - K.  - 
2 f i / 1 2 / i / 1 i f ^ 


- <1  - 


Vi;  t > hi 


Third  Equation  Set — Movement  From  the  Reserves ; 
Reserves  Total  Drawdown . 


*i,t  + d3,i  - d3,i  = Ri,0  - *i 


(3.38) 


Fourth  Equation  Set — Movement  From  the  Reserves : 


Reserve  Activation  Limits . 


yi,t  + d4,i,t  “ d4 , i , t " Yi , t ¥l,t 


(3.39) 


Fifth  Equation  Set — Inputs  Into  Technical  Traininc 


Z1}  . w.  . + dc  . - d+ 

i— 1 i , t 5ft  5ft 


"t-T-1 


E?  . w.  . 

1=1  ift 


t < T + 1 (3.40) 


6et-t-l  + d5,t 


- « 


t > x + 1 (3.41) 


Sixth  Equation  Set — Technical  Training  Instructors. 
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i / 1 6 f l f 1 v / 1 / J- 

= E . I . , -E.  (v.  „ - U.  n) 

l 1,-1  l 1,0  i,0 
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Vi;  t > 1 (3.43) 


Seventh  Equation  Set — Flow  into  Basic  Military 
Training. 


et  + d7,t  " d7 , t Lt 


¥t 


(3.44) 


Eighth  Equation  Set — Flow  of  Instructors  to  CONUS . 


u.  . + d_  . . - d_  . . 
x,l  8,1,1  8,i,l 


Pli,0 


Vi;  t = 1 (3.45) 


,t-l 


u.  . - pE  . . (v.  . - u.  .)  + d_ 

i,t  j=l  i,j  i,j  8,i,t 


d8,i,t  pIi , 0 


Vi;  t > 1 (3.46) 


Objective  Function  Formulation 

The  objective  function  is  constructed  from  deviation 
variables  relating  to  goals  in  various  preemptive  priority 
levels.  Because  of  the  complexity  of  the  contingency  plan- 
ning problem,  all  goals  as  specified  in  the  system  of  goal 
equations  normally  cannot  be  fully  satisfied.  Therefore, 
we  must  establish  a hierarchy  of  importance  among  goals 
and  express  this  hierarchy  in  terms  of  a preemptive  priority 
structure . 

Prior  to  establishing  a preemptive  priority  structure 
and  formulating  the  objective  function,  some  discussion 
of  our  contingency  planning  problem  is  warranted.  The 
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purpose  of  our  model  is  to  assist  us  in  planning  to  meet 
overseas  contingency  requirements.  Therefore,  we  should 
assign  high  priorities  to  goals  associated  with  meeting 
required  overseas  personnel  levels.  At  the  same  time, 
however,  we  should  also  consider  satisfying  mission  com- 
mitments for  personnel  in  the  CONUS.  Failure  to  meet  CONUS 
commitments  adequately  could  result  in  a possible  degrada- 
tion of  strategic  preparedness.  Thus,  we  should  also  assign 
high  priorities  to  goals  associated  with  meeting  required 
CONUS  personnel  levels.  Similarly,  a strong  argument  can 
be  made  for  assuring  that  an  adequate  instructor  pool  be 
provided  for  the  large  input  of  trainees  expected  with 
the  anticipated  initiation  of  the  draft  system.  Finally, 
certain  goals  within  the  problem  relate  to  resource  limita- 
tions which  are  normally  outside  of  the  control  of  a deci- 
sion maker.  There  is  a limit  on  the  availability  of  both 
the  civilian  labor  pool  and  the  reserve  force  which  cannot 
be  exceeded. 

While  all  of  the  goals  mentioned  above  are  important, 
it  is  unrealistic  to  assume  that  all  can  be  simultaneously 
satisfied.  It  is  also  unrealistic  to  pick  one  goal  (such 
as  the  achievement  of  overseas  personnel  requirements) 
and  use  that  goal  in  the  formulation  of  a linear  program- 
ming problem.  Thus,  we  will  establish  our  goals  in  a 
hierarchical  order  and  use  a preemptive  priority  structure 
in  preparing  our  objective  function.  The  hierarchy  of 
goals  reflects  the  needs  and  desires  of  a decision  maker 


at  the  time  he  formulates  the  problem.  Two  different  deci- 
sion makers  (or  the  same  decision  maker  formulating  the 
problem  at  different  times)  could  establish  quite  differ- 
ent goal  hierarchies.  A distinct  advantage  of  goal  program- 
ming is  the  fact  that  a decision  maker  may  want  to  compare 
, the  overall  satisfaction  of  goals  under  several  different 

goal  hierarchies.  Thus,  the  hierarchy  of  goals  and  pre- 
emptive priority  structure  which  we  establish  below  repre- 
sents just  one  of  several  reasonable  possibilities. 

First  Priority  Level . Our  first  priority  level  goals 
are  concerned  with  resource  limitations  which  are  normally 
outside  of  the  control  of  the  decision  maker.  These  goals 
have  to  be  achieved,  and,  therefore,  we  need  assign  no 
weighting  to  the  various  deviation  variables  in  the  first 
priority  level.  Since  we  are  concerned  with  not  exceeding 
certain  resource  limitations,  we  utilize  the  positive  devi- 
ation variables  found  in  the  appropriate  resource  goal 
equations.  Specifically,  we  include  positive  deviation 
variables  from  the  third  goal  equation  set  (Reserves  Total 
Drawdown) , fourth  goal  equation  set  (Reserves  Activation 
Limits) , and  seventh  goal  equation  set  (Flow  into  Basic 
Military  Training)  in  the  first  priority  level.  The  flow 
. rate  of  people  into  technical  training  is  determined 

directly  from  the  flow  rate  of  people  into  basic  training 
and  the  attrition  rate  in  basic  training.  Therefore,  we 
include  both  the  positive  and  negative  deviation  variables 


from  the  fifth  goal  equation  set  (Inputs  into  Technical 
Training)  in  the  first  priority  level. 

We  also  include  in  the  first  priority  level  our  con- 
cern for  maintaining  a stable  instructor  pool  by  limiting 
the  percentage  of  instructors  who  may  be  returned  to  the 
CONUS  force  during  any  time  period.  We  include  the  posi- 
tive deviation  variables  from  the  eighth  goal  equation 
set  (Flow  of  Instructors  to  CONUS)  in  the  first  priority 
level. 

Finally,  we  express  a desire  for  a stable  force  struc- 
ture after  the  contingency  is  concluded.  In  the  latter 
periods  of  the  contingency,  the  flow  rate  of  people  into 
basic  training  does  not  occur  in  time  to  help  satisfy  CONUS 
or  overseas  contingency  requirements  due  to  delays  required 
for  basic  and  technical  training.  However,  we  would  like 
to  maintain  a stable  force  structure  after  the  contingency 
in  order  to  help  satisfy  any  unforeseen  future  requirements. 
Therefore,  we  include  negative  deviation  variables  for 
the  appropriate  latter  time  periods  from  the  seventh  goal 
equation  set  (Flow  into  Basic  Military  Training)  in  the 
first  priority  level. 

Second  Priority  Level . Our  second  priority  level 
goals  are  concerned  with  providing  an  adequate  instructor 
force  to  meet  the  needs  of  the  large  numbers  of  trainees 
expected  to  support  our  contingency.  Our  main  concern 
is  that  our  desired  trainee  per  instructor  ratio  not  be 
exceeded.  Thus,  we  include  the  positive  deviation 
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variables  from  the  sixth  goal  equation  set  (Technical  Train- 
ing Instructors)  in  the  second  priority  level.  Addition- 
ally, we  assign  weights  to  those  deviation  variables  for 
the  particular  skill  levels  and  time  periods  that  we  favor 
over  other  skill  levels  and  time  periods.  We  may  use  unit 
weights  if  we  have  no  such  preferences. 

It  is  logical  to  place  the  meeting  of  the  desired 
trainee  per  instructor  ratio  in  a higher  priority  level 
than  the  meeting  of  overseas  and  CONUS  requirements  for 
at  least  two  reasons.  First,  personnel  who  enter  the  mili- 
tary service  to  help  offset  overseas  and  CONUS  requirements 
need  adequate  training.  Second,  the  size  of  an  adequately 
manned  instructor  force  is  small  compared  to  the  overall 
CONUS  and  overseas  force  levels.  Moreover,  each  instructor 
is  capable  of  training  several  trainees.  Therefore,  small 
reductions  in  the  CONUS  force  initially  in  order  to  build 
an  adequate  instructor  pool  will  result  in  very  significant 
gains  in  meeting  overseas  and  CONUS  requirements  later 
in  the  contingency. 

Third  Priority  Level . Our  third  priority  level  goals 
are  concerned  with  meeting  overseas  theater  personnel 
, requirements.  We  would  like  to  satisfy  all  of  our  overseas 

theater  requirements.  Thus,  we  include  the  negative  devi- 
ation variables  from  the  first  goal  equation  set  (People 
Available  in  the  Overseas  Theater)  in  the  third  pric rity 
level.  In  addition,  we  include  a weighting  structure 
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which  reflects  our  preferences  among  the  various  skill 
levels  and  time  periods. 

Fourth  Priority  Level . Our  fourth  priority  level 
goals  are  concerned  with  satisfying  all  of  our  CONUS  per- 
sonnel requirements.  We  include  the  negative  deviation 
variables  from  the  second  goal  equation  set  (People  Avail- 
able in  the  CONUS  Base)  in  the  fourth  priority  level.  We 
include  a weighting  structure  reflecting  our  preferences 
within  the  fourth  priority  level. 

I 

Fifth  Priority  Level.  Our  fifth  priority  level  goals 
are  very  similar  to  our  third  priority  level  goals  except 
that  now  we  desire  not  to  exceed  our  overseas  theater  per- 
sonnel requirements.  Thus,  we  include  the  positive  devia- 
tion variables  from  the  first  goal  equation  set  (People 
Available  in  the  Overseas  Theater)  and  a weighting  struc- 
ture reflecting  our  preferences  in  the  fifth  priority  level. 

Sixth  Priority  Level . Our  final  priority  level  is 
concerned  with  not  exceeding  jour  CONUS  personnel  require- 
ments.  We  include  the  positive  deviation  variables  from 
the  second  goal  equation  set  (People  Available  in  the  CONUS 
Base)  and  a weighting  structure  reflecting  our  preferences 
i in  the  sixth  priority  level. 

Preemptive  Goal  Structure  Summarized 
in  Objective  Function 

We  will  now  incorporate  our  goals  assigned  in  the 
six  preemptive  priority  levels  into  an  objective  function. 
Since  our  various  levels  represent  an  ordinal  ranking  among 
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goals,  we  should  not  combine  objective  function  terms  from 
one  level  with  terms  fi?om  a different  level.  Thus,  it 
is  convenient  to  think  of  the  objective  function  as  a vector 
a which  contains  one  component  a^  for  each  preemptive 
priority  level  (Ref  15:16-17).  With  six  preemptive 
priority  levels,  our  objective  function  has  the  following 
form : 

a = (alfa2,a3,a4,a5,a6)  (3.47) 

Each  term  in  the  objective  function  consists  of  a 
summation  of  applicable  deviation  variables  multiplied 
by  their  representative  weights.  As  a result,  we  can  easily 
identify  the  extent  of  goal  satisfaction  within  each  pri- 
ority level.  (For  example,  a^  = 0 implies  that  all  goals 
th 

at  the  l priority  level  have  been  completely  satisfied.) 
Based  on  our  earlier  formulation  of  six  preemptive  priority 
levels,  we  now  construct  our  objective  function. 

a = (a1,a2,a3,a4,a5,a6)  (3.47) 

where 


al = 

zn 

Li=l 

3,1 

+ E^1 
Lt=l 

E1?  ..  dt 
1=1  4 

/ i » t 

+ 

Ih 

2't=l 
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1=1 
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L 

d7,t 

(3.48) 

a2  = 

t=l 

zn 

Li=l 

W6.i.t 

• 4- 

6 , l , t 

(3.49) 
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and 


a_  = E^,  E1?  - w7  . . d7 

3 t=l  i=l  l,i, t 1,1, t 


(3.50) 


a.  = . E1?  . W_  . ^ d~ 

4 t=l  i=l  2,i,t  2,i,t 


(3.51) 


_ _ vh  n Ta7+  .+ 

5 t=l  i=l  l,i,t  l,i,t 


(3.52) 


. rh  n M+  ,+ 

a6  Zt=l  Ei=l  W2,i,t  d2,i,t 


(3.53) 


tT  = h - min  (A ■ + t);  i = l,...,n 
(tL  is  the  first  time  period  in  which  gains  in 
basic  training  can  no  longer  influence  the  out- 
come of  the  contingency . ) 

W+  . and  W . are  the  weiahts  assigned  to  the 

K/l/U  Kf 1/t 

positive  and  negative  deviation  variables,  respec- 
tively . 


Final  Linear  Goal  Programming  Model 

Our  final  model  consists  of  the  objective  function 
(equations  (3.47)  through  (3.53)),  the  goal  equations 
(equations  (3.32)  through  (3.46)),  and  the  requirement  that 
each  deviation  variable  be  nonnegative  (x.  , c.  .,  v.  , 

1 f t 1/1-  1 f t 

y.  +.  / u.  .,  w.  , e.  > 0,  ¥i,t). 

lft  lft  lft  L. 
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z: 


First  Numerical  Example 


In  the  last  part  of  this  chapter,  we  present  a small 
numerical  example  for  the  model  involving  three  skill  levels 
and  six  time  periods.  Since  we  have  already  developed 
the  general  forms  of  the  goal  equations  and  objective  func- 
tion, we  need  only  provide  the  specific  numerical  values 
for  the  various  goals,  initial  conditions,  objective  func- 
tion weights,  and  parameters  of  the  model.  The  solution 
of  our  example  problem  is  deferred  until  Chapter  Four  in 
which  we  discuss  the  simplex  computer  algorithm. 

Our  example  problem  concerns  the  assignment  of  enlisted 
personnel  in  three  skill  levels  for  six  time  periods.  Each 
time  interval  (t-l,t)  is  one  month,  and,  thus  our  overall 
contingency  lasts  for  six  months.  Our  three  skill  levels 
are  identified  below. 


D 

Enlisted  Personnel  Category 

1 

Mission 

2 

Direct 

3 

Indirect  Mission  Support 

Our  initial  CONUS  force  size  (not  including  trainees 
or  instructors)  is  provided  below  (in  thousands  of  per- 
sonnel) . 


i 

l 

2 

3 

o 

•H 

75 

152 

170 

i 

1 

2 

3 

A. 

1 

1 

2 

1 

X. 

l 

0 

0 

0 

(Cross  leveling  could  be  handled  by  using  a technical  train- 
ing attrition  rate  matrix  A allowing  flows  of  personnel 
among  jobs  and  out  of  the  system.  Minor  modifications 
to  the  second  goal  equation  set  would  be  required.) 

The  initial  technical  training  pipeline  consists  of 
people  who  had  entered  technical  training  prior  to  the 
contingency  but  have  not  yet  entered  the  CONUS  force  by 
the  beginning  of  the  contingency. 


w.  . w.  n 

i,-l  1,0 


Data  related  to  desired  trainee  per  instructor  ratio 
and  initial  instructor  levels  and  flows  is  provided  below 


E . 
l 

Ii,-1 

Ii,0 

Vi,0 

Ui,0 

8 

.209 

.209 

0 

0 

5 

.634 

.634 

0 

0 

10 

.333 

.333  j 

0 

0 

We  desire  that  the  CONUS  force  remain  at  initial 


peacetime  levels  throughout  the  contingency. 


i 

1 

2 

3 

*i,t  t>0 

75 

152 

170 

The  length  of  basic  training  is  assumed  to  be  one 
month  (t=1)  due  to  the  small  size  of  our  problem.  The 
basic  training  survival  rate  is  assumed  to  be  70  percent 
(6=.7).  The  initial  basic  training  pipeline  consists  of 
people  who  had  entered  basic  training  prior  to  the  con- 
tingency, but  have  not  yet  entered  technical  training  by 
the  beginning  of  the  contingency.  (Entry  levels  provide 
replacements  necessary  to  support  a normal  peacetime  force 
turnover  rate  of  once  every  five  years.) 


e-l 

o 

Q) 

11.67 

1 

11.67 

The  length  of  technical  training  is  dependent  upon 
the  skill  level.  We  assume  no  attrition  in  technical  train- 
ing in  order  to  avoid  the  need  for  cross  leveling  of  per- 
sonnel among  skills. 
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We  desire  that  no  more  than  fifty  percent  of  our 
on-hand  instructors  be  returned  to  the  CONUS  force  during 
any  time  period  (p  = .5) . 

We  assume  that  our  initial  reserves  are  equally 
balanced  for  the  three  skill  levels  and  that  no  residual 
reserve  forces  are  to  be  maintained.  (Units  are  thousands 
of  personnel.) 


i 

Ri,0 

R. 

l 

1 

30 

0 

2 

30 

0 

3 

30 

0 

We  also  assume  that  up  to  ten  thousand  available  reserves 
can  be  activated  each  month  in  each  skill  level  (Y^^  = 10 
¥i,t) . 

During  the  first  month  of  the  contingency,  only  peace- 
time levels  of  recruits  are  available.  The  draft  system  is 
initiated  by  the  beginning  of  the  second  month,  the  labor 
availability  increases  accordingly.  (Units  are  thousands 
of  personnel.) 


t 

Lt 

1 

11.67 

>2 

50 
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Our  initial  overseas  theater  force  is  at  its  peace- 


time levels.  (Units  are  thousands  of  personnel.) 


i 

1 

2 

3 

Ti  / 0 

25 

38 

30 

However,  our  contingency  requires  a rapid  growth  in  over- 
seas levels  such  that  our  total  combined  CONUS  and  overseas 
force  level  increases  by  approximately  fifty  percent  by  the 
end  of  the  contingency.  (Units  are  thousands  of  personnel.) 


Overseas  Force  Goals,  T. 

1 / U 


\ 

1 

2 

3 

1 

29 

42 

32 

2 

33 

60 

35 

3 

50 

83 

40 

4 

70 

113 

45 

5 

85 

138 

60 

6 

95 

165 

75 

Finally,  our  monthly  overseas  attrition  rates  depend 
upon  both  the  applicable  time  periods  and  skill  levels. 

We  assume  that  losses  are  fairly  significant  during  the 
initial  months  of  the  contingency. 


Monthly  Combat  Atrrition  Rates,  a. 

1 , t 


1 

.03 

.02 

CN 

O 

• 

2 

.03 

(N 

O 

• 

.02 

3 

.02 

.013 

.01 

D 

.01 

.005 

.001 

5 

.01 

.005 

.001 

6 

.01 

.005 

.001 

We  have  now  established  all  of  the  numerical  values 
needed  to  solve  our  problem  except  for  the  weights  in  the 
objective  function.  Presently,  unit  values  are  assigned 
to  all  weights  to  reflect  no  special  preferences  among 
the  various  skill  levels  and  time  periods.  We  may  modify 
the  weights  later. 

Our  present  example  problem  consists  of  114  decision 
variables  and  105  goal  equations  with  150  deviation  vari- 
able terms  in  the  objective  function.  This  problem  and 
larger  problems  are  solved  in  Chapter  Four. 
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IV.  Development  and  Testing  of  the 
Simplex  Algorithm 


Initial  Steps 

The  simplex  method  of  mathematical  programming  is  an 
effective  analytical  tool  for  determining  the  optimal  value 
of  an  objective  function  given  a set  of  constraints.  In 
a goal  programming  context,  the  objective  function  is  a 
vector  quantity  with  more  than  one  measure  of  effectiveness. 
The  simplex  method  for  linear  goal  programming  iterates 
progressively  on  one  priority  level  of  the  objective  func- 
tion (sometimes  referred  to  as  the  achievement  function) 
at  a time.  This  method  switches  to  the  next  highest  pri- 
ority level  when  improvement  is  no  longer  possible  at  the 
current  level. 

Although  a variety  of  simplex  computer  codes  is  avail- 
able for  linear  programming,  only  three  explicit  codes 
were  found  in  the  literature  for  the  linear  goal  program- 
ming context.  The  code  in  the  Marini  thesis  illustrates 
a multiphase  version  of  linear  programming  in  which  each 
priority  level  is  treated  as  a separate  phase  in  the  exe- 
cution process.  Marini  tries  to  solve  goal  programming 
problems  with  a linear  programming  approach — each  priority 
level  is  solved  as  a distinct  linear  programming  simplex 
problem.  The  automated  nature  of  his  program  allows  every 
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level  to  be  handled  during  one  run  without  manual  initia- 
tion of  each  phase  (Ref  20) . The  other  two  codes  from  Lee 
and  Ignizio  (Refs  19;  15)  are  written  in  a straightforward  man- 
ner and  have  been  tested  for  small  size  samples.  Because 
of  their  ease  of  adaptability,  simplicity,  and  flexibility, 
these  latter  two  codes  appear  attractive  for  the  analysis 
of  the  contingency  planning  model. 

One  criterion  that  we  consider  in  dealing  with  the 
contingency  model  is  the  number  of  goal  equations  used 
to  describe  the  problem.  One  aim  of  our  model  is  to  experi- 
ment with  the  sizes  of  military  forces  in  terms  of  physical 
numbers  of  individuals  and  of  different  job  skills  flowing 
over  several  time  periods  in  the  contingency.  We  want 
to  determine  how  effectively  a personnel  force  configura- 
tion can  support  a given  war  scenario;  we  want  to  determine 
what  size  of  force  is  necessary  to  handle  the  war  scenario 
and  what  size  of  war  (length  of  time,  attrition  factors) 
make  our  capabilities  with  human  resources  a limiting  fac- 
tor. In  order  to  accommodate  more  and  more  job  categories 
and  time  periods,  the  computer  code  needs  the  capability 
to  expand  and  store  more  variables  and  goal  equations. 

Ignizio 's  computer  code,  a modified  simplex  procedure, 
demonstrates  obvious  savings  in  storage  capacity  over  Lee's 
algorithm.  Both  codes  deal  with  the  same  technology  matrix; 
Ignizio  excludes  the  righthand-side  submatrix  of  negative 
deviation  variable  coefficients;  Ignizic's  algorithm  pivots 
on  the  remaining  technology  coefficients  and  places  the 
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departing  basic  variable  of  a particular  row  into  the  column 
of  the  variable  which  enters  the  basis.  This  switching 
of  variables  between  the  pivoting  row  and  column  avoids 
a duplication  of  variables  in  the  basis  and  in  the  column 
headings.  Also,  Ignizio's  code  provides  helpful  subprograms 
which  Lee's  does  not — subprograms  to  generate  alternate 
optimal  solutions  and  to  fix  numbers  as  integers  if  requested 
by  the  user.  (This  latter  subprogram  changes  numerical 
values  which  stray  from  remaining  as  integers  into  integers 
if  the  values  lie  within  a constant  tolerance  of  that  inte- 
ger. Precision  of  variables  and  limited  control  of  roundoff 
error  is  retained.)  Because  of  its  advantages  of  storage 
requirements  and  its  useful  subprograms,  the  computer  code 
in  Ignizio's  textbook  serves  as  the  foundation  of  the  goal 
programming  simplex  for  our  contingency  model. 

Figure  4-1  shows  Ignizio's  original  design;  the  func- 
tional description  of  each  subprogram  follows: 

GMAIN  reads  in  all  input  data  and  controls  the 

execution  phase  of  the  algorithm;  appropriate 
errors  are  printed. 

PLACE  assembles  the  initial  objective  function 

weights  in  the  matrices  of  top  stub  and  left 
stub  values.  In  Ignizio's  context,  these 
matrices  perform  the  same  functions  in  the 
simplex  method  as  does  the  c^  top  row  of  Lee's 
model  (Cj  row  defined  in  Chapter  Two) . 


Li 
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OUTPUT 


Figure  4-1.  Original  Design 
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CINDX 


computes  new  index  rows  of  the  (z^-Cj)  matrix 
and  the  latest  values  of  the  objective  function 
components . 

TEST  determines  the  departing  variable  row  and  the 
entering  variable  column  by  means  of  the  new 
index  row  values. 

PERM  pivots  the  technology  matrix  on  the  element  of 
intersection  of  the  column  and  row  which  are 
determined  in  TEST;  values  are  changed  not 
only  throughout  the  technology  matrix  but  also 
in  the  matrix  of  righthand-side  goal  values 
through  this  permutation. 

FIX  retains  the  precision  of  integer-valued  quanti- 

ties . 

POUT  prints  the  processed  information  of  each 
optimal  solution. 

ALTST  tests  the  final  technology  matrix  for  addi- 
tional alternative  solutions  if  any  exist. 

FCPL  forms  a cutting  plane  construct  for  the  prob- 
lem and  adds  the  necessary  cutting  plane  equa- 
tions to  the  tableau  (technology  matrix) 
whenever  an  integer  solution  is  requested. 

INTST  indicates  whether  or  not  a requested  variable 
is  integer-valued  (Ref  15:227-228,  234-242). 

The  arrows  in  the  figure  represent  the  lines  of  com- 
munication between  the  subprograms  and  between  the  input/ 
output  media  and  the  execution  program.  (The  original 
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design  and  updated  design  are  created  in  the  language  of 
FORTRAN  IV  Extended,  and  the  discussion  which  follows 
assumes  the  reader  to  have  a very  basic  knowledge  of  the 
language . ) 

For  our  contingency  model,  we  decided  early  in  the 
algorithm  development  stages  to  exclude  those  subprograms 
dealing  with  integer-valued  solutions.  Considering  mili- 
tary forces,  the  problem  that  we  are  modelling  deals  with 
thousands  of  personnel,  and  a scaling  of  the  problem  to 
fit  "units"  of  a thousand  men  in  both  input  and  output 
results.  Optimal  values  to  the  nearest  person  are  available 
when  the  problem  is  scaled.  Thus,  FCPL  and  INTST  are 
dropped  from  the  model  as  well  as  the  initial  input  control 
variable  that  specifies  that  an  integer-valued  solution 
be  printed.  Only  program  FIX  is  retained  for  its  ability 
to  help  control  the  effects  of  roundoff  error. 

The  original  design  is  fundamentally  modular,  for 
each  subroutine  performs  one  and  only  one  specific  function 
of  the  algorithm.  It  is  unnecessary,  though,  to  have  the 
main  program  read  the  initial,  unprocessed  data  and  write 
error  diagnostics;  instead,  other  subroutines  are  constructed 
to  handle  this  input/output  processing  and  still  retain 
a modular  structure.  The  addition  of  these  subroutines 
makes  a clear  distinction  between  the  initialization  phase 
and  the  execution  phase  of  the  computer  algorithm.  Further- 
more, the  main  program  now  acts  like  an  executive  module, 
only  controlling  the  program  execution.  If  difficulties 
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in  processing  the  input  and  the  output  arise,  the  main 
program  is  alerted  and  halts  the  process  so  that  the  user 
may  correct  the  error. 

The  initialization  phase  consists  of  the  newly  created 
routines  INIT  and  COEF  and  the  subroutine  PLACE,  whose 
function  remains  unchanged  from  the  original  design.  INIT 
directs  the  reading  of  input  data  so  that,  at  INIT's  con- 
clusion, every  initial  matrix  value  is  in  place,  and  the 
algorithm  is  ready  to  test  for  optimality  of  the  system 
at  the  first  priority  level.  Labels  of  twenty  characters 
maximum  for  every  decision  variable  are  inserted  into  the 
model  prior  to  INIT's  conclusion.  INIT  resembles  a minia- 
ture executive  module  that  does  minimal  input  processing 
and  directs  the  initialization  through  subroutine  calls. 

COEF  generates  initial  coefficients  of  the  goal  equa- 
tions and  inserts  them  into  proper  positions  of  the  tech- 
nology matrix  (referred  to  as  TE  in  the  code) . INIT  is 
responsible  for  reading  the  initial  righthand-side  values, 
but  COEF  is  responsible  for  generating  updated  values 
for  the  righthand-sides  of  the  goal  equations.  The  original 
design  reads  in  the  TE  coefficients  one  row  at  a time  with 
formatted  READ  statements.  Reducing  input /output  time 
in  the  updated  design  is  effected  by  the  calculation  of 

coefficients  with  the  initial  data  such  as  T.  n,  K.  n, 

i,0  i,0 

and  8 which  were  explained  in  the  previous  chapter.  The 
algorithm  requires  a sparse  TE  matrix  (with  few  nonzero 
elements)  where  coefficients  of  earlier  time  periods  of 


the  contingency  help  to  build  the  coefficients  of  subse- 
quent time  periods;  COEF  performs  this  arithmetic. 

The  original  design  reads  objective  function  values 
also  as  combinations  of  the  priority  level,  the  particular 
deviation  variable  subscript,  and  the  weighting  factor, 
with  each  combination  on  a single  card.  Input  processing 
time  is  reduced  even  further  in  the  updated  design.  Our 
model  allows  these  input  values  to  be  read  in  as  four  com- 
binations maximum  of  objective  function  values  on  a single 
card,  and,  hence,  requires  fewer  READ'S  to  the  input  file 
and  one-quarter  as  many  physical  cards. 

The  execution  phase  of  the  algorithm  consists  of  the 
originally  designed  CINDX,  TEST,  PERM,  POUT,  and  ALTST 
(with  slight  modifications  for  efficient  code)  with  an 
updated  version  of  FIX  and  the  addition  of  subroutine 
ERRORS.  ERRORS  is  called  by  the  MAIN  program  only  if  error 
diagnostics  are  needed.  Once  either  ALTST  or  ERRORS  ends, 
the  algorithm  is  complete  with  the  current  problem,  and 
another  is  begun. 

FIX  is  changed  to  permit  a variable  number  L of  deci- 
mal places  of  precision.  L is  an  integer  greater  than 
one.  In  the  original  design  FIX  assured  precision  to  the 
fourth  decimal  place  by  examining  the  input  value  and  con- 
firming whether  or  not  it  lies  in  an  interval  of  .005  from 
an  integer.  With  the  modification,  FIX  examines  the  input 
value  and  determines  if  it  lies  within  a distance  of 
5 x 10‘(L_1)  from  an  integer.  If  so,  the  integer  value 
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is  returned  in  the  function;  otherwise,  the  original  value 
is  returned  unchanged. 


The  design  and  logic  flow  of  the  updated  model  are 
depicted  in  Figure  4-2.  The  algorithm  begins  with  INIT 
directing  the  initialization  phase.  Objective  function 
values  are  in  place  as  are  technology  equation  coefficients 
Control  returns  to  the  MAIN  program  with  the  variable  NROW 
set  to  one.  NROW  serves  as  the  indicator  of  the  current 
priority  level  at  which  the  program  is  operating. 

CINDX  creates  the  index  rows  which  are  tested  later 
with  the  technology  equations  for  the  entering  variables. 
The  criterion  to  select  the  entering  variable  is  aptly 
stated  by  Ignizio:  "Select  the  largest,  positive  I, 

Kf  S 

for  which  there  are  no  negative  valued  index  numbers,  at 

a higher  priority,  in  the  same  column."  (I,  in  the  cita- 

K / S 

tion  is  the  value  in  the  index  matrix  at  the  x priority 

th 

level  (row)  and  the  sc  column  (Ref  15:47).)  The  departing 

variable  is  then  determined  by  finding  the  minimum  value 

of  the  quotients  b^/e^  g where  b^  is  the  righthand-side 

goal  equation  quantity  and  e.  is  the  corresponding  i*"*1 

1 , s 

row  element  of  the  entering  variable  column.  This  deter- 
mination of  variables  of  the  model  is  one  iteration  of 
the  algorithm  at  the  current  priority  level. 

If  the  entering  column  number  is  calculated  to  be 
zero,  the  priority  level  cannot  be  optimized  further,  and 
control  passes  to  the  next  priority  level.  Otherwise, 
permutation  of  the  TE  matrix  followed  by  another  call  to 
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INITIALIZATION  ■ EXECUTIVE  / EXECUTION 


Figure  4-2.  Updated  Model:  Modular  Design  and  Logic  Flow 
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Figure  4- 2 --Continued 


CINDX  prepares  the  process  to  iterate  once  more  to  deter- 
mine the  principal  row  and  column.  This  iterative  process 
continues  for  each  priority  level  until  one  of  two  condi- 
tions is  satisfied:  either  the  column  value  is  zero  or 
the  total  number  of  iterations  exceeds  an  initial  number 
called  IMAX.  (IMAX  for  our  problem  is  200.)  If  the  first 
condition  is  satisfied,  the  priority  level  NROW  is  incre- 
mented by  one,  and  the  process  continues  until  NROW  exceeds 
the  number  of  priority  levels  in  the  objective  function. 
When  NROW  does  exceed  that  number,  the  execution  phase 
becomes  one  of  printing  optimal  solutions  (the  first  opti- 
mal and  then  the  alternate  optimal  solutions)  and  of  stop- 
ping unless  another  problem  is  ready  in  the  data  stream. 
When  the  number  of  iterations  exceeds  IMAX,  subroutine 


ERRORS  has  control  over  execution,  prints  the  error  state- 
ment, and  returns  control  to  MAIN  for  additional  input 
contingency  problems. 

The  design  of  COEF  requires  further  description.  This 
subprogram,  like  the  other  subprograms  of  the  algorithm, 
is  developed  under  the  guideline  of  retaining  variable 
values  instead  of  constant  values  in  the  computer  code 
wherever  possible.  Expansion  to  other  contingency  planning 
problems  is  provided  for  by  these  variables  except  for 
changes  to  the  size  of  the  file  buffers  and  the  array  dimen- 
sions which  are  normally  at  a fixed  level  for  a specific 


run. 


The  contingency  model's  goal  equations  consist  of 
eight  distinct  constraint  sets  as  described  in  the  third 
chapter.  It  is  natural  to  generate  coefficients  of  the 
goal  equations  in  COEF  by  splitting  the  total  number  into 
eight  segments  and  developing  one  constraint  set  at  a time. 
The  number  of  rows  of  equations  generated  for  each  set 
varies;  one  constraint  set  produces  as  many  equations  as 
there  are  job  categories,  a number  which  we  designate  as 
JOBCAT;  two  sets  produce  as  many  as  there  are  time  periods, 
a number  which  we  designate  as  JTIMES;  the  remaining  sets 
require  an  equation  for  each  job  category  and  for  each 
time  period  and,  thus  produce  as  many  as  there  are  in  the 
product  of  JOBCAT  and  JTIMES,  which  we  designate  as  NEW. 

The  initial  TE  matrix  is  established  with  the  decision 
variables  heading  the  lef thand-side  columns  as  nonbasic 
variables  and  the  positive  deviation  variables  heading 
the  righthand-side  columns  as  nonbasic  variables.  (A  vari- 
able heading  a column  means  that  the  coefficients  in  suc- 
cessive rows  of  the  column  belong  to  the  variable  in  the 
goal  equations.)  The  way  that  the  decision  variables  are 
ordered  across  the  top  of  TE  is  arbitrary;  our  setup,  upon 
which  COEF  is  oriented,  is  the  following  block  structure, 
left  to  right:  x's,  c's,  v's,  y's,  u's,  w's,  and  e’s,  where 
these  symbols  are  defined  in  the  third  chapter.  Each  block 
of  variables,  except  for  the  e's  consists  of  a variable 
for  each  job  category  and  each  time  period  and  is  arranged 
thus:  the  variables  of  the  first  time  period  from  the  first 


job  skill  to  the  last  job  skill,  the  variables  of  the  second 
time  period  from  the  first  job  skill  to  the  last,  . .., 
the  variables  of  the  last  time  period  from  the  first  job 
skill  to  the  last.  (This  arrangement  is  left  to  right 
in  the  block.)  The  block  of  e's  consists  of  an  e variable 
for  each  time  period  arranged  from  left  to  right  from  the 
first  to  the  last  time  period;  no  job  skill  shredout  is 
applied  to  these  variables,  the  flows  from  the  national 
labor  supply  into  basic  training.  The  positive  deviation 
variables  are  ordered  left  to  right  in  succession  from 
the  variable  of  the  first  goal  equation  to  the  variable 
of  the  last  goal  equation. 

Zeros  dominate  the  initial  technology  matrix.  The 
righthand-side  submatrix  of  positive  deviation  variables 
is  identified  as  a negative  identity  matrix  of  order  equal 
to  the  number  of  goal  equations.  The  lef thand-side  sub- 
matrix is  neither  upper  nor  lower  triangular;  however, 
a pattern  is  evident.  The  repeated  use  of  coefficients 
of  early  time  periods  in  a given  constraint  set  appear 
to  expand  more  nonzero  coefficients  into  rightward  and 
higher- row-number  positions  of  the  TE  block  for  that  set. 

In  these  constraint  sets  that  produce  NEW  equations,  the 
setup  in  TE  conforms  to  rows  of  equations  for  time  1, 
time  2,  ...,  time  JTIMES  and  within  each  of  these  smaller 
sets  for  a time  period  from  job  category  1 to  job  category 
JOBCAT.  With  this  setup,  each  column  block  of  variables 
included  in  a constraint  set  shows  this  pattern  to  be  lower 


triangular.  (See  Figure  4-3.)  Each  new  subdiagonal  below 


the  block's  main  diagonal  is  positioned  a total  of  JOBCAT 
units  away  from  the  next  diagonal;  also,  a total  of  JTIMES 
diagonals  (including  the  main  diagonal  and  every  sub- 
diagonal) exists  for  each  block  of  decision  variables  per- 
tinent to  the  constraint  set. 

COEF  is  designed  to  produce  this  lower  triangular 
pattern  for  the  constraint  sets  of  NEW  equations  by  two 
DO  loops,  one  nested  inside  the  other  with  the  same 
terminal  statement  for  both  loops.  The  outer  loop  of  the 
form  DO  sn  J=l,  NEW  iterates  on  one  column  of  the  block 
of  variables  at  a time;  J increments  by  one  only  after 
all  coefficients  of  the  constraint  set  are  in  position 
in  that  column.  (sn  in  the  DO  loop  is  the  terminal  state- 
ment of  the  loop.)  The  inner  loop  of  the  form  DO  sn 
I = J,  NEW,  JOBCAT  iterates  on  a particular  starting  row 
and  increments  in  steps  of  JOBCAT,  reflecting  the  relative 
position  of  the  subdiagonals.  The  lower  triangular  pattern 
is  guaranteed  by  the  fact  that  the  inner  loop  never  iterates 
on  a value  less  than  the  column  parameters.  It  is  essen- 
tial to  realize  that  J refers  to  horizontal  movement  in 
the  block  and  I to  the  vertical  movement  in  the  block. 
Increasing  J implies  dealing  with  variables  of  differing 
job  categories  and  increasing  time  periods;  increasing 
I implies  dealing  with  goal  equations  of  later  time  periods. 
For  a further  discussion  of  this  nested  structure,  we 
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Figure  4-3.  One  Block  of  Technology  Coefficients  (Three  Skill,  Four  Period  Example) 
(y . represents  (1-a.  ) and  blanks  represent  zero  coefficient  values  in  the  block) 


advise  the  reader  to  examine  the  comments  of  COEF  of  the 
computer  code  located  in  Appendix  B. 

COEF  has  the  capability  to  modify  any  negative  right- 
hand-side  values  with  the  variable  SWITCH.  Certain  policy 
configurations  may  generate  negative  righthand-side  goal 
equations  values  in  the  second  constraint  set.  As  Ignizio 
warns,  "However,  the  modified  simplex  algorithm  requires 
that  all  righthand-side  values  be  nonnegative  [Ref  15:57]." 
If  the  analyst  or  decision  maker  recognizes  that  this  limita 
tion  has  been  observed  with  the  policy  that  he  has  formu- 
lated, he  sets  a nonzero  value  for  the  input  variable  SWITCH 
A nonzero  SWITCH  forces  the  algorithm  to  multiply  the  right- 
hand-side values  of  the  goal  equation  and  the  lefthand- 
side  TE  values  by  -1.  A zero  SWITCH  bypasses  this  multi- 
plication. Where  negative  righthand-side  values  develop, 
care  is  necessary  in  the  interpretation  of  the  negative 
and  positive  deviation  variables;  the  original  positive 
deviations  of  the  second  goal  set  are  initially  in  the 
basis . 

The  subprograms  of  the  updated  model  conform  to  cer- 
tain standards.  All  subroutine  calls  are  calls  by  refer- 
ence in  which  the  address  of  the  argument  is  passed'  to 
the  subroutine  through  the  argument  list  of  the  CALL  state- 
ment. No  matrix  is  an  argument  of  any  subprogram.  The 
essential  matrices  for  the  algorithm  lie  in  blank  COMMON 
or  labeled  COMMON  blocks  which  provide  the  majority  of 
communications  between  different  program  units. 
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The  design  of  this  source  code  is  also  an  effort  in 
efficiency  with  features  for  saving  time  in  central  process- 
ing as  well  as  in  input/output  processing  and  for  conserving 
the  use  of  central  memory.  Problems  in  running  the  program 
develop  as  the  model  expands  to  accommodate  more  time  periods 
and  job  categories.  More  memory  is  needed;  hence,  ways 
to  reduce  memory  requirements  are  employed.  The  buffer 
size  of  every  file,  inserted  in  the  initial  PROGRAM  state- 
ment, is  normally  set  at  the  default  value  of  2002_  = 1026 

O 

words  in  the  CDC  6600  computer.  This  quantity  is  ample 
for  normal  input/output  processing,  and  optimization  by 
the  FORTRAN  compiler  helps  to  determine  an  optimal  size 
as  well.  A smaller  buffer  size,  however,  decreases  overall 
memory  allocation  but  at  the  expense  of  increased  input/ 
output  processing  time.  The  buffer  size  is  strictly  depen- 
dent on  the  size  of  the  particular  contingency  problem 
to  be  solved  and  on  the  number  of  records  of  information 
processed  at  any  instant  of  time  during  execution.  A trade- 
off in  the  conservation  of  computer  resources  must  be 
resolved,  because  as  the  buffer  size  decreases  central 
memory  allocated  to  the  run,  the  input /output  time  tends 
to  increase.  The  converse  is  likewise  true.  Considering 
our  file  buffers,  we  find  that  a smaller  input  file  buffer 
and  output  file  buffer  of  300g  words  (192  words)  for  small 
contingency  problems  (such  as  a test  case  of  two  skills, 
three  time  periods)  save  space  with  little  increase  in 
input /output  processing  time. 
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The  use  of  ENTRY  statements  (nonstandard  FORTRAN  IV 
features  available  with  the  CDC  6600  machine)  helps  to 
reduce  computer  memory.  As  subroutines  TEST  and  PERM  have 
identical  subroutine  arguments,  a single  subroutine  with 
two  entry  points  satisfactorily  performs  the  same  work 
as  the  individual  routines.  This  structure  results  in 
using  18  fewer  computer  words  with  both  parts  of  the  routine 
having  access  to  the  same  COMMON  block  which  is  declared 
only  once.  Similarly,  subroutine  ERRORS  is  an  alternate 
entry  point  for  the  program  unit  with  entry  point  POUT; 
both  routines  require  one  argument  in  the  CALL  statement. 

Another  consideration  to  reducing  the  number  of  com- 
puter words  is  word  packing.  This  technique  aims  to  use 
more  available  space  in  a computer  word  by  the  storage 
of  more  than  one  value.  The  word-packing  and  bit- 
manipulation  capabilities  with  the  CDC  6600  machine  deal 
only  with  integer  variables.  The  structure  of  the  computer 
word  used  for  integer  values  is  the  complete  60-bit  con- 
figuration without  an  exponent.  For  example,  in  array 

g 

of  positive  integer  values  not  greater  than  512  = 2 , each 
value  requires  no  more  than  10  bits,  and  six  values  can 
be  stored  in  one  60-bit  word.  The  array  space  required 
for  this  simple  example  is  reduced  to  one-sixth  of  the 
usual  unpacked  array  space.  Unfortunately,  the  simplex 
algorithm  does  not  benefit  from  this  technique,  for  most 
arrays  of  the  algorithm  are  real-valued.  The  four  integer 
variables  in  blank  COMMON  seldom  change,  and  although  they 
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may  be  packed  into  one  word,  the  savings  of  three  words 
of  storage  is  hardly  substantial. 

Masking  expressions  with  the  use  of  the  SHIFT  and 
MASK  functions  (both  nonstandard  FORTRAN  features)  are 
beneficia]  for  reducing  execution  time  when  repetitious 
multiplication  and  division  by  powers  of  2 or  the  repeated 

lr 

use  of  the  modulus  function  of  the  form  MOD(n,  2 ) occurs. 
The  potential  for  saving  memory  by  this  method  is  not  evi- 
dent. 

In  a test  of  the  modified  algorithm  on  a small  model, 
the  following  statistics  were  generated: 

36576g  words  alone  to  execute  the  program 

40500g  words  to  load  and  execute  the  program 

14.992  seconds  central  processing  time,  including 
6.487  seconds  of  time  to  compile  the  source 
code 

22.921  seconds  input/output  processing  time 
The  final  solutions  match  those  provided  by  the  original 
algorithm  developed  by  Ignizio.  Coefficients  in  the  orig- 
inal model  are  read  into  TE  (hence,  influencing  input/output 
time)  and  are  calculated  in  the  updated  version  (hence, 
influencing  central  processing  time) . In  the  updated  ver- 
sion, MAIN,  INIT , and  COEF  accomplish  what  the  original 
GMAIN  \ased  to  perform,  yet  the  number  of  words  of  code 
generated  during  compilation  for  the  three  subprograms 
is  greater  than  that  for  GMAIN.  The  two  models  differ 
noticeably  in  design. 


In  order  to  reduce  computer  memory  further,  we  use 
segmentation  loading.  In  some  computer  codes,  the  program 
units  need  not  reside  in  memory  simultaneously.  Two  dis- 
tinct phases  exist  in  our  source  code.  When  initialization 
occurs,  testing  for  entering  and  departing  basic  variables 
is  not  necessary.  Initialization  routines  are  unnecessary 
once  the  execution  phase  begins  operations.  Segmentation 
is  an  efficient  way  of  loading  compiled  subprogram  modules 
into  central  memory  so  that  compiled  program  units  which 
are  not  involved  with  the  current  processing  are  not  loaded 
into  memory;  hence,  those  units  uninvolved  are  not  occupy- 
ing extra  space  in  the  machine.  The  binary  subprograms 
are  arranged  by  segmentation  directive  statements  into 
"tree  structures."  These  subtrees  share  memory. 

Segmentation  is  capable  of  significantly  reducing 
the  amount  of  central  memory  needed  by  large  programs. 
It  is  much  more  powerful,  efficient,  and  flexible 
than  overlay  for  large  core  programs  and  is  therefore 
recommended  for  this  class  of  problem  [Ref  2:5-13]. 

Communication  among  the  segments  is  two-way  unlike 

the  one-way  communication  of  overlay  loading.  Labelled 

COMMON  blocks  are  considered  ?s  well  by  the  segmentation 

directives  GLOBAL.  The  segments  include  a "root  segment" 

which  contains  the  MAIN  program.  Out  of  this  root  segment 

are  branches  containing  other  modules  of  the  particular 

program.  The  segmentation  load  for  the  updated  model  is 

depicted  in  Figure  4-4.  The  directives  that  create  this 

structure  are  also  presented. 
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Figure  4-4.  Segmentation:  Tree  Structure  and  Directives 


For  the  updated  design,  the  segmentation  structure 
includes  two  subtrees  and  a branch  stemming  from  the  root 
segment.  One  subtree  T1  deals  with  the  initialization. 
During  this  phase,  COEF  is  called,  provides  TE  coefficients 
and  resumes  control  to  INIT.  INIT  directs  the  reading 
and  placement  of  objective  function  values  with  PLACE. 

The  split  between  COEF  and  PLACE  is  inherent;  hence,  the 
tree  structure  involves  two  branches  out  of  the  segment 
SI  where  INIT  resides. 

The  execution  phase  is  divided  into  two  logical  parts: 
(1)  a testing/permuting  phase,  and  (2)  an  output  phase. 

This  division  appears  in  the  structure  by  segments  S3  and 
T2 . 

It  is  important  to  understand  that  each  branch  of 
the  structure  (S0-S1-S5,  S0-S1-S6,  S0-S3,  S0-S4-S7,  or 
S0-S4-S8)  is  mutually  exclusive  from  another  branch.  Thus, 
only  one  branch  occupies  central  memory  at  one  time.  The 
blank  COMMON  block  which  is  first  declared  in  MAIN  is  acces- 
sible to  all  segments  for  it  resides  in  the  root  segment. 

\ 

Communication  is  preserved  as  it  is  in  the  normal  loading 
scheme . 

The  same  small  case  above  is  now  run  with  a segmenta- 
tion load  and  reveals  improvement  in  the  use  of  core  memory; 
however,  the  disadvantage  in  this  run  lies  with  increased 
input /output  processing  time. 
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SAVINGS 


35751g  words  alone  to  execute 
34300g  words  to  load  and  execute 


625g  words 
4200g  words 


-.354  sec. 


-7.392  sec. 


15.346  seconds  central  processing 

time  including  5.859  seconds 
to  compile  the  source 
30.313  seconds  input/output  process- 
ing time 

This  run  requires  a minimum  of  35751_  words,  while  the 

O 

normal  load  needs  40500g  words,  a total  savings  of  2527g= 
1367  words  of  central  memory.  Whether  or  not  an  analyst 
employs  segmentation  to  the  updated  model  is  a choice 
involving  the  tradeoff  between  decreased  memory  versus 
increased  input/output  time.  As  the  problem  size  expands, 
the  tradeoff  widens. 


Use  of  the  Simplex  Algorithm 
on  the  Case  of  Three  Skills 
and  Six  Time  Periods 

The  computer  algorithm  is  developed  to  the  point  where 
we  can  solve  the  first  numerical  example  outlined  in  Chapter 
Three.  Initial  conditions  specified  in  that  chapter  are 
read  into  the  program  and  the  algorithm  is  executed  in 
normal  and  segmented  loading  schemes.  Unit  weights  are 
assigned  to  every  element  of  the  objective  function. 

The  results  indicate  that  all  goals  have  been  satis- 
fied except  for  the  fourth  priority  level  goals  (a^  = 
171.758)  and  sixth  priority  level  goals  (a^  = 11.290). 
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We  have  not  fully  met  our  CONUS  requirements  during  some 
of  the  time  periods  and  simultaneously  have  exceeded  CONUS 
requirements  during  other  time  periods.  All  overseas  theater 
requirements  have  been  fully  achieved.  Table  IV-1  summar- 
izes our  resulting  CONUS  and  overseas  levels. 

CONUS  requirements  are  exceeded  slightly  during  the 
first  and  second  time  periods  for  skill  types  1 and  2. 

Our  CONUS  levels  drop  below  CONUS  requirements  during  time 
period  4 for  skill  type  1 and  during  time  periods  3-6  for 
skill  type  2.  CONUS  requirements  are  fully  satisfied  for 
skill  type  3.  At  time  t = 6 only  fifty-three  and  six- 
tenths  percent  (53.6%)  of  skill  type  2 CONUS  requirements 
are  met. 

Obviously,  our  solution  has  resulted  in  a poor  achieve- 
ment of  skill  type  2 (direct  mission  support)  CONUS  goals. 
This  may  seem  unusual  because  we  assigned  no  extra  weights 

in  the  objective  function  in  favor  of  the  other  two  skill 

< 

levels.  However,  our  problem  structure  is  implicitly 
weighted  away  from  the  second  skill  type  for  three  reasons. 
First,  we  desire  that  skill  type  2 levels  increase  by  a 
much  larger  amount  overseas  than  the  amounts  required  for 
the  other  skills.  Second,  the  length  of  technical  training 
for  the  second  skill  type  is  twice  as  long  as  the  length 
of  technical  training  for  the  other  skills,  creating  a 
longer  delay  in  satisfying  skill  type  2 goals.  This  delay 
has  direct  influence  on  goal  quotas  to  be  met  later  in 
the  contingency.  Third,  while  the  desired  rate  of  increase 


Table  IV- 1 


Resulting  CONUS  and  Overseas  Levels 
For  First  Numerical  Example  (in  Thousands) 


t i 

Scenario 

Requirements 

Solution 

Levels 

CONUS 

O/S 

CONUS 

O/S 

1 

1 

75 

29 

75 

29 

2 

152 

42 

157.8 

42 

3 

170 

32 

170 

32 

2 

1 

75 

33 

78.8 

33 

2 

152 

60 

153.7 

60 

3 

170 

35 

170 

35 

3 

1 

75 

50 

75 

50 

2 

152 

83 

143.9 

83 

3 

170 

40 

170 

40 

4 

1 

75 

70 

58.8 

70 

2 

152 

113 

122.1 

113 

3 

170 

45 

170 

45 

5 

1 

75 

85 

75 

85 

2 

152 

138 

104.9 

138 

3 

170 

60 

170 

60 

6 

1 

75 

95 

75 

95 

2 

152 

165 

81.4 

165 

3 

170 

75 

170 

75 

t represents  the  time  period. 


i represents  the 


skill  type. 


for  skill  type  2 overseas  is  much  higher  than  for  the  other 
skills,  the  allocation  of  reserves  among  the  three  skills 
does  not  reflect  this  disproportionately.  Reserves  are 
allocated  equally  among  the  three  skills. 

Table  IV-2  provides  the  various  flows  and  instructor 
levels  resulting  from  the  simplex  solution.  Flows  into 
technical  training  w.  , during  the  time  interval  (0,3) 
reflect  the  peacetime  and  pre-conscription  training  pipe- 
line. Significant  flows  into  skill  type  2 technical  train- 
ing occur  initially  because  of  rapid  skill  type  2 growth 
requirements.  Also,  at  least  during  the  initial  time 
periods,  CONUS  requirements  for  other  skills  are  met  by 
the  reserves  and  the  initial  training  pipeline.  However, 
during  later  time  periods  (t  4)  the  growth  requirements 
for  skills  1 and  3 also  become  more  significant.  Because 
their  technical  training  times  are  shorter,  skills  1 and 
3 begin  receiving  the  majority  of  trainees.  The  shorter 
training  times  enable  quicker  goal  satisfaction.  Moreover, 
beginning  at  t = 5,  gains  into  skill  type  2 technical  train 
ing  have  no  influence  on  meeting  CONUS  and  overseas  require 
ments . 

It  is  also  noteworthy  that  the  third  skill  type 
instructor  pool  appears  to  be  a sink  for  absorbing  person- 
nel who  are  excess  to  CONUS  needs.  The  excesses  are  per- 
sonnel who  are  conscripted  after  CONUS  requirements  are 
met  to  the  maximum  extent.  Recall  that  we  do  impose  the 
requirement  that  all  available  labor  be  utilized  during 
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the  last  two  time  periods.  Most  gains  into  technical  train- 
ing during  the  last  two  time  periods  are  in  skill  type 
3 because  more  skill  type  3 instructors  may  be  obtained 
from  the  CONUS.  While  we  desired  that  CONUS  requirements 
not  be  exceeded,  we  placed  no  upper  limit  on  the  number 
of  instructors. 

The  solution  discussed  is  one  of  several  (25)  alternate 
solutions  resulting  from  the  problem.  Each  optimal  solu- 
tion is  basically  the  same  result  as  presented  above  with 
only  minor  differences  in  certain  decision  variables.  Thus, 
for  each  alternate  solution,  the  objective  function  satis- 
faction levels  are  identical,  and  the  resulting  CONUS  and 
overseas  levels  are  identical.  The  difference  between 
successive  solutions  occurred  with  the  direct  flows  of 
personnel  from  the  CONUS  force  to  the  instructor  pool, 
from  this  pool  to  the  CONUS  force,  and  the  activation  of 
reserve  forces  for  use  in  the  CONUS  force.  In  all  cases 
the  differences  between  solutions  were  compensating  and 
yielded  basically  the  same  result.  For  instance,  one  solu- 
tion differed  from  the  first  printed  solution  in  the  vari- 
ables of  the  second  time  period.  During  that  time  period, 
1200  fewer  persons  left  the  CONUS  to  become  instructors. 
Simultaneously,  2400  fewer  persons  in  the  reserves  were 
activated  coupled  with  1200  more  instructors  leaving  their 
technical  training  positions  for  assignments  to  the  CONUS. 
The  CONUS  force,  hence,  remained  unchanged  from  the 
original  strength  allocation  of  the  first  optimal  solution; 
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the  composition  of  that  force  was  created  differently  by 
each  successive  solution.  The  fact  that  no  upper  limit 
was  placed  on  the  instructor  pool  contributed  to  the  genera- 
tion of  alternate,  but  basically  the  same,  solutions. 

The  segmentation  run  showed  identical  results  as  anti- 
cipated. The  comparison  of  computer  source  utilization 
follows  in  Table  IV-3.  The  creation  of  the  segments  for 
segmentation  loading  contributes  to  the  higher  input/output 
tim . Savings  in  memory  over  increasing  time  is  an  evident 
tradeoff . 


Table  IV-3 

Comparison  of  the  Two  Loads 


Normal  Load 

Segmentation  Load 

Memory  (words) 
to  execute 

107050g 

105654g 

to  load  and  execute 

35700g 

31100g 

Processing  Time  (sec.) 

central 

522.355 

524.791 

input /output 

23.385 

50.233 

Now  we  attempt  to  improve  the  achievement  of  skill 
type  2 CONUS  requirements  by  utilizing  appropriate  weights 
in  the  objective  function.  Specifically,  we  assign  weights 
of  10  to  skill  type  2 CONUS  negative  deviation  variables 
for  3 <_  t <_  6 . Additionally,  we  assign  weights  of  5 
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to  skill  type  1 CONUS  negative  deviations  for  the  same 
time  periods.  We  retain  unit  weights  for  skill  type  3. 


Table  IV-4  provides  a comparison  of  the  solution  using 
these  different  weights  with  the  first  solution.  In  both 

* cases,  the  overseas  goal  requirements  were  fully  met. 

Our  second  solution  provides  sixty-two  and  four-tenths 
percent  (62.4%)  achievement  (as  compared  to  fifty-three 
and  six-tenths  percent  (53.6%)  for  the  first  solution). 

9 

CONUS  requirements  are  underachieved  now  at  t = 4 and 
t = 5 for  skill  type  1 (as  compared  with  underachievement 
for  this  skill  type  only  at  t = 4 in  the  initial  solution) . 
Even  with  unit  weights  assigned  to  skill  type  3 deviation 
variables,  implicit  weighting  is  still  in  favor  of  the 
third  skill  type.  This  circumstance  is  probably  due  to 
the  fact  that  the  attrition  rate,  the  instructor  needs, 
and  the  buildup  requirements  are  at  the  minimum  for  the 
third  skill  type.  Effective  measures  for  reducing  the 
number  of  redundant  alternate  solutions  are  provided  later 
in  the  thesis. 

Further  Algorithm  Development : 

Employment  of  Core-Swapping 

As  we  attempt  to  treat  larger  problems  and  expand 

* the  contingency  to  a situation  involving  three  job  cate- 
gories over  twenty  time  periods,  difficulties  in  size  of 
storage  capacity  result.  For  the  earlier  cases,  adapta- 
tion of  array  dimensions  and  file  buffer  sizes  pose  no 
problems.  The  blank  COMMON  block,  though,  for  the  new 
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Table  IV- 4 


Comparison  of  Solutions:  One  Without  Weights  and 


One 

With 

Different  Weights 

Among  Skill 

Types 

CONUS 

Levels 

t 

i 

Requirements 

Test  1 

Test  2 

1 

1 

75 

75 

75 

2 

152 

157.8 

154.3 

3 

170 

170 

170 

2 

1 

75 

78.8 

76.6 

2 

152 

153.7 

152 

3 

170 

170 

170 

3 

1 

75 

75 

75 

2 

152 

143.9 

141.2 

3 

170 

170 

170 

4 

1 

75 

58.8 

57.7 

2 

152 

122.1 

120.7 

3 

170 

170 

173.2 

5 

1 

75 

75 

60.8 

2 

152 

104.9 

104.2 

3 

170 

170 

170 

6 

1 

75 

75 

75 

2 

152 

81.4 

94.8 

3 

170 

170 

170 

t represents  the  time  period, 
i represents  the  skill  type. 
Test  1 had  no  weights  applied 
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situation  enlarges  significantly.  For  the  smallest  test 
case,  the  COMMON  block  is  a storage  area  of  12706  words; 
overall  field  length  to  load  and  execute  is  16704  words 
which  is  reasonable.  For  the  new  situation,  the  COMMON 
block  resembles  the  following  and  requires  261968  words: 

COMMON  TL (343,6) , TT(6,723),  TE(343,723),  TI(6,723), 
TB (343) , TA(6) , JCOL(723,2),  JROW(343,2), 

XI ( 380 , 2 ) , NOB J , NPRI , NVAR,  NCOL 

The  total  capacity  of  the  CDC  6600  computer  is  131072  words. 
It  is  impossible  to  continue  using  the  current  updated 
model  with  this  excessive  demand  for  storage. 

The  majority  of  the  space  in  blank  COMMON  is  occupied 
by  the  TE  matrix  (247989  words) . This  matrix  is  sparse 
but  the  position  of  coefficients  is  crucial  for  the  pivot- 
ing routines  of  the  method.  Knowing  which  variables  leave 
and  enter  the  basis  during  a given  iteration  at  a priority 
level  must  be  retained.  Dealing  with  the  formidable  size 
of  TE  requires  the  capability  to  separate  TE  into  smaller 
blocks  and  then  use  each  block  where  necessary.  The  program 
must  adapt  the  capability  of  reading  and  writing  new  values 
of  TE  after  each  call  to  the  permuting  subprogram.  Now 
the  question  we  impose  is:  How  should  the  TE  matrix  be 
stored  and  retrieved? 

Three  forms  of  available  input/output  processing  are 
possible  answers:  (1)  mass  storage,  (2)  buffer  input/out- 
put, and  (3)  unformatted  input/output . Mass  storage  is 
beneficial  whenever  a program  requires  that  information 

90 


be  accessed  in  a random  fashion.  This  process  tends  to 
be  expensive  and  prohibitive,  for  the  mass  storage  sub- 
routines require  excessive  input /output  time.  As  we  observe 
the  input/output  requirements  of  the  subprograms  that  access 
TE,  the  data  is  needed  in  a sequential  manner  (from  the 
first  row  of  TE  to  the  final  row  of  TE) , and  mass  storage 
becomes  a less  attractive  option  for  our  solution. 

The  second  alternative  saves  central  memory  allocation 
by  eliminating  the  space  required  of  the  file  buffer.  The 
comparison  of  input/output  processing  time  of  this  approach 
with  that  of  unformatted  input/output  leads  to  the  abandon- 
ment of  buffer  input/output  for  the  solution  of  the  problem. 
The  choice  of  unformatted  input /output  allows  us  to  specify 
variable  lengths  of  computer  records  to  be  processed  as 
well.  Because  each  constraint  set  has  one  of  three  pos- 
sible lengths  of  rows  (JOBCAT,  JTIMES,  or  their  product, 

NEW) , unformatted  processing  permits  variable  length  records 
and  can  process  a specific  constraint  set  more  easily  than 
buffer  input/output.  Not  every  element  of  the  TE  array 
in  blank  COMMON  has  to  be  accessed  with  the  use  of  unformat- 
ted processing  whereas  buffer  input/output  requires  every 
element  of  TE  to  be  processed  at  once.  Storage  of  informa- 
tion is  larger  when  using  buffer  input/output  statements 
than  when  using  unformatted  files. 

The  updated  version  of  separate  phases  is  modified 
to  accept  core-swapping  with  unformatted  input /output . 
(Core-swapping  is  the  technique  of  keeping  an  entire 


portion  of  a matrix  on  file  but  not  in  central  memory  at 
once;  parts  of  the  matrix  are  processed  one  at  a time  so 
that  the  memory  space  allocated  to  the  array  is  small.) 

The  single  main  program  of  the  updated  model  is  divided 
into  three  main  programs.  The  first  main  program  includes 
the  initialization  phase  without  the  establishment  of  the 
arrays  of  initial  objective  function  values.  These  arrays 
are  excluded  from  the  blank  COMMON  block  and  reduce  the 
allocation  of  this  program.  Only  TE  and  the  array  of  right- 
hand-side  values  exist  in  this  program. 

As  the  largest  number  of  rows  in  a constraint  set 
which  are  generated  at  a time  by  COEF  is  NEW,  the  space 
for  TE  in  COMMON  must  have  that  number  as  the  row  dimen- 
sion in  the  COMMON  statement.  The  number  of  necessary 
columns  of  TE  may  be  reduced  further  to  be  either  the  num- 
ber of  decision  variables  (NVAR)  or  the  number  of  goal 
equations  (NOBJ) , whichever  is  larger.  Thus,  TE  is  split 
up  and  shortened  into  blocks  as  large  as  JOBCAT*JTIMES 
rows  and  max (NVAR, NOBJ)  columns,  either  the  lef thand-side 
of  TE  or  the  righthand-side  of  TE.  (See  Figure  4-5.)  For 
the  large  case  under  consideration,  the  TE  block  assignment 
in  COMMON  is  60  rows  by  380  columns  for  a total  of  22800 
words . 

An  additional  feature  allows  for  more  than  these  16 
blocks  of  TE  to  be  generated  during  this  initialization 
program.  The  integer  variable  INDEX  is  read  in  to  allow 
the  program  to  divide  the  generated  block  of  a constraint 


Figure  4-5.  Breakdown  of  Split  Technology  Matrix 

cet  into  smaller  blocks  of  a maximum  of  INDEX  rows  each. 

For  instance,  if  INDEX  = 20  and  we  deal  with  this  "3  by  20" 
case,  the  first  constraint  set  generates  a block  of  60 
equations,  but  the  program  divides  this  block  further  into 
three  blocks  of  20  rows  each  prior  to  writing  the  block 
to  the  unformatted  file.  Other  constraint  sets  follow 
similarly,  and  for  instances  like  the  third  constraint 
set  where  only  three  equations  are  created,  the  block 
generated  in  COEF  is  written  to  the  file  as  it  is.  The 
program  tries  to  split  a generated  block  into  as  many  sub- 
blocks as  possible  with  an  INDEX  number  of  rows.  Hence, 
the  execution  phase  of  the  algorithm  (the  second  and  third 
main  programs)  requires  a TE  matrix  allocated  in  blank 
COMMON  equal  to  INDEX  rows  and  max (NVAR,NOBJ)  columns. 
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Three  new  subroutines — BEGIN,  FINISH,  and  RHSTE — are 
constructed  to  generate  the  blocks  on  the  unformatted  file. 
As  this  phase  operates,  the  entire  lef thand-side  of  TE 
is  written  on  the  file.  Each  generated  subblock  is  a record 
of  the  file.  After  the  eighth  constraint  set  is  complete, 
the  righthand-side  of  the  TE  matrix  is  similarly  divided 
into  associated  blocks  like  those  of  the  lef thand-side 
and  written  on  the  file.  Once  the  TE  matrix  is  processed, 
the  array  of  righthand  goal  values  becomes  the  final  record 
on  the  file. 

Deciding  what  value  of  INDEX  to  select  takes  experi- 
mentation. There  is  no  fixed  rule  for  the  decision.  We 
must  remember  that  the  smaller  the  INDEX,  the  smaller  the 
TE  dimensioning  and  the  more  blocks  processed  on  the  file. 
More  blocks  processed  on  the  file  implies  greater  input/ 
output  processing  time  in  the  execution  phase.  The  larger 
the  INDEX  becomes,  the  greater  amount  of  storage  in  blank 
COMMON  is  required  but  at  a decrease  in  time  for  the  input/ 
output  functions. 

For  the  "3  by  20"  case  with  an  INDEX  of  20,  36  blocks 
of  initial  TE  coefficients  are  created  on  the  file.  The 
matrix  of  righthand-side  goal  equations  coefficients  which 
are  calculated  during  initialization  is  written  as  the 
final  block  of  data  on  the  unformatted  file.  The  first 
main  program  then  requires  12.909  seconds  execution  time 
and  31.660  seconds  input/output  time.  The  field  length 
necessary  to  load  and  execute  the  program  is  112731  words. 


Concerning  the  retrieval  of  TE  values  from  the  file 
and  using  the  values  correctly  with  other  matrix  values, 
the  model  needs  a method  to  calculate  the  correct  sub- 
scripts of  the  original  TE  matrix  and  the  matrices  always 
resident  in  core.  Subscripting  the  variables  grows  more 
complex  since  core-swapping  is  the  way  of  handling  the 
storage  capacity  problem.  (Core-swapping  is  the  technique 
of  processing  information  in  a large  array  without  having 
the  entire  array  in  memory  at  one  time.  This  array  is 
kept  complete  on  an  external  file  on  magnetic  disk  or  tape, 
and  the  records  of  the  file  which  are  subblocks  of  the 
array  are  read  and  processed  one  at  a time.)  Once  the 
program  determines  which  TE  coefficient  is  requested,  it 
must  know  whether  or  not  that  the  current  TE  block  in  cen- 
tral memory  contains  the  requested  value.  If  the  current 
block  is  not  the  "right"  one,  that  block  which  is  "right" 
must  be  read  into  central  memory,  and  the  coefficient  can 
be  retrieved  directly  from  core. 

Concerning  the  updating  of  the  TE  matrix  on  file, 
the  program  generates  new  values  each  time  that  PERM  is 
called.  Once  PERM  ends,  the  file  must  be  correct  so  that 
the  program  can  continue  searching  for  the  optimal  solution. 
It  is  impractical,  though,  to  construct  a function  sub- 
program to  handle  this  retrieval  and  updating  process  for 
only  one  value  at  a time.  Instead,  whole  blocks  are  updated 
at  a time.  The  elements  of  the  TE  matrix  are  not  randomly 
chosen;  subroutine  CINDX  requires  every  value  in  row  after 


row  of  TE  in  order  to  compute  new  index  rows;  subroutine 
PERM  pivots  and  simultaneously  updates  each  TE  element 
in  a row-by-row  fashion.  Instead  of  relying  on  a function 
subprogram  to  retrieve  TE  values  singly,  we  must  develop 
the  structure  in  these  extant  subprograms  to  read  each 
record  of  TE,  modify  or  access  it,  and  continue  sequen- 
tially until  the  final  record  of  TE  is  processed.  The 
structure  includes  reformulating  double  summations,  espe- 
cially where  the  inner  summation  iterates  over  TE  row 
subscripts.  (Sometimes  this  untangling  of  indexes  of  double 
summations  is  impossible  to  reverse,  but  the  process  is 
successful  with  CINDX,  PERM,  TEST,  and  ALTST.) 

To  accomplish  this  development,  a new  array  referred 
to  as  IN  is  introduced  into  COMMON  to  keep  a historical 
set  of  information  on  the  initial  creation  process  of  the 
unformatted  file.  IN  is  dimensioned  eight  rows  by  three 
columns,  where  each  of  these  eight  rows  refers  to  a par- 
ticular constraint  set  of  TE.  The  first  column  indicates 
the  original  number  of  equations  generated  in  the  constraint 
set,  the  second  indicates  the  number  of  rows  contained 
in  the  final  subblock  of  the  set,  and  the  third  indicates 
the  number  of  subblocks  created  for  the  set.  The  computer 
code  is  adapted  through  reordering  of  the  DO  loops  to  handle 
both  the  lefthand-side  of  the  TE  matrix  (first  half  of 
the  file)  and  the  righthand-side  of  the  matrix  (second 
half  of  the  file)  in  that  order.  Each  block  is  read  in 
and  processed.  Correct  limits  on  DO  loops  for  reading 


the  unformatted  file  are  specified  by  integer  variables 
and  change  only  when  the  final  block  of  a constraint  set 
is  about  to  be  read. 

The  program  is  modified  to  handle  the  updating  of 
the  TE  matrix  in  PERM  by  reading  from  one  unformatted  file 
to  a second  one.  The  integer  variable  NTAPE  introduced 
in  blank  COMMON  refers  to  the  file  ready  for  reading  and 
the  local  integer  variable  NSCRAT  refers  to  the  file  ready 
for  writing  the  corrected  TE  records.  Once  PERM  is  near 
the  conclusion  of  permutation,  the  variable  NTAPE  is  set 
to  equal  NSCRAT;  in  this  manner,  the  other  subprograms 
use  the  most  up-to-date  information  of  the  TE  coefficients. 
Input/output  time  is  saved  by  not  having  to  rewind  the 
original  file  and  rewrite  the  new  information  onto  the 
same  file. 

To  maintain  the  pivoting  structure  as  before,  communi- 
cations between  PERM  and  the  subroutines  TEST  and  ALTST 
require  a vector  of  the  coefficients  of  the  TE  column  of 
the  entering  basic  variable.  Both  TEST  and  ALTST  create 
this  vector  prior  to  their  subsequent  calls  to  PERM.  (ALTST 
must  create  this  vector  twice  because  of  the  second  call 
to  PERM  in  its  program  logic.)  PERM  uses  this  vector  and 
creates  a second  vector  of  coefficients  of  the  row  of  the 
departing  basic  variable.  PERM  requires  both  vectors  of 
coefficients  in  permuting  the  TE  matrix  and  prevents  the 
unwanted  alteration  of  these  original  values  before  the 
permutation  process  is  complete.  The  original  values  are 


necessary  to  achieve  permutation.  The  problem  does  not 
exist  in  the  updated  model  where  the  entire  technology 
matrix  resided  in  central  memory. 

Now  that  the  storage  capacity  of  TE  is  reduced,  it 
is  possible  to  reduce  the  required  central  memory  further. 
The  index  rows  are  always  created  by  a call  to  CINDX.  Once 
the  index  rows  are  in  memory,  the  top  stub  values  which 
produce  this  set  of  index  rows  are  not  needed  until  the 
flow  of  control  passes  to  PERM.  We  observe  that  the  top 
stub  values  and  the  index  rows  have  identical  dimensions 
in  COMMON  and  are  mutually  exclusive  matrices;  when  one 
is  required  in  memory  by  a subprogram,  the  other  is  not 
referenced.  Thus,  we  can  eliminate  the  extra  space  declared 
for  one  of  the  matrices,  allow  both  sets  of  values  to  occupy 
the  same  locations  in  core  when  needed,  and  refer  to  the 
mutual  array  as  TT  (formerly  the  reference  to  top  stub 
values) . 

After  initialization  has  positioned  the  objective 
function  weights  in  TT,  it  is  necessary  to  copy  these  values 
to  a file.  The  choice  of  buffer  input/output  on  this  file 
saves  extra  memory  by  eliminating  the  file  buffer.  After 
the  copy  is  complete,  the  integer  variable  NTT  signifies 
which  set  of  values  is  in  memory  at  the  present  time.  This 
switch  indicates  whether  or  not  top  stubs  are  in  memory, 
for  only  these  values  are  read  in  or  written  on  the  file. 

In  subprogram  CINDX  the  index  row  values  evolve  through 
computation  with  the  current  top  stub  values.  The  switch 


is  set  according]  In  subroutine  PERM,  the  too  stub  values 
are  changed  because  of  the  departure  of  and  entrance  of 
variables  in  the  TE  matrix  into  the  basis.  The  file  of 
current  top  stub  values  is  rewound  and  read  into  memory. 
While  in  central  memory,  the  changes  are  enacted,  and  the 
new  information  is  buffered  out  to  the  file,  ready  for 
a subsequent  call  to  CINDX. 

The  second  main  program  of  the  new  algorithm  begins 
with  an  initialization  routine  to  process  the  objective 
function  values  and  to  read  in  variable  labels.  Because 
central  processing  time  increases  as  a result  of  the  pro- 
cessing of  the  large  technology  matrix,  the  main  program 
uses  the  SECOND  function  to  determine  if  a length  of  pro- 
cessing time  has  elapsed.  Once  this  length  of  time  is 
reached,  control  flows  to  a newly  designed  subroutine  STORE 
which  stores  the  current  TE  values  on  one  file  and  all 
remaining  COMMON  information  along  with  the  current  priority 
level  on  a second  file. 

The  third  main  program  of  the  new  algorithm  differs 
from  the  second  main  program  in  a few  respects.  Instead 
of  an  initialization  routine  to  read  initial  values  into 
matrices,  this  third  program  begins  by  restoring  the  arrays 
and  integer  variables  in  COMMON  to  the  values  that  they 
assumed  during  the  run  of  the  second  main  or  in  a previous 
run  of  this  program.  Thus,  subroutine  PLACE  is  deleted 
from  the  third  main  program.  Subroutine  RESTOR  is  an 
alternate  entry  point  for  subroutine  STORE.  When  the 


computations  in  MAIN  dealing  with  elapsed  time  indicate 
that  a predetermined  quantity  (3000  sec)  has  been  reached, 
the  third  main  program  calls  subroutine  STORE  to  prepare 
the  files  with  the  latest  values.  Repetition  of  this  pro- 
gram is  necessary  until  no  new  optimal  entering  variable 
results  for  the  last  priority  level. 

The  final  design  of  the  core-swapping  algorithm  is 
depicted  in  Figure  4-6. 


Running  the  Final  Model 

In  order  to  execute  this  "3  by  20"  contingency  model 
on  the  CDC  6600  computer,  we  are  led  to  use  magnetic  tape 
as  a medium  of  storage.  Two  tapes  are  necessary  for  hold- 
ing the  TE  matrix  and  the  COMMON  block  data. 

During  the  testing  of  the  model,  we  encountered 
several  hardships  in  dealing  with  the  computer — other 
users'  priority  runs,  computer  downtime,  and  unforeseen 
coding  errors.  One  significant  problem  is  the  storage 
of  our  initial  and  subsequent  TE  matrices  on  magnetic  disk. 
The  TE  matrix  for  this  model  requires  nearly  490  blocks 
of  computer  allocation.  Trying  to  make  this  data  file 
permanent  for  later  runs  is  an  impossible  task  when  we 
consider  magnetic  disk  space.  The  data  file  of  TE  elements 
is  large  enough  to  begin  dumping  earlier  permanent  file 
information.  The  problem  does  not  materialize,  however, 
with  the  use  of  magnetic  tape. 
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Figure  4-6.  Final  Model  for  Large-Scale  Contingencies 
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With  the  first  main  program,  one  tape  is  necessary 
for  storing  the  initial  TE  coefficients  and  righthand- 
side  goal  values.  This  local  file  is  referred  to  as  TAPE 1 . 
During  the  second  main  program  the  tape  file  is  called 
TAPE1  again,  and  by  the  time  that  STORE  prepares  to  retain 
the  most  current  TE  data,  the  tape  file  is  mounted  on  the 
tape  drive  during  this  program;  by  the  time  subroutine 
STORE  ends,  the  remaining  COMMON  block  data  is  written 
on  the  second  tape,  and  the  file  is  referred  to  as  TAPE4. 

The  third  main  program  requires  both  tape  files  to 
restore  the  data  and  begin  anew.  During  this  program's 
execution,  TAPE8  is  written  to  TAPE1,  a magnetic  disk  file 
which  is  not  a permanent  file.  TAPE1  and  TAPE2  are  the 
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two  local  files  which  alternate  in  subroutine  PERM  as  the 
tape  with  the  most  current  TE  coefficients. 

Similarly  TAPE4  initializes  the  magnetic  disk  file 
TAPE3  with  TT  values.  TAPE3  and  TAPE4  are  processed  by 
buffer  input/output  commands.  Keeping  TAPE3  as  well  as 
TAPE1  and  TAPE2  as  local  files  on  magnetic  disk  saves  input/ 
output  processing  time  from  having  to  rewind  magnetic  tapes 
and  write  on  the  slower  medium  of  magnetic  tape. 

Fortunately,  the  modified  simplex  algorithm  allows 
the  user  to  obtain  a global  optimum  solution  according 
to  the  implicit  policy  structure.  Infeasible  intermedi- 
ate basic  solutions  are  cause  for  program  termination, 
a safeguard  in  the  algorithm.  Running  the  case  of  three 
job  skills  over  twenty  time  periods  is  a slow  process 
despite  the  conservation  features  implemented  in  the  over- 
haul of  the  updated  model.  The  effects  of  enormous  storage 
capacity  are  felt  in  the  other  computer  resources  of  time. 

A numerical  test  example  is  run  with  an  objective 
function  similar  to  that  in  the  earlier  example  of  3-skills, 

6-time  periods.  Now,  policy  objectives  include  a buildup 
of  the  force  to  double  the  theater  force  size  over  the 
last  eight  months  of  the  conflict.  With  an  INDEX  value 
set  at  20  (allowing  no  more  than  twenty  rows  of  TE  matrix 
to  be  in  central  memory  at  one  time) , the  three  main  pro- 
grams require  a maximum  field  length  of  164046g  words  for 
loading  and  execution.  The  savings  in  computer  memory 
is  offset,  though,  by  several  hours  of  central  processing 

J 


r 


104 


time.  Five  hundred  total  iterations  are  necessary  to  test 
and  permute  the  TE  matrix  for  the  six  priority  levels 
assigned.  Because  of  the  sequential  access  of  the  data 
on  file,  the  permutation  routine  handles  a quarter  of  a 

i 

million  calculations  on  every  iteration.  The  routine  CINDX 
which  computes  index  rows  needs  to  access  the  file  for 
every  iteration  from  the  first  priority  level  to  the  cur- 
rent priority  level.  The  iterations  at  the  sixth  priority 
level  take  longer  to  accomplish  because  of  the  necessity 
to  access  the  data  more  frequently  (from  the  first  to  the 
sixth  priority  level) . 

The  final  solution  seems  cogent  with  roundoff  error 
kept  at  a minimum.  Observing  intermediate  solutions  failed 
to  detect  roundoff  error  beyond  the  eighth  decimal  place 
(a  value  equivalent  to  .00001  man).  The  final  solution, 
though,  shows  some  roundoff  error  in  the  fourth  and  fifth 
decimal  places  but  not  an  alarming  result. 

Because  our  programming  with  simplex  consumes  more 
computer  time  than  permissible  for  accurate  and  adequate 
analysis,  we  decided  to  try  an  alternate  approach,  the 
pattern  search  techniques  for  goal  programming.  These 
techniques  become  the  focus  of  our  next  chapter. 


V.  Development  and  Testing  of  the 
Pattern  Search  Algorithm 


Introduction 

In  this  chapter  we  develop  and  use  an  alternate  solu- 
tion procedure  for  our  model  which  is  derived  from  James  P 

4 

Ignizio's  pattern  search  method  for  nonlinear  goal  program 
k ming  (Ref  15:166-175).  Ignizio's  method  is  based  on  an 

extension  of  the  direct  search  method  of  Hooke  and  Jeeves 
(Ref  14:202-229)  and  is  a fairly  easy  method  to  employ 
or  to  program.  The  method  is  selected  primarily  because 
it  requires  a great  deal  less  computer  memory  and  computa- 
tion time  than  does  the  simplex  procedure  for  large  prob- 
lems . 

The  Pattern  Search  Procedure 

Ignizio's  pattern  search  procedure  consists  of  three 
major  components:  the  Hooke  and  Jeeves  search  algorithm, 

f 

the  system  of  goal  equations  with  righthand-side  values, 

- and  the  objective  function.  We  established  the  general 

* expressions  for  the  last  two  components  earlier  in  this 

research.  Prior  to  initiation  of  the  search  algorithm, 
two  more  items  are  needed.  We  need  an  initial  guess  for 
each  decision  variable  so  that  we  may  establish  a base 
point  from  which  to  begin  our  search.  We  also  need  a step 
size  for  each  decision  variable  to  identify  how  each 
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variable  is  incremented  during  the  search  procedure.  Nor- 
mally, the  better  our  choice  of  an  initial  base  point, 
the  faster  the  search  will  converge  to  an  optimal  solution. 

The  Hooke  and  Jeeves  search  algorithm  begins  at  the 
initial  base  point.  For  our  particular  model  the  base 
point  Xjjggg  consists  of  an  initial  guess  for  each  decision 
variable  cit,  vit,  yljt,  uit,  wl<t,  et  Vi,t). 

Each  priority  level  component  of  the  objective  function 
is  evaluated  at  the  base  point  to  provide  the  initial  satis- 
faction of  goals  at  each  priority  level.  This  action  is 
accomplished  by  the  following  procedures. 

1.  The  goal  equations  are  evaluated  in  terms  of  the 
initial  guess  for  the  decision  variables. 

2.  The  values  of  the  goal  equations  are  compared 
with  the  righthand-side  values,  and  the  deviation  vari- 
ables are  evaluated. 

3.  Each  priority  level  component  of  the  objective 
function  is  evaluated  in  terms  of  the  applicable  deviation 
variables  multiplied  by  the  assigned  weights. 

The  initial  values  of  the  various  priority  level  components 
of  the  objective  function  provide  the  initial  base  criterion 
for  determining  whether  or  not  improvements  are  made  in 
the  objective  function  during  the  search. 

Having  established  the  initial  base  point  and  initial 
evaluation  of  the  achievement  function,  we  now  begin  the 
search  procedure.  The  initial  decision  variable  (in  our 
model,  x.  .)  is  increased  by  an  amount  equal  to  the 


After  this  first  step  we  reevaluate  our  goal  equations 
and  the  objective  function  to  determine  if  an  improvement 
in  the  objective  function  has  occurred.  Examples  of  improve- 
ment include: 

1.  An  improvement  in  the  first  priority  level  goal 
achievement. 

2.  No  change  in  the  achievement  of  first  priority 


level  goals,  but  an  improvement  in  the  achievement  of 
second  priority  level  goals. 


Examples  of  no  improvement  include: 

1.  No  change  in  achievement  of  any  priority  level 
goals . 

2.  Improvement  in  the  achievement  of  second  priority 
level  goals,  but  a worsening  in  the  achievement  of  first 
priority  level  goals. 

Moreover,  our  preemptive  priority  goal  structure  is 
in  effect.  We  attempt  to  satisfy  higher  priority  level 
goals  prior  to  considering  lower  priority  level  goals. 

We  attempt  to  improve  lower  priority  level  goal  achievement 
only  if  we  do  not  worsen  higher  priority  level  goal  achieve- 
ment. 

If  an  improvement  in  the  objective  function  has 
resulted  from  our  step  increase  in  x.  . , then  x..  . will 

1 r 1 if  J- 

retain  its  new  increased  value.  If  no  improvement  has 
resulted,  we  will  then  decrease  the  original  value  of  x^  ^ 
by  an  amount  equal  to  the  step  size  specified  for  x^  ^ 
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(x,  , = x,  , - 6 ) , reevaluate  our  goal  equations  and 

the  objective  function,  and  determine  if  an  improvement 
has  resulted  from  this  action.  If  there  has  been  an  improve- 
ment over  our  initial  evaluation  of  the  objective  function 
(based  on  our  initial  base  point) , then  x.  . will  retain 
its  new  decreased  value.  If  no  improvement  has  resulted, 
x^  ^ will  be  returned  to  its  original  base  value.  The 
value  of  the  objective  function  against  which  further  com- 
parisons of  attempts  for  improvement  will  be  made  is  the 

value  based  on  the  "best"  value  of  x.  . to  this  point  (and 

1/1 

the  remaining  initial  base  values  of  our  other  decision 

variables) . If  our  initial  increase  in  x.  . has  resulted 

in  an  improvement  in  the  objective  function,  then  we  would 

not  test  for  an  improvement  based  on  a decrease  in  the 

initial  value  of  x^  Instead,  we  would  keep  x^  ^ at 

its  new  increased  value  and  test  for  improvements  based 

on  incrementing  our  second  variable  x0  . . 

With  x^  ^ at  its  new  "best"  value  and  the  objective 

function  at  its  value  based  on  x.  . and  the  initial  base 
» J- » J- 

values  of  the  remaining  decision  variables,  we  now  incre- 
ment x2  ^ and  test  for  improvements  in  the  objective  func- 
tion. We  then  hold  x„  . at  its  new  value  (x.  . remains 

Z / 1 if  J- 

at  its  value  found  earlier)  and  increment  the  third  deci- 
sion variable.  This  process  continues  until  all  decision 
variables  have  been  incremented.  Once  all  variables  have 
been  incremented  and  tested  for  improvements,  the  new  values 
of  our  decision  variables  xngw  provide  an  improvement  in 
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our  objective  function  as  compared  to  our  initial  base 

values  x, 

base 

One  method  of  continuing  our  search  would  be  to  let 
xnew  become  our  new  base  point  and  begin  our  step  incre- 
ments of  each  decision  variable  over  again.  However,  we 
use  a method  developed  by  Hooke  and  Jeeves  (Ref  15:168- 
169)  instead  to  provide  acceleration  to  the  search  pro- 
cedure. We  find  a test  point  xtest  which  is  located  in 

the  direction  of  x from  x.  , but  twice  the  distance 

new  base' 

between  x and  x, 

new  base 

xtest  - ^xnew  xbase 


We  then  increment  each  decision  variable  from  its  value 


in  xtest  and  attempt  to  find  improvements  in  the  objective 

function.  We  define  x as  our  base  point  and  call  it 

new 

x.  now.  The  new  values  of  our  decision  variables  as 
base 

modified  from  (and  reflecting  improvements  over)  xfc 

are  now  called  x . We  then  compare  evaluations  of  the 
new  c 

objective  function  based  on  x and  x.  . If  x pro- 
J new  base  new  c 

vides  improvement  in  goal  achievement  over  the  goal 


achievement  provided  by  xbase • we  continue  to  accelerate 

the  search.  (That  is,  x.  . = 2x  _ - x , and  we  then 

test  new  base 

define  x as  our  base  point  and  call  it  x.  „ .)  If, 
new  base 

however,  xnew  provides  no  improvement  in  goal  achievement 
over  the  goal  achievement  provided  by  xj3ase/  we  return 
to  xj3age  and  reduce  our  step  size  for  each  decision 
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variable.  We  then  begin  new,  smaller  increments  of  each 


decision  variable  from  its  value  in  x^gg.  With  the  smaller 
step  sizes  it  is  possible  that  we  may  continue  making 
improvements  through  searching  and  accelerating,  or  we 
may  find  that  we  need  to  reduce  our  step  sizes  even  more. 

The  pattern  search  continues  until  one  of  the  follow- 
ing events  occurs: 

1.  We  reach  an  optimum  achievement  of  goals  after 
which  no  improvements  can  be  made  (within  certain  prescribed 
tolerances  for  each  priority  level) . 

2.  We  terminate  the  search  after  reaching  a pre- 
determined maximum  number  of  search  patterns. 

3.  We  terminate  the  search  after  reaching  a predeter- 
mined maximum  number  of  step  size  reductions. 

4 . We  terminate  the  search  due  to  computer  computa- 
tion time  limitations. 


Minor  Modification  of  Problem  Structure 

Rather  than  continuing  dicussion  of  the  pattern  search 
procedure  at  this  time,  it  is  convenient  to  introduce  a 
minor  modification  in  our  problem  structure.  With  the 
present  goal  equations  and  objective  function  formulation 
we  encounter  difficulties  when  trying  to  support  a rapid 
buildup  of  forces  overseas.  Because  satisfying  overseas 
requirements  is  at  a higher  priority  level  than  satisfying 
CONUS  requirements,  overseas  requirements  are  satisfied 
at  all  costs.  Thus,  it  is  possible  for  CONUS  force  levels 
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to  be  totally  depleted  while  we  attempt  to  satisfy  extremely 
large  overseas  needs.  This  problem  cannot  be  solved  by 
simply  exchanging  the  priority  levels  between  the  CONUS 
and  overseas.  Such  an  action  would  result  in  the  CONUS 
needs  being  fully  met,  but  with  large  underachievements 
of  overseas  goals. 

The  basic  problem  facing  us  is  that  we  want  to  meet 
overseas  requirements,  but  not  seriously  deplete  our  CONUS 
force.  This  problem  is  solved  by  introducing  a new  set 
of  goal  equations  and  a new  preemptive  priority  level. 

Our  ninth  goal  equation  set  is  concerned  with  maintain- 
ing minimum  allowable  CONUS  force  levels.  The  minimum 

levels  represent  a percentage  of  CONUS  requirements  K.  . . 

i , r 

The  equations  are  identical  to  the  second  equation  set 
(People  Available  in  the  CONUS  Base) , except  that  now  the 
CONUS  goal  t is  multiplied  by  the  minimum  fraction  of 
goal  achievement  authorized  ck. 

Ninth  Equation  Set — Minimum  Allowable  CONUS  Levels . 


n (c.  . + u.  . + y.  .-x.  . - v.  .) 

j=l  i,]  i,]  i,j  if]  if] 


+ d9,i,t  d9,i,t 


= a . K . . - K.  . 

l i , t i,0 


(1  - X. ) E . . w.  . . 
l ]=1  i,j-A 


¥i;  t <_  Ai 
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Vi;  t > A±  (5.2) 


We  introduce  a new  preemptive  priority  level  in  our 
objective  function  concerned  with  meeting  minimum  CONUS 
requirements.  A reasonable  place  for  the  new  priority 
level  is  between  our  present  second  and  third  levels.  The 
new  level  becomes  the  third  priority  level  with  each  remain- 
ing level  being  shifted  down  one  unit. 


Old  Priority 
Structure 


New  Priority 
Structure 


a,.  = t W~  . d” 

3 t=l  i=l  9 , l , t 9 , l , t 


(5.3) 


Additionally,  we  make  a minor  change  to  our  second 
priority  level  goals  related  to  meeting  desired  trainee 
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per  instructor  ratios.  We  now  express  the  desire  that 
trainee  per  instructor  ratios  be  met  exactly,  rather  than 
simply  not  exceeded. 


a_  = E*?  , E?  , (W*  . . dt  . . + w“  . . d"  . .)  (5.4) 

J -1-  D ( 1 f U D / 1 / L-  D / 1 / L.  O / 1 f t 


This  change  allows  added  flexibility  because  we  can  return 

to  our  original  second  priority  level  structure  by  setting 

W,  . . = 0 ¥i,t.  However,  unless  otherwise  stated,  all 

o , l , t 

weights  in  our  objective  function  will  be  at  the  unit  level. 
We  now  have  an  effective  method  for  preventing  unnecessary 
growth  of  the  instructor  pool  and  eliminating  the  genera- 
tion of  redundant  alternate  solutions. 


General  Comments  About  Pattern 
Search  and  Our  Model 

The  pattern  search  method  has  been  recommended  by 
Ignizio  as  a useful  approach  for  solving  nonlinear  goal 
programming  problems  (Ref  15:167).  This  is  partly  because 
a simplex  procedure  does  not  exist  for  the  nonlinear  case. 
Ignizio  does  state  that  the  pattern  search  method  cannot 
guarantee  a solution  which  is  a global  optimum  (Ref  15:171). 
Moreover,  the  procedure  may  converge  to  a local  optimum 
^,id  is  dependent  upon  the  initial  guess  for  the  decision 
variables.  Depending  upon  the  problem  structure  and  the 
initial  guess  for  xj-,a~e'  the  solution  procedure  may  exceed 
the  number  of  patterns  allowed  or  exceed  the  available 
computation  time  before  converging  to  a solution. 
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Therefore,  it  is  important  that  we  thoroughly  understand 
our  problem  structure  and  make  a reasonable  first  guess 
for  our  decision  variables.  We  also  need  to  identify 
methods  for  avoiding  the  convergence  of  our  search  to  a 
local  optimum.  Instead,  we  desire  that  our  procedure  con- 
verge to  the  true  global  optimum  solution. 

We  are  attempting  to  use  a nonlinear  procedure  to 
solve  a linear  problem.  Although  we  have  a simplex  pro- 
cedure for  linear  problems,  we  are  unable  to  use  the  simplex 
for  realistic  size  formulations  of  our  problem.  This  is 
due  to  computer  memory  and  computation  time  limitations. 
However,  by  using  the  pattern  search  method,  we  are  able 
to  reduce  considerably  the  amount  of  computer  memory  needed. 
While  the  simplex  procedure  requires  computer  space  for 
a large  simplex  tableau  (with  many  zero  entries) , the  pat- 
tern search  procedure  requires  space  only  for  a small  num- 
ber of  vectors  (such  as  x,  , x , x.  . ) . The  goal 
equations  need  not  be  accounted  for  in  a large  matrix  as 
required  by  the  simplex  procedure.  Instead,  the  values 
for  the  goal  equations  are  generated  into  a vector  Y by 
using  FORTRAN  arithmetic  statements  which  take  advantage 
of  the  recursive  nature  of  the  equations.  Therefore,  com- 
puter space  requirements  are  relatively  small  for  the  pat- 
tern search  method. 

The  computation  time  problem  encountered  while  using 
the  simplex  procedure  can  also  be  alleviated  to  a large 
extent  when  using  the  pattern  search  method.  In  a large 
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size  linear  goal  programming  problem  involving  several 
decision  variables  and  goal  equations,  a great  deal  of 
time  is  expended  in  performing  the  tableau  manipulations 
necessary  to  find  the  global  optimum  solution.  The  pat- 
tern search  method  can  also  expend  a great  deal  of  time 
searching  for  a solution  (and  perhaps  find  only  a local 
optimum) . However,  if  we  are  familiar  with  our  problem 
structure  and  can  make  a "good"  initial  guess  for  our  deci- 
sion variables,  then  we  can  reduce  the  amount  of  computa- 
tion time  needed  considerably. 

For  example,  it  is  reasonable  to  assume  for  our  prob- 
lem that  the  first  three  priority  level  goals  will  be  fully 
satisfied.  We  know  that  the  first  priority  level  goals 
must  be  satisifed.  Nothing  should  prevent  the  achievement 
of  the  second  priority  level  goals  related  to  maintaining 
desired  trainee  per  instructor  ratios.  Likewise,  our  third 
priority  level  goals  concerned  with  keeping  the  CONUS  force 
within  a certain  percentage  of  requirements  should  also 
be  met.  However,  we  may  not  be  able  to  satisfy  the  fourth 
priority  level  goals  (overseas  requirements)  without  also 
violating  the  achievement  of  our  third  priority  level  goals. 
Using  such  reasoning  based  on  the  characteristics  of  our 
problem  structure,  we  can  determine  a reasonable  first 
guess  for  our  decision  variables  which  satisfies  several 
of  the  initial  priority  level  goals.  Hopefully,  this  guess 
will  greatly  reduce  the  search  time  needed  for  the  problem. 
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While  the  simplex  procedure  has  the  disadvantages 
of  requiring  larger  amounts  of  computer  memory  and  computa- 
tion time,  simplex  does  have  one  distinct  advantage  over 
the  pattern  search.  Given  enough  computer  memory  and  com- 
putation time,  simplex  will  find  a global  optimum  solution 
(subject  to  possible  roundoff  error) . However,  the  problem 
of  converging  to  a local  optimum  when  using  the  pattern 
search  can  be  alleviated  somewhat  if  we  do  three  things. 

1.  Develop  a thorough  understanding  of  the  problem 
structure . 

2.  Determine  a "good"  first  guess  for  the  decision 
variables  based  on  the  problem  structure. 

3.  Modify  the  search  procedure  slightly  to  take 
advantage  of  the  problem  structure. 

We  now  discuss  the  implications  of  attempting  to  solve 
our  problem  using  the  pattern  search  procedure  discussed 
thus  far  in  this  chapter.  The  discussion  demonstrates 
the  need  for  accomplishing  the  three  actions  listed  above. 

For  our  problem  structure  the  decision  variables  are 
listed  in  the  following  order,  x = (x.  . , c.  , v.  , y.  . , 

lft  1 f t 1 f t 1 f t 

u.  , , w.  , e , ¥i,t),  and  this  ordering  affects  how  the 

lft  1ft  t 

search  is  accomplished.  Let  us  assume  that  we  decide  to 

set  our  initial  guess  at  the  origin  (xkase  = (0,0,...0)), 

since  this  is  where  we  begin  our  solution  procedure  using 

the  simplex.  We  begin  incrementing  followed  by 

and  so  forth,  until  we  have  incremented  each  x.  . variable 

i,t 

by  one  step  size.  We  then  increment  the  c's,  v's,  y's, 
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u's,  w's,  and  e's  similarly.  With  each  incrementation 
we  try  to  improve  goal  achievement. 

Let  us  assume  that  there  is  a requirement  to  build 
forces  in  the  overseas  theater  (fourth  priority  level) 
and  that  our  initial  CONUS  force  begins  at  our  CONUS  require- 
ments level  (fifth  priority  level) . Furthermore,  we  desire 
that  the  CONUS  force  remain  at  initial  peacetime  levels 
throughout  the  contingency.  Additionally,  we  will  allow 
some  reductions  in  the  CONUS  force  to  within  a certain 
prescribed  percentage  of  our  CONUS  requirements  (third 
priority  level) . With  the  above  problem  structure,  move- 
ments of  personnel  from  the  CONUS  to  the  overseas  theater 

will  result.  (x.  . values  will  take  positive  steps  from 
1 , t 

zero.)  A small  step  increase  in  ^ from  xbase  = 

(0,0,...,0)  provides  no  change  in  the  achievement  of  first, 
second,  or  third  priority  level  goals.  But,  there  is  an 
improvement  in  achievement  of  the  fourth  priority  level 


goals . 

After  the  initial  incrementing  of  the  x.  variables, 

1 f t 

we  then  increment  the  c.  . variables.  A positive  step 

1 , t 

for  c.  . negates  the  overseas  gains  made  by  x.  , so  this 
x , t i , t 

event  is  not  allowed.  A negative  step  for  c.  from  zero 

1 f t 

results  in  a negative  value  for  c.  , so  this  action  is 

l , t 

not  permitted  either.  Initially,  if  the  x.  variables 

i ,t 

assume  positive  values,  the  c.  variables  will  remain 

l , t 

at  zero. 
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The  next  variables  to  be  incremented  from  zero  values 
are  the  v.  variables.  There  are  no  improvements  in  the 
achievement  of  any  goals  to  be  gained  from  allowing  people 
to  flow  into  the  instructor  pool  from  the  CONUS.  The 
instructor  pool  begins  at  its  initial  strength  1^  and 
no  students  so  far  are  flowing  into  technical  training 
(w.  = 0)  ¥i,t) . Moreover,  removing  people  from  the 

CONUS  worsens  the  achievement  of  CONUS  goals  (fifth  pri- 
ority level)  at  this  point.  The  v.  , variables  will  ini- 

1 9 t 

tiallv  remain  at  zero. 

Each  y.  variable  takes  one  positive  step  so  long 

l , t 

as  improvements  result  in  the  fifth  priority  level  goal 
achievement  and  limitations  on  the  availability  of  reserves 
(first  priority  level)  are  not  exceeded. 

Each  u^  t variable  takes  a positive  step  within  the 
limitations  we  have  set  on  maintaining  a stable  instructor 
force  (first  priority  level)  so  long  as  improvements  are 
made  in  meeting  CONUS  requirements  (fifth  priority  level). 
If  the  y.  values  have  already  ensured  the  satisfaction 

1 r t 

of  CONUS  requirements,  then  the  u.  variables  remain  at 

l , t 

zero. 

The  next  variables  to  be  incremented  from  zero  values 

are  the  w.  variables.  Each  w.  variable  for  t < t + 1 
l , t i , t 

can  be  incremented  positively  as  long  as  we  stay  within 
earlier  gains  into  basic  training  before  the  contingency 
began  (first  priority  level)  and  also  do  not  worsen  the 
meeting  of  desired  trainee  per  instructor  ratios  (second 
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priority  level).  However,  for  t _>  t + 2 , no  further  posi- 
tive incrementations  for  the  w.  variables  can  result. 

1 1 1 

Such  actions  result  in  violations  of  first  priority  level 

goals  relating  to  earlier  gains  into  basic  training.  We 

have  not  yet  incremented  the  efc  variable  from  zero  values. 

When  e^  = 0,  then  w^  must  also  be  zero. 

The  final  variables  to  be  incremented  from  zero  values 

are  the  efc  variables.  We  cannot  increment  any  e^  variable 

without  violating  first  priority  level  goals.  A positive 

value  for  e^  implies  that  we  also  need  a positive  value 

for  w.  for  at  least  one  skill  level  i.  However,  we 

i,t+2 

found  above  that  w.  .=0  Vt  > T + 2;  i = 1,  ...,n. 

r,t  — 

In  our  second  search  pattern  (after  we  have  accelerated 
to  xtest) i we  notice  for  the  most  part  the  same  variables 
increasing  and  the  same  variables  unable  to  change.  We 
need  the  capability  of  increasing  the  efc  and  w^  t variables 
in  order  to  find  an  optimum  solution.  This  is  because  the 
labor  pool  is  normally  the  greatest  source  of  personnel 
resources  available  to  help  satisfy  overseas  and  CONUS 
requirements.  However,  our  problem  structure  and  initial 
guess  prevent  us  from  taking  advantage  of  our  labor  pool. 
Thus,  any  convergence  is  to  a local  optimum  solution  only. 

In  summary,  we  are  unable  to  converge  to  a global 
optimum  solution  using  the  present  pattern  search  pro- 
cedure for  the  following  reasons; 
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1.  Our  initial  guess  at  the  origin  does  not  take 


' 

advantage  of  the  problem  structure  in  that  we  cannot  use 
the  available  civilian  labor. 

2.  The  search  procedure  allows  incrementation  of 

] 

only  one  decision  variable  at  a time,  although  simultane- 
ous variable  incrementations  might  be  useful. 

3.  The  sequence  of  our  search  is  restricted  to  the 
ordering  of  our  decision  variables. 

Modification  of  Search  Procedure 
for  Our  Problem 

In  order  that  we  may  find  an  improved  (and  hopefully 
global  optimal)  solution  to  our  problem,  we  need  to  modify 
both  our  initial  guess  for  x]3ase  and  our  search  procedure. 

We  want  to  allow  the  realization  of  gains  into  basic  train-  ■ 

ing  from  the  labor  pool.  While  several  possibilities  exist 
for  improvements  in  our  solution  procedure,  we  focus  on 

one  possibility  below,  which  appears  to  be  useful. 

— 

We  first  discuss  an  initial  guess  for  x^ase  which 
supports  a change  in  the  search  procedure  to  be  presented 
later.  In  our  initial  guess  we  now  assign  values  to  the 
x.  t and  c.  . decision  variables  so  that  our  overseas  force 

1 , u l , t 

requirements  are  met  exactly.  These  values  may  be  calcu- 
lated directly  from  the  first  equation  set  (People  Avail- 
able in  the  Overseas  Theater) . As  a result,  the  fourth 
and  sixth  priority  level  goals  are  totally  satisfied.  We 
assign  values  to  the  y.  variables  so  that  the  reserves 

1 / t 

are  activated  into  the  CONUS  force  as  rapidly  as  possible 
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within  the  reserve  activation  limitations.  The  first  pri- 
ority level  goals  relating  to  reserve  activation  are  still 
fully  satisfied. 

We  also  make  a change  in  the  first  priority  level 
goals  relating  to  the  flow  rate  of  people  into  technical 
training  from  basic  training.  We  now  include  only  the 
positive  deviation  variables  from  the  fifth  equation  set 
in  the  first  priority  level  of  the  objective  function. 

This  is  equivalent  to  stating  that  our  gains  into  tech- 
nical training  can  go  below,  but  not  above,  the  amounts 
of  people  who  successfully  complete  basic  training.  This 
is  not  realistic,  because  people  who  successfully  complete 
basic  training  must  enter  one  of  the  technical  schools. 
However,  the  reason  for  this  approach  becomes  apparent 
when  the  change  in  the  pattern  search  method  is  discussed 
later . 

As  part  of  our  initial  guess  for  xj3ase»  we  assign 
the  maximum  possible  values  to  the  et  variables  within 
the  limitations  of  our  civilian  labor  pool.  The  first 
priority  level  goals  relating  to  civilian  labor  limita- 
tions and  flow  rate  of  people  into  technical  training 
remain  fully  satisfied.  Additionally,  we  assign  values 
to  w.  . Vi  so  that  the  trainee  per  instructor  ratios  are 
totally  satisfied  (w^  ^ Q)  and  gains  into  technical 

training  equal  the  basic  training  completion  rate 

(Ei=l  wi , 1 = 6s-t) ’ 
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For  reasons  which  become  apparent  later,  we  assign 
initial  values  to  the  u.  , variables  so  that  all  initial 

1 9 -1* 

instructors  I . n are  returned  to  the  CONUS  force  during 
the  time  interval  (0,1).  We  also  let  p = 1 so  that  one 
hundred  percent  (100%)  of  our  instructors  may  be  returned 
to  the  CONUS  any  time.  We  set  all  remaining  decision  vari- 
ables equal  to  zero. 

Before  discussing  changes  in  the  search  procedure, 
it  is  useful  to  review  our  current  extent  of  goal  satisfac- 
tion based  on  our  new  initial  guess  for  xbase* 

First  Priority  Level.  All  first  priority  level  goals 
are  fully  satisfied. 

Second  Priority  Level.  All  second  priority  level 
goals  (trainee  per  instructor  ratios)  are  fully 
satisfied. 

Third  Priority  Level.  Because  we  now  meet  our  over- 
seas goals  exactly,  we  probably  do  not  have  our 
CONUS  force  at  or  above  its  predetermined  per- 
centage of  CONUS  requirements.  Therefore,  third 
priority  level  goals  are  not  likely  to  be  satisfied. 

Fourth  Priority  Level.  Overseas  force  requirements 
are  completely  satisfied. 

Fifth  Priority  Level.  CONUS  force  requirements  are 
probably  not  satisfied. 

Sixth  Priority  Level.  Overseas  requirements  are  not 
exceeded,  so  sixth  priority  level  goals  are  com- 


pletely satisfied. 
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Seventh  Priority  Level.  Because  we  have  activated 
the  reserves  as  quickly  as  possible,  our  CONUS 
force  probably  exceeds  CONUS  requirements  in  the 
early  periods  of  the  contingency.  (The  opposite 
case  is  likely  to  be  true  in  later  periods.) 

Therefore,  seventh  priority  level  goals  are  not 
likely  to  be  fully  satisfied. 

Our  solution  procedure  is  now  modified  in  order  to 
help  avoid  the  problems  we  encountered  earlier.  Key  aspects 
of  the  modification  include  selectively  deciding  which 
decision  variables  to  iterate  and  determining  useful  com- 
binations of  variables  to  iterate  simultaneously. 

First  Step.  We  now  begin  our  search  procedure  by 
"freezing"  all  decision  variables  at  their  initial  base  . p 

values  in  x,  , except  for  the  u.  , and  w.  , variables 
base  i,t  i,t 

for  t > 2 and  the  v.  . variables  for  t > 1. 

- i,t  - 

Second  Step.  We  now  begin  the  simultaneous  incremen- 
tations of  v^  w^  2'  anc^  ui  2'  W1  2 '*'s  ;*-ncrementec*  by 
its  normal  step  size,  but  v.  . and  u.  „ are  incremented 

1 f X X f z 

by  the  step  size  of  w.  _ divided  by  the  desired  trainee 
per  instructor  ratio  for  skill  level  1,  E^.  All  instructors 
were  initially  returned  to  the  CONUS  during  the  time  inter- 
val (0,1)  with  our  initial  guess  for  u^  Therefore, 

in  order  to  bring  any  students  into  technical  training 
during  the  time  interval  (1,2)  using  w.  ,,  we  need  to  have 
a sufficient  number  of  instructors  available  at  t = 1. 

This  is  accomplished  by  moving  just  enough  instructors 
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into  the  instructor  pool  using  v..  . to  satisfy  our  desired 
trainee  per  instructor  ratio.  Additionally,  since  we  are 
not  yet  incrementing  w.  , from  0,  we  also  return  all 
instructors  to  the  CONUS  during  the  time  interval  (1,2) 
using  u^  This  avoids  worsening  the  achievement  of  second 

priority  level  goals.  While  our  second  priority  level 
goals  remain  fully  satisfied,  we  realize  improvements  in 
our  third  and  fifth  priority  levels  (CONUS  requirements) 
as  a result  of  simultaneously  incrementing  v.  1 , w.  „, 
and  u^  2* 

After  incrementing  v^  w^  ^ , and  u^  ^ • we  then 
simultaneously  increment  w ^ 2'  anc^  u2  2‘  Thi-S  is 

followed  by  the  combinations  v.  . , w.  u.  _ for 

1 / 1 1 f £ 1 / z 

i = 3,...,n.  We  then  accomplish  a second  cycle  of  incre- 
mentations for  v.  . , w.  0,  u.  _ ¥i,  a third  cycle,  and 

so  forth.  We  continue  this  process  as  long  as  we  continue 
to  make  improvements  in  achieving  CONUS  goals  within  the 
limitations  imposed  by  the  labor  pool  and  training  pipe- 
line (first  priority  level) . Thus  far  in  our  modified 
search  procedure,  we  have  allowed  only  the  decision  vari- 
able combinations  v.  .,  w.  „,  u.  _ Vi  to  be  incremented. 

i,l  i,2  i,2 

One  way  to  interpret  the  limited  search  procedure  is  to 
say  that  we  are  trying  to  satisfy  CONUS  goals  as  much  as 
possible  for  the  remainder  of  the  contingency  by  the  flow 
of  trainees  into  technical  training  in  the  time  interval 
(1,2)  . 
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After  maximum  goal  satisfaction  has  resulted  from 
the  incrementations  of  v.  , , w.  _,  and  u.  0 ¥i,  we  then 
freeze  these  variables  in  their  present  values.  We  then 
perform  a similar  simultaneous  search  procedure  for 

v.  „,  w.  , , and  u.  _ Vi,  followed  by  a search  for  v.  _, 

1/4  1/  J 1 f J 

w.  . , and  u.  . Vi,  and  so  forth.  Our  last  simultaneous 

1 f 4 1/4 

search  is  made  for  v.  , . , w.  , and  u.  , . Since  we  earlier 

i / n“l  i / h / l / n 

set  our  values  for  efc  Vt  to  their  maximum  possible  levels 
based  on  civilian  labor  limitations,  we  encounter  no  prob- 
lems initially  in  incrementing  our  w^  values.  Since 
we  no  longer  require  that  our  gains  into  technical  training 
at  least  equal  our  completion  rate  for  basic  training  (first 
priority  level) , then  we  use  the  numbers  of  trainees  enter- 
ing technical  training  needed  to  provide  improvements  in 
our  CONUS  goals. 

Third  Step.  After  maximum  goal  satisfaction  has 
resulted  from  the  incrementations  of  the  combinations  of 
v-  . , w.  , u.  as  discussed  above,  we  now  "freeze" 

1 / t If  t'  1 1/  1 

all  decision  variables  at  their  latest  values.  At  this 
point  we  should  notice  that  significant  improvements  have 
been  made  in  meeting  CONUS  requirements  (third  and  fifth 
priority  levels) . The  first,  second,  fourth,  and  sixth 
priority  level  goals  are  still  fully  satisfied.  Some 
worsening  in  the  seventh  priority  level  has  probably 
resulted  because  we  have  further  exceeded  CONUS  require- 
ments during  the  initial  phases  of  the  contingency.  Even 
though  we  have  made  improvements  in  meeting  CONUS 
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requirements,  we  probably  still  have  not  fully  met  our 
prescribed  minimum  CONUS  force  levels  (third  priority  level) 
throughout  the  contingency.  There  is  an  apparent  conflict 
here,  because  at  the  same  time  all  overseas  theater  force 
requirements  have  been  met  (fourth  priority  level) . 

Fourth  Step.  We  now  attempt  to  rectify  the  conflict 
between  the  third  and  fourth  priority  level  goals.  Our 
approach  is  to  postpone  the  movement  of  people  to  the  over- 
seas theater  from  the  CONUS  until  our  CONUS  force  is  at 
least  at  its  minimum  authorized  levels. 

The  problem  facing  us  is  represented  by  Figure  5-1. 
There  are  no  goal  conflicts  in  the  time  intervals  (0,  t^) 
and  (t2,20)  for  the  20-period  contingency  represented. 
However,  during  the  time  interval  (t^,t2)  our  CONUS  force 
has  dropped  below  its  minimum  authorized  levels  for  skill 
level  i because  we  have  satisfied  overseas  requirements 
instead.  To  rectify  this  conflict,  we  need  to  insure  that 
our  actual  CONUS  levels  never  drop  below  the  minimum  auth- 
orized CONUS  levels. 

A method  for  solving  this  problem  is  to  institute 
a search  procedure  which  moves  the  actual  CONUS  force  levels 
to  the  minimum  authorized  levels  during  the  interval 
(ti , t2) . As  a result  we  still  have  maximum  achievement 
of  overseas  requirements  within  the  limitations  imposed 
by  the  third  priority  level. 

We  now  begin  our  search  by  simultaneously  incrementing 
i by  a negative  step  and  x^  2 by  a positive  step.  This 


* 


Type  i 


Figure  5-1.  Conflict  Between  Achievement  of  Third 
and  Fourth  Priority  Level  Goals 


procedure  has  a net  effect  of  postponing  the  flow  of  per- 
sonnel overseas  while  we  attempt  to  improve  third  priority 
goals.  Normally,  no  improvements  will  result  for  the 
x.  . , x.  0 combination.  We  attempt  the  same  procedure 
for  the  remaining  combinations  x.  .,  x.  „ for  i = 2,...,n. 

1/1  1 t 

When  no  improvement  occurs,  the  variables  are  returned 
to  their  original  values.  We  next  attempt  the  same  pro- 
cedure for  the  x.  _,  x.  _ combinations.  Again,  no  improve- 

1 / c.  1 / 

ments  probably  result.  Let  us  assume  that  after  time  t^, 
the  CONUS  force  for  skill  level  1 begins  to  drop  below 


minimum  CONUS  requirements.  The  simultaneous  negative 
incrementation  of  x^  and  positive  incrementation  of 

X1  t +2  result  in  an  improvement  in  the  third  pri- 

ority level  goals  for  time  t^+1.  At  the  same  time  no 
changes  will  occur  in  the  third  priority  level  for  other 
time  periods  for  skill  level  1.  A similar  point  in  time 
is  also  reached  for  the  other  skill  types  after  which 
improvements  are  made  in  the  third  priority  level  goals 
by  the  search  procedure.  There  may  also  be  a time  later 
in  the  contingency  t 2 when  further  improvements  in  the 
third  priority  goals  are  not  required  because  minimum  CONUS 
force  levels  are  met. 

The  search  procedure  described  above  continues  until 
each  variable  combination  x.  . , x.  . has  been  allowed 
to  increment  by  a single  step.  The  net  result  to  this 
point  will  be  an  improvement  in  the  achievement  of  minimum 
CONUS  requirements  and  a worsening  in  achievement  of  over- 
seas requirements.  Referring  to  Figure  5-1,  the  actual 
CONUS  level  for  skill  type  i has  increased  by  a step  size 
for  t^<t<t2»  We  have  begun  the  movement  of  our  actual 

I 

CONUS  levels  to  the  minimum  desired  CONUS  requirements 
for  t^<  t<  t2  • 

We  now  continue  our  search  process  by  repeating  the 
procedure  described  above  until  the  third  priority  level 
goals  are  fully  satisfied.  For  each  skill  type  and  time 
period  in  the  interval  (t^,t2>,  we  attempt  to  just  satisfy 

I 
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third  priority  level  goals,  so  that  we  still  meet  over- 
seas goals  to  the  maximum  extent. 

Fifth  Step.  We  now  attempt  an  improvement  in  the 
achievement  of  seventh  priority  level  goals  relating  to 
our  desire  not  to  exceed  CONUS  requirements.  For  those 
initial  time  periods  in  which  the  CONUS  force  is  over- 
strength, we  simply  reduce  the  values  of  the  y.  vari- 

i > t 

ables  accordingly  and  allow  the  differences  in  reserves 
to  be  activated  later  when  needed. 

Sixth  Step.  We  now  correct  the  problem  of  certain 
gains  into  basic  training  being  in  excess  of  the  needs 
of  our  contingency.  Rather  than  use  the  maximum  values 
of  efc  speficied  by  our  initial  guess,  we  instead  calculate 
the  levels  of  efc  needed  to  provide  the  personnel  flows 
into  technical  training  found  in  our  earlier  search  pro- 


cedure. 

Seventh  Step.  During  the  second  step  of  our  modified 
search  procedure,  we  allowed  frequent  movements  of 
instructors  back  and  forth  between  the  CONUS  and  instructor 
pools.  We  did  this  so  that  desired  trainee  per  instructor 
ratios  (second  priority  level)  were  met  exactly.  While 
such  action  was  convenient  for  the  search  procedure,  we 
do  not  want  this  resulting  instructor  instability  reflected 
in  our  final  solution.  Thus  we  calculate  net  flows  and 


assign  all  net  flows  to  CONUS  as  u.  variables  and  all 

1 / t 

net  flows  to  the  instructor  pool  as  v^  variables. 


Our  problem  is  now  solved  by  our  modified  search  pro- 
cedure. While  there  is  no  guarantee  that  we  have  found 
a global  optimum  solution,  we  have  found  a solution  that 
is  both  useful  and  reasonable.  Our  result  is  much  better 
than  the  results  found  by  the  earlier  pattern  search  method 
because  we  have  taken  advantage  of  the  problem  structure. 


Pattern  Search  Computer  Algorithm 

Appendix  C provides  a listing  of  the  pattern  search 
computer  algorithm.  The  algorithm  is  derived  from  a pat- 
tern search  computer  program  provided  by  Ignizio  (Ref  15: 
228,231-234,242-247)  and  incorporates  the  modified  search 
procedures  discussed  in  this  chapter. 

The  design  of  the  pattern  search  algorithm  is  depicted 
in  Figure  5-2,  and  the  various  subprograms  and  their  related 
functions  are  described  below: 

PATSH  controls  the  initialization  and  execution 

phases  of  the  algorithm  (the  MAIN  program) . 

DATAIN  reads  in  the  initial  data  and  prepares  the 
matrices  for  later  computation. 

HJALG  begins  the  pattern  search  and  accelerates 

toward  a solution  if  achievement  of  goals  is 
possible . 

VALUE  evaluates  the  objective  function  for  improve- 
ment . 

YVALUE  creates  current  values  of  the  goal  equations 
using  the  latest  decision  variables. 


DECIDE  makes  the  actual  decision  whether  or  not  an 


improved  solution  results  from  the  last  call 
to  DECIDE. 

FINISH  stops  the  algorithm  if  the  maximum  number  of 
pattern  search  cycles  is  reached. 

Testing  of  the  Modified 
Search  Procedure 

We  now  conduct  testing  of  the  modified  search  pro- 
cedure using  the  computer  algorithm  for  numerical  problem 
formulations . 

Test  One.  Our  first  test  is  based  on  a three-skill 
type,  twenty-time  period  model.  The  problem  formulation 
includes  the  additional  set  of  goal  equations  and  modifica- 
tion of  the  objective  function  established  in  this  chapter. 

Tables  V-l  through  V-4  summarize  parameters  and  ini- 
tial conditions  for  the  problem.  The  three  skill  types 
considered  are  mission,  direct  mission  support,  and 
indirect  mission  support.  One  time  period  equals  one  month. 

Weights  of  100  are  assigned  to  the  negative  deviation 
variables  concerned  with  meeting  CONUS  requirements  for 
skill  type  2 for  6 <_  t <_  20  (fifth  priority  level)  . Weights 
of  10  are  assigned  to  negative  deviation  variables  con- 
cerned with  meeting  CONUS  requirements  for  skill  type  1 
for  6 <_  t <_  20.  It  is  anticipated  that  the  required  con- 
scription rate  should  be  fairly  low  near  the  end  of  the 
contingency.  Therefore,  zero  weights  are  assigned  to  the 
negative  deviation  variables  concerned  with  using  all 
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Table  V-l 


Parameters  and  Initial  Conditions 


Skill  Type,  i 

1 

2 

3 

Initial  CONUS  Level,  K.  _ 

(thousands)  1 ' 

75 

152 

170 

Minimum  Fraction  of  CONUS 

Goal  Achievement  Authorized,  o ^ 

.8 

. 8 

.8 

Length  of  Technical  Training, 

4 

6 

2 

(months) 

Technical  Training  Attrition 

Rate,  X ^ 

0 

0 

0 

Desired  Trainee  per  Instructor 

Ratio,  E. 

l 

8 

5 

10 

Initial  Instructor  Levels, 

(thousands)  1^  _ ^ 

.209 

.634 

. 333 

Zi , 0 

.209 

.634 

. 333 

Initial  Instructor  Flow  Rates, 

(thousands)  v^  ^ 

0 

0 

0 

U • n 
i,0 

0 

0 

0 

Initial  Reserve  Level, 

(thousands)  R.  n 

1 , u 

30 

30 

30 

Desired  Residual  Reserve, 

(thousands)  R^ 

0 

0 

0 

Reserves  Activation  Limitation, 

(thousands)  ^ 

10 

10 

10 

' Yi  t 2 

20 

20 

20 

30 

30 

30 

Initial  Overseas  Levels 

(thousands)  T.  _ 

1 , u 

25 

38 

30 

1 


’ 

I 


I 

i 


Overseas  Attrition  Rates 


a . . 1<  t<  3 


03 

02 

01 


.02 

.013 

.005 


02 

01 

001 


Table  V-2 


Initial  Basic  Training  Pipeline 
(in  Thousands) 


CN 

1 

<U 

e-l 

CD 

O 

11.67 

11.67 

11.67 

Table  V-3 

Initial  Technical  Training  Pipeline 
(in  Thousands) 


i 

w . c 
i,-5 

w . . 
1,-4 

w.  _ 
1,-3 

w . „ 

1,-2 

Wi,-1 

o 

•H 

> 

1 

1.67 

1.67 

1.67 

1.67 

2 

3.17 

3.17 

3.17 

3.17 

3.17 

3.17 

3 

3.33 

3.33 

Table  V-4 

Additional  Parameters  and  Initial  Conditions 


Length  of  Basic  Training,  x (months)  2 

Basic  Training  Survival  Rate,  6 .7 

Maximum  Fraction  of  Instructors 

Permitted  to  Return  to  CONUS,  p 1 

Civilian  Labor  Available 
(thousands) 

Lx  11.67 


wr  » 

available  labor  in  the  latter  time  periods  of  the  contin- 
gency (first  priority  level) . Unit  weights  are  assigned 
to  all  other  deviation  variables  in  the  objective  function. 

Tables  V-5  through  V-7  summarize  results  of  the  modi- 
fied pattern  search  solution  procedure  for  Test  One.  The 
first,  second,  and  third  priority  level  goals  are  fully 
satisfied.  The  sixth  and  seventh  priority  level  values 
are  close  enough  to  zero  so  that  we  may  consider  the  associ- 
ated goals  also  effectively  satisfied.  However,  neither 
I t 

overseas  requirements  nor  CONUS  requirements  are  fully 
met. 

Table  V-5 

Extent  of  Goal  Achievement  for  Test  One 

0 
0 
0 

1526.08 
45336.92 
.40 


1.75 


Table 


G O 

0 U) 

•H  «H 

-P  0 

G > 

r H 0 W 

OGD 

cn  2 

o 

u 


G P 
0 -H  CO 
O G D 
W IT'S 
0 O 
« CJ 


-P  0 
G > 

•H  0 in 
OPD 
W 2 

o 

L> 


co  in 

hHO 


o co  a\ 
vo  cm  co 


lo  co  in 

(N  H O 
HCN  H 


m (\  o 

md 


a>  o 

00  .H  rH 


iH  CM  VO 
VO  CM  CO 


O 00  O 
CO  CM 
H CM  i — I 


in  cm  o 

MO 


VO  CO  CM 
(T*  CM  H 
rH  rH 


rH  CO  VO 
VO  CM  CO 


in  CO  O 
in  vo  co 

rH  CM  rH 


uo  cm  o 
Mn 


oo  co 

O rH 


rH  CO  VO 
VO  CM  co 


O CM  rH 

co 

rH  CM  rH 


in  cm  o 
r-  in 


in  cr>  vo 

HVOH 


rH  CM  VO 
VO  CM  CO 


CO  VO  ^ 

co 

rH  CM  rH 


in  cm  o 


vo  f-H  00 

CM  O'  H 


CM  CM  VO 
VO  CM  CO 


CO  VO  ^ 
t r-*  co 

rH  CM  rH 


in  cm  o 

MAh* 


r-  cm  vo 

CO  rH  CM 
rH  CM  rH 


rH  co  r- 

VO  CM  CO 


CO  VO  ^ 

r-  co 

rH  CM  rH 


in  cm  o 
Mn 


0 H W 
O G D 
cn  trss 
0 o 
« a 


CT>  CM  CM 
CM  M1  CO 


in  CM  rH 

Mn 


00  CM  CM 
CM  ^ CO 


uo  cm  o 
Mn 


co  o in 
co  vo  co 


m cm  o 
r-  m 


coo  in 
co  vo  co 


uo  cm  o 
r^m  r*» 


o co  o 

IT  CD  ^ 


in  co  o 

r-  r^- 


o co  o 
in  oo  h* 


in  cm  o 

r-  m r^* 


o co  uo 
o ^ 


rH  CM  O 
vo  cm 


o co  in 
r*  h h* 


in  cm  o 
r*  lo  r- 


rH  VO  O 

o vo 


rH  rH  O 
VO  CM  l 


U0  00  O 
00  CO  VO 


UO  CM  O 

r-  uo 


co  uo 
o r- 


H CO 
VO  CM  VO 


U0  CO  U0 

o>  vo 


UO  CM  O 

uo 


vo  o o 
r-  rH  cr» 


O CM  H 

vo  co  in 


o oo  o 

H CO  O' 


UO  CM  O 

r"  uo  Is" 


L3 


F 


i 

> 


x> 

(0 

Eh 


U 

O 

Mh 

w 

1 — I 
Q) 
> 
<11 
P) 

)H 

O 

-P 

U 

3 

P 

-P 

c/i 

c 


TS 

c 

(0 

01 

s 

o 


t7> 

c 

•H 

■p 

rH 

3 

01 

<D 


-P 

•* 

•H 


vo 

o 

o 

O 

o 

o 

w 

-P 

• 

• 

0 

• 

• 

• 

T3 

CD 

tH 

o 

o 

o 

o 

o 

C 

rH 

in 

in 

in 

in 

in 

«3 

W 

3 

o 

CN 

CN 

CN 

£ 

r- 

r- 

co 

r- 

r* 

r- 

in 

m 

in  in 

Eh 

■P 

vo 

i — i 

co 

in 

o 

vo 

in 

o vo 

in 

O 

vo 

r-  o 

CN 

r-  o CN 

C 

•rl 

iH 

co 

CO 

CN 

in 

CN 

in 

CN 

in 

O CN 

CN 

O CN  CN 

•H 

tH  CN 

rH  CN 

CD 

C 

O 

CO 

co 

co 

r- 

in 

in 

in 

■P 

■p 

i — 1 

VO 

rH 

vo 

rH 

vo 

<N 

CN 

^ CN 

W 

* 

co 

o 

o 

co 

o 

o 

CO 

o o 

CO 

CN 

co 

CN 

CO  N1  CN 

a) 

•rH 

0 

• 

• 

• 

• 

• 

* 

• • 

• 

• 

• 

0 0 

• 

• • • 

Eh 

> 

1 — 1 

rH 

rH 

tH 

rH 

•H 

CO 

r* 

co 

CO 

H1  in 

in 

H1 

in 

•p 

tH 

vo 

rH 

vo 

tH 

VO 

^ CN 

CN 

CN 

CO 

o 

o 

CO 

o 

o 

CO 

o 

o 

CO  ^ CN 

CO 

CN 

CO 

CN 

H 

tH 

rH 

rH 

tH  ^ 

rH 

rH 

a*  ^ oo 
o co  co 
cm  vo  oo 


co  co 

rH  VO  tH  VO 

co  o O CO  o o 


co 

rH  VO 

co  o o 


m 

^ CN 

CO  (N 


rH  ^ 


ooo  ooo  ooo  ooo 


in  vo 
r-  r*  vo 


r-  rr 

oo  oo  vo 


<y> 

OWN 


Tj*  ^ CN 


rH  CN  CO 


^ 00  CO 


r-  in 

tH  CN 


o 
• • 
rH  m 


ooo 


r-  in 
m o 


H CN  CO  rH  CN  CO 


CN  CN 


rH  CN  CO 


rH  CO  in 


rH  CN  CO 


CN 


CO 


in 


in 

^ CN 
CO  T^r  CN 


tH  H1 


cr> 

o r- 

tH  VO  CTl 

G\ 

00 

in 

•p 

iH  CJ>  O 

CN  rH  O 

CO  CO  VO 

CO  O 00 

o o co 

ooo 

* 

• • 

• • 

• • • 

• • 

• 

• 

•rH 

CO  rH 

co  in 

VO  CN  tH 

G\ 

•H 

rH  tH 

tH 

ooo 


in  <y>  vo 
00  vo  o 


CN 


in 


H CN  CO 


vO 


139 


> I 





Table  V-7 — Continued 
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represents  the  time  period;  i represents  the  skill  type. 


From  Table  V-6  it  can  be  seen  that,  although  CONUS 
force  requirements  are  not  fully  achieved  in  most  cases, 
the  CONUS  levels  never  fall  below  eighty  percent  (80%) 
of  CONUS  goals  (third  priority  level) . Also,  CONUS  levels 
exceed  CONUS  requirements  only  during  the  first  two  time 
periods  for  the  third  skill  type.  Moreover,  these  excesses 
are  relatively  small. 

Because  the  CONUS  levels  must  meet  or  exceed  eighty 
percent  (80%)  of  CONUS  requirements,  overseas  theater 
requirements  are  not  fully  achieved  for  a large  portion 
of  the  contingency.  While  overseas  goals  are  effectively 
met  during  the  first  three  and  last  four  time  periods, 
significant  goal  underachievements  occur  during  the  interim 
time  periods.  For  example,  at  time  t = 12  the  following 
percentages  of  overseas  goal  achievement  result. 

Skill  Type  1 - Sixty-seven  percent  (67%) 

Skill  Type  2 - Sixth-one  and  three-tenths  percent 
(61.3%) 

Skill  Type  3 - Eighty-six  and  six-tenths  percent 

(86.6%) 


Skill  type  3 (indirect  mission  support)  has  the  best 
overall  achievement  of  overseas  goals,  while  skill  type  2 
(direct  mission  support)  has  the  poorest  overall  achieve- 
ment. This  result  has  occurred  even  though  higher  weights 
were  assigned  to  skill  type  2 objective  function  terms 
and  even  though  larger  relative  step  sizes  were  assigned 
to  skill  type  2 decision  variables  during  the  pattern 
search.  (During  the  pattern  search,  step  sizes  of  1, 
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, respec- 


.5,  and  .1  were  assigned  to  w,  , w„  , and  w 

1 / L A f t 

tively.)  The  poor  results  for  skill  type  2 can  be  explained 
by  at  least  four  factors: 

1.  Skill  type  2 has  the  largest  overseas  growth 
requirements . 

2.  Skill  type  2 has  the  largest  technical  training 
time  with  the  smallest  trainee  per  instructor  ratio. 

3.  The  initial  reserve  force  mix  does  not  reflect 
the  larger  overseas  growth  requirements  for  skill  type  2. 

4.  The  initial  CONUS  force  mix  also  does  not  reflect 
skill  type  2 overseas  growth  needs. 

The  various  flows  shown  in  Table  V-7  support  the 
results  found  for  the  CONUS  and  overseas  forces.  Reserves 
are  activated  on  an  "as  needed"  basis,  and  gains  into  tech- 
nical training  reflect  the  relative  step  sizes  used  during 
the  search  procedure.  Maximum  utilization  of  available: 
labor  occurs  during  the  first  half  of  the  contingency  in 
an  effort  to  build  a wartime  training  pipeline.  But,  dur- 
ing the  second  half  of  the  contingency,  little  use  of 
civilian  labor  occurs.  It  seems  reasonable  to  expect  that 
significant  gains  into  basic  training  should  occur  during 
the  time  period  t = 11,  because  of  the  poor  overall  achieve- 
ment of  skill  type  2 overseas  goals  during  the  middle  time 
periods.  However,  such  gains  do  not  arrive  in  the  CONUS 
force  until  time  period  t = 20. 

Because  we  have  not  specified  that  all  available  labor 
be  utilized  during  the  latter  time  periods  (t  _>  16) , our 


solution  provides  for  only  achievement  of  requirements 
during  the  span  of  the  contingency.  No  provisions  exist 
for  anticipating  requirements  after  the  twentieth  period. 

In  order  to  properly  anticipate  training  pipeline  require- 
ments to  support  time  periods  after  the  twentieth  period, 
we  would  have  to  add  additional  structure  to  our  model. 

This  has  not  been  accomplished  in  the  interest  of  simplic- 
ity. An  alternative  approach  for  handling  this  problem 
would  be  to  test  the  model  for  a figure  like  30  time  periods, 
but  to  use  results  for  only  the  first  20  time  periods. 

Thus,  the  training  pipeline  and  instructor  pool  for  the 
latter  time  periods  of  our  20-period  contingency  would 
be  filled  to  levels  required  to  support  future  needs  during 
the  immediate  time  after  the  twentieth  period. 

The  values  for  the  u.  , and  v.  . variables  result 

1 , t l , t 

from  the  second  step  of  the  modified  search  procedure, 
and,  thus,  reflect  highly  unstable  movements  of  the 
instructor  force.  This  instability  may  be  corrected  by 
calculating  and  using  the  net  flows  between  the  CONUS  and 
instructor  pool.  The  frequent  back  and  forth  flows  are 
only  a solution  technique  for  satisfying  second  priority 
level  goals. 

The  present  solution  requires  55300g  words  of  central 
memory  storage,  8.776  CP  seconds  compilation  time,  191.127 
CP  seconds  execution  time  and  22.609  seconds  input/output 
time. 
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Test  Two . The  model  formulated  for  the  first  test 


is  now  evaluated  under  four  additional  scenarios.  The 
overseas  theater  and  CONUS  requirements  remain  the  same 
for  each  scenario,  but  different  initial  force  mixes  are 
now  examined.  Table  V-8  provides  specific  conditions  used 
for  each  scenario.  Parameters,  initial  conditions  and 
goals  not  specifically  addressed  remain  unchanged  from 
Test  One. 

Test  2A  differs  from  Test  One  in  that  lower  achieve- 
ments of  CONUS  requirements  are  now  authorized  for  skills  1 
and  2 (o^  = a ^ = . 7) . 

Test  2B  differs  from  Test  2A  only  in  the  initial  train- 
ing pipeline.  While  the  total  number  of  people  initially 
in  technical  training  is  unchanged,  larger  numbers  are 
now  in  skills  1 and  2 training.  The  new  initial  basic 
training  level  supports  the  technical  training  levels. 

Test  2C  differs  from  Test  2B  only  in  the  initial 
reserve  force  mix.  While  the  total  initial  reserve  force 
size  is  unchanged,  larger  portions  of  reserves  are  now 
in  skills  1 and  2. 

Test  2D  differs  from  Test  2C  only  in  the  initial  CONUS 
force  mix.  While  the  total  initial  CONUS  force  size  is 
unchanged,  larger  numbers  of  people  are  now  in  skills  1 
and  2.  Initially  the  skill  3 CONUS  level  is  at  eighty- 
two  and  four-tenths  percent  (82.4%)  of  CONUS  contingency 
requirements . 
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Table  V-8 


Testing  Scenarios 


Scenario 

Test  1 

Test  2 

A 

Test  2 

B 

Test  2 

C 

Test  2 
D 

Initial  BMT  Level 
(thousands) 

1(2  classes) 

23.34 

23.34 

18.10 

18.10 

18.10 

Initial  Technical 
Training  Level 
(thousands) 

Skill  1,  E(4  classes) 

6.68 

6.68 

8.68 

8.68 

8.68 

Skill  2,  E(6  classes) 

19.02 

19.02 

23.02 

23.02 

23.02 

Skill  3,  E(2  classes) 

6.66 

6.66 

.66 

.66 

.66 

Initial  Reserve 
(thousands) 

Skill  1 

30 

30 

30 

40 

40 

Skill  2 

30 

30 

30 

45 

45 

Skill  3 

30 

30 

30 

5 

5 

Initial  CONUS  Level 
(thousands) 

Skill  1 

75 

75 

75 

75 

85 

Skill  2 

152 

152 

152 

152 

172 

Skill  3 

170 

170 

170 

170 

140 

Initial  Overseas 

Level  (thousands) 

Skill  1 

25 

25 

25 

25 

25 

Skill  2 

38 

38 

38 

38 

38 

Skill  3 

30 

30 

30 

30 

30 

Minimum  Permitted 
CONUS  Levels 
(fraction  of  CONUS 
goals) 

Skill  1 

.8 

.7 

.7 

.7 

.7 

Skill  2 

.8 

.7 

.7 

.7 

.7 

Skill  3 

.8 

.8 

.8 

. 8 

.8 

— 


Table  V-9  provides  testing  results  under  each  scenario 
for  the  CONUS  and  overseas  theater.  (Even  numbered  periods 
are  shown.) 

It  is  reasonable  that  we  find  improvements  in  overseas 
goal  achievement  for  skills  1 and  2 in  Test  2A  because 
lower  CONUS  levels  for  those  skills  are  now  authorized. 

In  Test  2B  there  are  only  moderate  improvements  in 
overseas  goal  achievement  for  skills  1 and  2 over  Test 
2A.  There  is  also  a somewhat  poorer  achievement  of  skill 
3 overseas  requirements.  Redistribution  of  the  initial 
technical  training  pipeline  has  little  influence  on  improv- 
ing the  achievement  of  overseas  goals. 

Test  2C  provides  significant  overseas  improvements 
for  skills  1 and  2 over  Test  2B.  The  composition  of  the 
initial  reserve  force  plays  a key  role  in  determining  the 
outcome  of  our  contingency.  There  is  also  a poorer  achieve 
ment  of  skill  3 overseas  requirements  than  found  for  Test 
2B. 

Similarly,  Test  2D  provides  even  better  achievement 
of  overseas  needs  for  skills  1 and  2 than  Test  2C.  Thus, 
the  composition  of  the  initial  CONUS  force  is  also  an  influ 
ential  factor  for  the  contingency.  The  achievement  of 
skill  3 overseas  requirements  continues  to  worsen. 

Table  V-10  summarizes  the  percentage  achievement  of 
force  level  requirements  in  periods  8 and  12. 
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Five  Scenario  Testing  Results  (in  Thousands) 
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Table  V-9 — Continued 
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rx 


Test  2C  appears  to  provide  the  best  overall  balance 
of  CONUS  and  overseas  forces  for  the  three  skill  types. 
Better  achievements  for  skills  1 and  2 result  in  Test  2D 
only  at  the  added  expense  of  quite  low  skill  3 overseas 
levels.  In  Test  2D  some  improvements  in  skill  3 overseas 
levels  can  occur  if  minimum  CONUS  levels  for  skill  3 are 
lowered.  However,  any  reductions  in  minimum  CONUS  levels 
affect  our  ability  to  respond  in  a second  theater,  if 
required. 

Moreover,  the  pattern  search  procedure  appears  to 
be  a reasonable  and  useful  solution  algorithm  for  our  model. 
It  is  quite  capable  of  evaluating  a series  of  different 
scenarios  and  demonstrating  the  advantages  and  disadvantages 
of  each.  Because  it  requires  relatively  small  amounts 
of  computer  memory  and  computation  time,  the  pattern  search 
method  should  be  useful  for  solving  problems  involving 
several  job  categories  and  time  periods. 
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VI.  Conclusions  and  Recommendations 

Goal  programming  is  a valuable  analytic  tool  for 
handling  a model  with  multidimensional  goals.  The  objec- 
tives underlying  this  research  have  been  to  broaden  our 
understanding  of  the  art  and  science  of  goal  programming 
and  to  develop  a successful  goal  programming  structure 
for  our  contingency  model.  These  objectives  have  certainly 
been  achieved.  We  hope  that  the  work  which  is  initiated 
with  this  study  advances  the  usefulness  and  development 
of  goal  programming  solutions  to  such  problems. 

The  two  approaches  to  problem  solution  in  our  research 
have  strengths  as  well  as  weaknesses.  Generally,  a weak- 
ness of  one  method  is  a strength  of  the  other.  The  weak- 
nesses and  strengths  that  we  have  encountered  and  considered 
relate  to  the  usefulness,  the  conservation  of  available 
computer  resources,  and  the  applicability  of  the  technique 
to  the  problem. 

The  simplex  algorithm  for  solving  linear  goal  program- 
ming problems  is  independent  of  the  problem  structure; 
its  general  nature  poses  few  restrictions  on  the  type  of 
problem  to  which  it  can  be  applied.  A fundamental  strength 
for  relying  upon  the  simplex  method  is  the  provision  of 
a "best"  possible  solution  to  a formal  policy;  furthermore, 
that  solution  is  guaranteed  to  be  "best" — the  global 
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optimum  over  all  optimal  solutions.  The  algorithm  examines 
every  variable  which  is  not  in  the  current  basis  and  con- 
siders it  for  possible  entry  into  the  basis  at  every  itera- 
tion of  every  priority  level. 

The  principal  weakness  of  the  simplex  method  is  its 
use  and  accommodation  of  computer  storage.  The  method 
requires  large-scale  matrices  for  many  job  categories  and 
many  time  periods  of  interest.  The  position  of  coefficients 
within  the  technology  matrix  is  important  for  the  iterative 
procedure  which  is  the  essence  of  the  method.  As  the  num- 
ber of  job  categories  (n)  or  time  periods  (h)  increases 
by  one,  the  number  of  words  of  central  memory  required 
for  the  technology  matrix  increases  by  a polynomial  of 
degree  3 in  terms  of  both  n and  h.  (With  n job  categories 
and  h time  periods,  the  number  of  decision  variables  equals 
6nh  + h,  and  the  number  of  constraint  equations  equals 
5nh  + 2h  + n.  The  total  number  of  words  required  is  the 
products  of  (5nh  + 2h  + n)  and  (6nh  + h + 5nh  + 2h  + n) 
or  55n2h2  + 37nh2  + 16n2h  + 5nh  + 6h2  + n2.  If  a job  cate- 
gory is  added  to  this  model,  the  increase  in  words  is 

llOnh2  + 92h2  + 32nh  + 21h  + 2n  + 1.  If  an  extra  time 

2 

period  is  included,  the  increase  in  words  is  llOnh  + 

2 

71n  + 74nh  + 42n  + 12h  + 6.)  The  sparse  nature  of  this 

technology  matrix  requires  a minimum  of  forty-seven  per- 
cent (47%)  of  all  TE  coefficients  to  be  zero.  Because 
of  the  storage  requirements,  the  program  has  lengthy  execu- 
tion time  and  processing  time  for  input  and  output. 


Another  weakness  of  the  simplex  algorithm  is  the 
inability  to  adapt  to  nonlinear  constraints.  The  simplex 
method  restricts  all  goal  equations  to  be  strictly  linear; 
although  our  research  is  a development  of  linear  equa- 
tions, the  adaptation  of  nonlinear  equations  to  our  contin- 
gency model  is  not  permitted  with  the  simplex  method.  Even 
logarithmic  or  other  special  transformations  of  the  goal 
equations  are  meaningless  in  explaining  goal  deviations. 

The  pattern  search  technique  offers  several  strengths. 
This  technique  yields  a reasonable  solution  especially 
when  the  computer  algorithm  is  tailored  to  fit  the  goal  struc- 
ture for  the  particular  problem.  The  savings  in  the  use 
of  computer  storage  is  evident,  simply  by  the  fact  that 
zero  values  need  not  be  kept  for  the  goal  equations.  Com- 
putation time  as  well  as  input/output  processing  time  is 
reduced  for  handling  large  cases.  The  "3  by  20"  numerical 
example  presented  in  the  previous  chapter  finishes  in  less 
than  200  seconds  of  central  processing  time  whereas  the 
simplex  version  finishes  in  hours. 

The  weaknesses  in  pattern  search  stem  from  the  flexi- 
bility to  be  tailored  to  a formal  policy  structure.  The 
advanced  search  technique  is  totally  dependent  on  the  par- 
ticular policy.  This  accommodation  of  policy  to  the  design 
of  the  pattern  search  subroutines  requires  the  analyst 
to  examine  every  possible  flow  of  personnel  resources 
thoroughly.  Without  careful  investigation  of  policy  objec- 
tives and  these  flows  of  variables,  the  analyst  discovers 
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that  the  computer  program  iterates  and  searches  blindly, 
seeking  improvement  without  guarantee  of  a mathematical 
global  optimum.  One  final  disadvantage  of  this  approach 
is  the  sensitivity  of  the  algorithm  to  incremental  changes 
in  step  size.  The  computer  program  tries  to  find  the  mini- 
mal value  of  the  objective  function  and  depends  heavily 
on  the  physical  direction  of  the  search  in  Euclidean 
n-space. 

This  research  represents  a first  step  in  the  develop- 
ment and  implementation  of  a goal  programming  model  for 
contingency  planning.  As  a result,  the  model  structure 
has  remained  relatively  simple.  Future  enhancements  can 
be  made  if  one  considers  planning  in  other  resource  areas 
(for  instance,  materiel  resources  in  lieu  of  personnel 
resources) . 

At  this  stage  the  pattern  search  procedure  appears 
to  provide  results  of  an  immediate  and  useful  nature.  Cer- 
tainly, gains  can  be  realized  from  future  research  involv- 
ing this  method.  Moreover,  even  greater  potential  gains 
are  possible  from  enhancement  of  the  goal  programming  sim- 
plex method.  By  incorporating  the  notions  of  revised  sim- 
plex procedures  and  a decomposition  structure  (Ref  29:224- 
285) , it  may  be  possible  to  realize  significant  improve- 
ments in  computational  efficiency  while  also  finding  a 
global  optimum  solution.  Although  Taha  speaks  to  the  linear 
programming  context  for  revised  simplex  and  decomposition 
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appendix  a 

THE  RECRUITING  REQUIREMENTS  MODEL 
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The  model  attempts  to  meet  manpower  requirements  by 
job  category  for  a particular  organization  over  several 
time  periods  and  requires  knowledge  of  the  initial  popula- 
tion per  job  category  as  well  as  total  salary  budget  and 
total  manpower  capability.  Provisions  within  the  model 
account  for  hiring,  firing,  attrition,  and  individual  move- 
ment among  job  categories.  The  model  is  expressed  in  words 
as  shown  in  Figure  A-l  (Ref  7:1-10). 

The  variables  used  in  the  mathematical  formulation 
of  the  model  are  explained  below  (Ref  21:5). 

Ek+(t),  Ek  (t)  - positive  or  negative  deviation, 

respectively,  for  the  k*"*1  manpower 
category  in  time  t 

x^(t)  - manpower  on-board  (in  place)  in  the 

th 

x manpower  category  in  period  t 

y^(t)  - hires  in  the  kfc^  manpower  category 

in  period  t 

z^(t)  - fires  or  reduction-in-force  (RIFs) 

in  the  kfc^  manpower  category  in 
period  t 

G^(t)  - manpower  requirement  (goal)  for  the 

kfc^  manpower  category  in  period  t 


Figure  A-l.  Recruiting  Requirements  Model 
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- the  weight  which  is  applicable  to  the 
positive  goal  discrepancy  for  the  k 
manpower  category  in  period  t 

- the  weight  which  is  applicable  to  the 
negative  goal  discrepancy  for  the  k 
manpower  category  in  period  t 

- the  weight  on  the  hires  in  the  k 
manpower  category  in  period  t 

- the  weight  on  the  fires  or  RIFs  in  the 

th 

k manpower  category  in  period  t 

- Markov-like  transition  matrix  depict- 
ing the  probabilities  of  internal 
personnel  movement  (nu  ^ is  the  frac- 
tion of  personnel  in  job  category  j 
who  move  to  category  i after  one 
period) 

- vector  of  initial  inventory  of  per- 
sonnel in  all  job  categories 

T 

- sum  vector  (i.e.,  e = [1,1, ...1]) 

- manpower  ceiling  for  period  t 

- average  salary  vector  in  period  t 

- salary  budgetary  ceiling  for  period  t 


The  model  may  now  be  expressed  mathematically  (Ref 
21:6).  (The  objective  is  to  minimize  the  weighted  sum  of 
deviations  from  gross  manpower  goals  and  weighted  number 
of  hires  and  fires  (RIFs).) 
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Minimize 


Subject  to 


t [\tE* 


<tl  + !ktV,w  + ^ktykt + {ktzk,t) 


(Manpower  Transition  Conditions) 


x (0)  = a 


- [Mx  ( t-1 ) ] , + x,  (t)  - y,  (t)  + z,  (t)  = 0 


(Manpower  Ceiling  Constraints) 
(Salary  Budget  Constraints) 
(Nonnegativity  Constraints) 


e x ( t)  <_  C(t) 
sT(t)  x ( t ) < B ( t) 


■] 


(Goal  Constraints)  x^(t)  - (t)  + E^  (t)  * G^(t) 


x,(t)  , yv(t)  » z,  (t)  , E,  (t)  , E (t)  _>  0 


The  objectives  of  the  model  are  accomplished  by  using 
relative  weighting  factors.  Weights  are  set  so  that  a 
penalty  is  paid  whenever  the  manpower  requirement  (goal 
for  each  job  category)  is  not  met.  A penalty  is  also  paid 
for  hiring  or  firing.  The  weighting  factors  may  be  adjusted 
so  that  on-board  manpower  and  new  hires  are  favored  over 
excess  personnel  and  RIFs . A given  set  of  manpower  require- 
ments need  only  be  met  as  closely  as  possible,  and  an 
increasing  penalty  (size  of  the  objective  function)  is 
paid  as  one  moves  away  from  the  goals.  Other  planning 
policies  may  be  examined  by  changing  the  weighting  scheme 
(Ref  7:1-9) . 

This  model  differs  somewhat  from  the  example  of  a 
typical  goal  programming  problem  presented  earlier  in  that 
different  preemptive  priority  levels  are  not  shown. 

Instead,  only  relative  weighting  of  goals  is  provided. 


If  weights  used  differ  by  large  enough  amounts,  then  solu- 
tions similar  to  those  obtained  by  using  preemptive  priori- 
ties will  result.  However,  in  general,  the  use  of  pre- 
emptive priorities  and  weighting  is  preferred  to  the  use 
of  weighting  alone  in  order  to  insure  that  higher  priority 
goals  are  attained  first  and  that  satisfaction  of  lower 
goals  will  not  violate  already  attained  higher  goals.  Devi- 
ation variables  are  not  provided  in  the  manpower  ceiling 
constraints  and  salary  budget  constraints.  In  the  manpower 
transition  condition  equation,  the  hires  and  fires  variables 
may  be  considered  as  deviation  variables.  Since  there 
is  only  one  priority  level,  the  problem  can  be  solved  with 
linear  programming. 
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APPENDIX  B 


THE  SIMPLEX  COMPUTER  ALGORITHM 
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Appendix  B.  The  Simplex  Computer  Algorithm 


The  computer  programs  that  follow  are  the  updated 
version  of  the  original  design  (before  the  use  of  core- 
swapping) and  the  successor  to  that  version  in  which  core- 
swapping of  the  technology  matrix  is  handled.  These  pro- 
grams were  designed  to  run  on  the  CDC  6600  computer  system 
associated  with  the  Aeronautical  Systems  Division,  Wright- 
Patterson  Air  Force  Base,  Ohio.  The  version  of  FORTRAN 
Extended,  Version  IV  is  applicable  to  that  machine.  Thus, 
the  user  of  these  codes  ought  to  be  aware  of  the  nonstandard 
FORTRAN  features  employed: 

1.  The  use  of  the  alternate  return  on  subroutine 
CALLS.  This  feature  allows  the  program  to  return  control 
to  a statement  label  in  the  calling  program  unit  if  condi- 
tions within  the  subroutine  warrant  an  alternate  return 
rather  than  a normal  return  to  the  statement  following 
the  CALL . 

2.  The  use  of  the  STOP  "cc...c"  statement  where 
cc...c  is  an  alphabetic  string  of  1-70  characters. 

3.  The  use  of  list-directed  read  and  print  statements 
in  which  no  FORMAT  is  referenced. 

4.  The  use  of  more  than  one  entry  point  to  a sub- 
routine with  the  ENTRY  statement  (notice  ENTRY  PERM  in 
SUBROUTINE  TEST  and  ENTRY  ERRORS  in  SUBROUTINE  POUT) . 


5.  The  use  of  BUFFER  IN  and  BUFFER  OUT  statements 


in  the  core-swapping  algorithms  for  the  top  stub  values. 

6.  The  use  of  multiple  equal  signs  for  one  statement 
as  N = NBLK  = 0.  Otherwise,  this  statement  would  be  written 
as  two : 

NBLK  = 0 
N = NBLK 

7.  The  use  of  the  SECOND  and  EOF  functions  in  order 
to  determine  the  current  system  time  and  to  test  for  an 
end-of-file  mark  on  a file,  respectively. 

The  first  main  program  deals  with  the  three- job-skill , 
six-time-period  case  described  in  Chapter  Four.  The  COMMON 
block  dimensions  should  be  enlarged  accordingly  for  prob- 
lems dealing  with  more  skills  or  more  time  periods.  The 
remaining  three  main  programs  utilize  the  core-swapping 
algorithm  for  a three- job-skill , twenty-time-period  model 
and  should  be  run  in  the  order  presented  in  the  appendix. 
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PROGRAM  M AIN( I SPJT= 10 2B, OUTPUT =70 08,TAPE5=30 08/60, TAPE6«0UTPUT) 
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IMAX  IS  THE  MAXIMUM  NUMBER  OF  TABLEAU  GENERATIONS 
ALLOWED  BEFORE  A RECURSIVE  LOOP  IS  ASSUMED.  NPRT 
IS  THE  NUMBER  OF  PRINTOUTS  POSSESSED. 
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DO  11  J=1,JTIMES 
DO  11  I=i,JORCAT 


THIS  STRUCTURE  OF  TWO  DO  .OOPS  WITH  THE  SAME 
TERMINAL  STATEMENT  PREPARES  THE  TE  COEFFICIENT 
IN  THEIR  LOWER  TRIANGULAR  rDRM  FOR  THE  X*S  AND  C'S 
J=1,NEW  IMPLIES  THE  OUTER  DO  LOOP  WORKS  ON 


ON1;  COLUMN  AT  A TIME 


V < T » J)  IS  THE  rL  OW  RATE  FROM  CDMJS  TO  INSTRUCTORS  IN 
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James  Anthony  Moreno  was  born  on  12  April  1947  in 
Mount  Clemens,  Michigan.  He  graduated  from  high  school 
in  Norton,  Ohio,  in  1965  and  attended  the  University  of 
Akron  from  which  he  graduated  in  1970  with  a Bachelor's 
degree  in  Mechanical  Engineering  (magna  cum  laude)  and 
a regular  commission  in  the  United  States  Army.  He  was 
initially  assigned  to  the  Second  Armored  Division,  Fort 
Hood,  Texas,  as  an  infantry  platoon  leader.  This  was  fol- 
lowed by  assignments  as  operations  officer  and  company 
commander  in  the  15th  Ordnance  Battalion,  60th  Ordnance 
Group  in  Germany  during  the  period  1972  through  1976.  He 
completed  the  Ordnance  Officer  Advanced  Course  prior  to 
entry  into  the  Air  Force  Institute  of  Technology  in  1977. 
He  is  married  to  the  former  Cheryl  Lynn  Schilling  of 
Barberton,  Ohio.  They  have  two  daughters,  Kimberly  and 
Andrea,  and  a son,  James. 

Permanent  address:  676  E.  Paige  Avenue 

Barberton,  Ohio  44203 


Bradley  Wayne  Utz  was  born  21  July  1951  in  Hanover, 
Pennslyvania . He  graduated  from  high  school  in  Hanover, 
Pennsylvania,  in  1969,  and  attended  Lehigh  University  where 
he  received  the  degree  of  Bachelor  of  Arts  in  Mathematics 
(summa  cum  laude)  with  interdepartmental  honors  in  1973. 

291 


I 


While  at  Lehigh,  he  received  a reserve  commission  in  the 
United  States  Air  Force  in  May  1973.  On  his  initial  active 
duty,  he  was  assigned  to  the  20th  Surveillance  Squadron 
(ADCOM)  to  work  as  a space  systems  analyst  and  trainer. 

He  received  his  regular  commission  in  1977  and  entered 
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